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STRESS  CORROSION  OF  CERAMIC  MATERIALS 


SUMMARY 

This  program  consists  of  three  portions:  (1)  Investigation  into  the 
chemistry  of  stress  enhanced  crack  growth  in  glasses  and  ceramics;  (2)  A  study 
into  the  basic  mechanisms  which  govern  crack  growth  and  fracture  in  piezo¬ 
electric  materials,  and  (3)  A  study  into  lifetime  prediction  procedures  for 
ceramic  capacitors  using  fracture  mechanics  techniques. 

(1)  A  molecular  model  was  developed  explaining  stress  enhanced  crack 
growth  in  vitreous  Si02.  During  the  past  year,  work  on  alkali  silicate 
glasses  has  shown  that  the  presence  of  Na+  ions  In  the  glass  may  modify  this 
mechanism,  but  that  under  most  conditions  fracture  is  still  controlled  by  the 
stress  aided  rupture  of  Si-0  bonds.  Further,  it  has  been  shown  that  the  model 
also  explains  crack  growth  in  A1203  (sapphire).  The  same  environments  which 
produced  enhanced  crack  growth  In  S102  also  did  so  in  A1203.  It  has  been 
shown  that  more  ionically  bound  solids  such  as  MgF2  show  somewhat  modified 
crack  growth  behavior  in  these  environments.  Infrared  reflection  spectroscopy 
is  now  being  used  to  establish  correlations  between  subcritlcal  crack  growth 
and  corrosion  behavior  of  glasses  In  a  number  of  liquid  environments. 

In  a  parallel  effort,  an  electrostatic  model  for  crack  growth  in  glasses 
has  been  developed.  This  model  Is  consistent  with  the  mechanism  discussed 
above.  In  the  absence  of  strong  chemical  reactions  at  crack  tips,  crack 
growth  rates  are  predicted  to  depend  on  the  dielectric  constant  of  the 
environment  and  the  stress  dependence  of  the  dielectric  constant  in  the  glass. 


This  electrostrictive  model  is  consistent  with  crack  growth  data  obtained  on 
soda-1 ime-silica  glass;  data  on  vitreous  silica  is  being  collected. 

(2)  Indentation  fracture  techniques  are  being  used  to  determine  effects 
of  the  cubic  to  tetragonal  phase  transformation  on  both  the  strength  and 
fracture  toughness  of  BaTiO^.  It  has  been  shown  that  a  £  50%  decrease  in  Kjj. 
occurs  in  going  from  the  ferroelectric  to  the  paraelectric  state.  This 
decrease  occurs  monotonically  with  an  increase  in  temperature  up  to  the  Curie 
point  120°C).  This  effect  is  reversible,  the  stress  required  to  propagate 
an  indentation  induced  flaw  returning  to  its  original  room  temperature  value 
after  heating  up  to  150°C  and  then  cooling.  This  reversibility  indicates  that 
the  toughness  variation  is  intrinsic  to  the  material  and  suggests  that  no 
change  in  the  indentation  induced  stress  field  occurs  due  to  heating  above  the 
Curie  point.  Deviations  in  strengths  from  indentation  theory  predictions  were 
observed  at  small  indentation  loads  and  flaw  sizes.  This  behavior  is 
interpreted  in  terms  of  the  influence  of  internal  stresses  induced  by  the 
cuMc  to  tetragonal  transformation  and  is  consistent  with  earlier  studies  on 
similar  materials.  Studies  are  underway  to  determine  the  effect  of  these 
localized  internal  stresses  on  subcritical  crack  growth  behavior  in  BaTiOj. 

(3)  The  objective  of  this  part  of  the  program  is  to  determine  whether 
certain  electrical  breakdowns  in  ceramic  capacitors  occur  due  to  cracks  which 
grow  between  two  electrodes  under  combined  mechanically  and  electrically 
induced  stress  fields.  Observations  of  indentation  induced  cracks  in  both  BX 
and  NPO  ceramic  capacitor  compositions  has  shown  that  there  is  a  strong 


interaction  between  cracks  and  the  metal  electrode  layers  in  the  capacitor, 
leading  to  a  significantly  increased  apparent  toughness  of  the  material. 
Dynamic  fatigue  studies  of  both  the  BX  and  NPO  capacitors  also  suggest 
significant  electrode-crack  interactions.  Low  stressing  rate  studies  on  this 
material  containing  indentation  Induced  cracks  on  which  an  electrical  field  of 
1400  v/cm  was  imposed,  show  no  effects  of  this  field  on  the  strength.  Also  no 
change  in  crack  length  for  an  Indentation  induced  crack  was  observed  on 
capacitors  under  the  above  electrical  field.  Work  on  effects  of  combined 
mechanical  loads  and  electrical  fields  on  crack  extension  is  being  continued. 
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A  Molecular  Mechanism  for  Stress  Corrosion  in  Vitreous  Silica 

TERRY  A.  M1CHALSKE*  and  STEPHEN  W.  FREIMAN*  * 

Ceramic  Development,  Sandia  Laboratories,  Albuquerque,  New  Mexico  871  IS 


The  mechanical  strength  of  moot  glasses  aad  ceramics  de¬ 
creases  with  time  under  static  loading  In  an  ambient  environ¬ 
ment.  This  strength  lorn  is  associated  with  slow  growth  of 
preexisting  surface  flaws  due  to  stress  corrosion  by  water  from 
the  surrounding  environment.  We  studied  stress  corrosion  in 
vitreous  silica  exposed  to  water  and  several  nonaqueons  envi¬ 
ronments;  environments  which  enhance  stress  corrosion  crack 
growth  in  silica  contain  active  groups  with  electron  donor  sites 
on  one  end  and  proton  donor  sites  at  the  other.  These  results 
suggest  a  detailed  chemical  model  for  the  interaction  of  the 
environment  with  mechanically  strained  bonds  in  the  solid  at 
the  tip  of  a  crack.  The  proposed  model  for  stress-corrosion 
crack  growth  also  has  Implications  for  the  long-term  strength 
behavior  of  a  wide  variety  of  brittle  materials. 

I.  Introduction 

A  decrease  in  strength  with  time  under  load  in  ambient  environ¬ 
ments  is  observed  for  most  glasses  and  ceramics .  This  phenome¬ 
non,  commonly  referred  to  as  static  fatigue  or  delayed  failure, 
results,  at  least  in  part,  from  slow  growth  of  preexisting  flaws  in 
the  materials.  The  small  flaws  grow  until  they  are  large  enough  to 
result  in  catastrophic  failure.  Up  to  now,  no  conclusive  evidence 
existed  for  stress  corrosion  in  glasses  or  ceramics  due  Co  any 
species  other  than  water.  In  addition,  models  for  crack  growth  due 
to  water,  such  as  those  proposed  by  Orowan,'  Charles  and  Hillig,2 
and  Wiederhom  el  a/.,  are  phenomenological  in  nature.  None  of 
these  theories  provides  an  understanding  of  why  water  is  especially 
effective  in  increasing  the  rate  of  slow  crack  growth,  nor  do  they 
provide  a  means  of  predicting  the  effect  of  nonaqueous  environ¬ 
ments  on  fracture  of  glass. 

In  this  paper,  we  present  a  detailed  view  of  the  chemical  inter¬ 
action  between  strained  crack-tip  bonds  in  vitreous  silica  and  water 
molecules  from  the  environment.  We  show  how  such  a  molecular 
interpretation  of  the  stress-corrosion  process  leads  to  an  under- 

Presented  «  the  Buk  Science  Division  Joint  Fall  Meeting  with  TMS-AJME, 
Louisville,  Kentucky.  October  14,  1981  (Paper  No.  106-B-8IF).  Received  June  21, 
1982;  revised  copy  received  October  4,  1982;  approved  October  21.  1982. 
Supported  in  put  by  the  Office  of  Naval  Research  and  the  Department  of  Energy. 
'Member,  the  American  Ceramic  Society. 
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Fig.  I.  Representation  of  proposed  reaction 
between  water  and  strsined  Si-O-Si  bond  u  crack 
tip.  Reaction  steps  involve  (/)  adsorption  of 
water  to  Si-O  bond,  (2)  concerted  reaction  in¬ 
volving  simultaneous  proton  and  electron  trans¬ 
fer,  and  ( 3)  formation  of  surface  hydroxyl  groups. 


standing  of  the  effectiveness  of  water  as  a  stress-corrosion  agent 
along  with  the  ability  to  qualitatively  predict  the  effects  of  various 
nonaqueous  stress-corrosion  agents.  Experimental  results  in  the 
form  of  crack  velocity  vs  stress  intensity  diagrams  will  be 
presented  for  crack  growth  in  vitreous  silica  in  support  of  this 
interpretation.  Finally,  a  brief  section  is  included  to  demonstrate 
the  coherence  of  such  a  molecular  interpretation  of  stress  corrosion 
with  current  atomistic  models  of  the  fracture  process. 


II.  Chemical  Bond  Rupture  in  Silica 

Network  silicates  are  composed  of  [SiOJ'4  tetrahedral  units 
linked  at  their  comers  by  bridging  Si-O-Si  bonds.  On  the  atomic 
scale,  it  is  the  bridging  Si-O  bond  which  supports  stress  in  the 
material  and  its  rapture  is  important  to  the  fracture  process.  At  the 
crack  tip  in  a  stressed  solid,  highly  concentrated  strain  fields  are 
produced.  Continuum  approximations  indicate  that  the  bridging 
bond  experiences  strains  >20%.  The  effect  of  this  strain  on  the 
bridging  Si-O  bond  can  be  discussed  in  terms  of  decreased  bond¬ 
ing  overlap  and.  thus,  an  increased  availability  for  Si  and  O  atoms 
to  bond  with  other  species. 

The  structure  and  bonding  in  the  water  molecule  are  established 
in  die  literature.  Oxygen  atomic  orbitals  2s.  2px,  and  2py  form 
three  hybrid  orbitals,  two  of  which  tr-bond  with  hydrogen  atoms. 
The  remaining  tone-pair  orbitals  (one  hybrid  and  the  remaining  pi) 
are  directed  away  from  the  hydrogen  atoms.  This  arrangement 
result!  in  a  positive  charge  center  opposite  the  lone-pair  orbitals. 

The  interaction  between  a  strained  bridging  bond  at  the  crack  tip 
in  silica  and  a  water  molecule  from  the  environment  can  be  repre¬ 
sented  by  a  three-step  process  (Fig.  1). 

Step  I.  A  water  molecule  from  the  environment  attaches  to  a 
bridging  Si-O-Si  bond  at  the  crack  tip.  The  water  molecule  is 
aligned  by:  (1)  formation  of  the  hydrogen  bond  with  the 
atom  and  (2)  interaction  of  the  lone-pair  orbitals  from  Olwmrr)  with 
the  Si  atom.  The  lone-pair  orbital  interaction  may  involve  either 
Van  der  Waals  attraction  or  some  covalent  bonding  with  unoccu¬ 
pied  orbitals  of  Si. 

Step  2.  A  conceited  reaction  occurs  in  which  proton  transfer  to  the 
0(*„  is  accomplished  simultaneously  with  electron  transfer  from 
the  0(.,  to  the  Si  atom.  As  a  result  of  this  reaction,  two  new 
bonds  are  formed,  one  between  0<„>  and  Si,  and  one  between 
hydrogen  and  Ottr);  the  original  bridging  bond  between  0,*,>  and  Si 
is  destroyed. 

Step  3 .  Rupture  of  the  hydrogen  bond  between  Olw,  and  transferred 
hydrogen  occurs  to  yield  surface  Si-O-H  groups  on  each  fracture 
surface.  Since  the  hydrogen  bond  is  weak,  this  step  is  expected  to 
occur  immediately  after  proton  transfer.  (Budd4  proposed  a  similar 
mechanism  for  the  dissolution  of  silica  in  water;  however,  it  should 
be  noted  that  the  present  model  for  crack  growth  is  not  contingent 
on  the  removal  of  material  from  the  fracture  surface.) 

It  is  worth  noting  that  the  chemical  mechanism  used  here  to 
describe  rapture  of  highly  strained  crack-tip  bonds  is  in  agreement 
with  reactions  observed  at  strained  siloxane  (SiOSi)  bonds  present 
on  dehydrated  silica  surfaces.9  Spectroscopic  results  show  that 
dehydrating  a  silica  surface  at  7<400°C  creates  strained  siloxane 
surface  bonds.  When  water  is  readmitted,  the  strained  surface 
bonds  act  as  sites  to  dissociatively  adsorb  water  molecules,  thus 
forming  two  silanol  (SiOH)  surface  groups. 

This  model  for  chemical  bond  rupture  has  a  number  of  inter¬ 
esting  implications.  Note  first  that  there  is  no  requirement  for  prior 
dissociation  of  the  water  molecule,  nor  do  any  reaction  products 
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Thbte  1.  Properties  of  Environments  litMUpted 


Dielsctric 

constant. 

Environment  «' 

Dinote  moment. 
(iXlO*  (Cm) 

Water  contcal 

Water.  H20 

78 

6.17 

Deuterated  water.  D20 

78 

6.20 

Ammonia,  NH> 

22 

4.90 

1-5%  rh 

Hydrazine,  N2H, 

52 

5.83 

5%  max 

Formamide,  CHjNO 

109 

12.4 

365  ppm* 

Acetonitrile,  CHiCN 

39 

13.1 

125  ppm* 

Carbon  monoxide,  CO 

1 

0.33 

1-5%  th 

Nitrobenzene.  CH6HJ'J02 

34 

14.1 

Nitrogen,  N2 

1 

0 

1-5%  rh 

•By  Karl  Fisher  technique 


need  be  removed  from  the  surface.  This  is  quite  a  different  picture 
than  that  presented  by  Charles  and  Hillig.2  The  model  is  consistent 
with  the  concepts  of  Orowan,'  who  suggested  that  the  driving  force 
for  crack  growth  is  the  reduction  in  surface  energy  by  the  environ¬ 
ment.  Second,  the  model  suggests  that  environments  other  than 
water  should  enhance  crack  growth  if  the  species  possess  structural 
and  bonding  features  similar  to  those  of  water,  i.  e.  proton  donor 
sites  at  one  end  of  the  molecule  (or  groups)  and  lone-pair  orbitals 
at  the  other.  The  environmental  molecule  (or  group)  must  also  be 
of  comparable  size  with  the  Si-O  bond  (*>0.163  nm). 

III.  Experimental  Procedure 

To  evaluate  the  stress-corrosion  reaction  proposed  in  Section  11, 
crack-growth  studies  were  conducted  on  vitreous  silica*  in  water, 
Nj  gas  of  low  relative  humidity,  and  several  nonaqueous  environ¬ 
ments  (Table  I).  The  glass  was  cut  into  double-cantilever-beam 
specimens  (50  by  12  by  1  mm)  having  a  center  groove  of  about 
one-half  the  specimen  thickness.  The  specimens  were  annealed  at 
I090°C  after  machining.  A  sharp  crack  was  introduced  into  one  end 
by  pressing  a  sharp  point  onto  the  side  opposite  the  groove.  Alumi¬ 
num  arms  were  epoxied  to  the  specimen  and  crack-growth  data 
obtained  using  the  applied-moment  technique*  by  applying  dead¬ 
weight  loading  and  using  a  40X  traveling  microscope  with  a  filar 
eyepiece  to  monitor  crack  motion.  All  tests  other  than  those  in  HjO 
and  deuterated  water  were  conducted  inside  a  Plexiglass*  chamber. 
Prior  to  testing,  the  chamber  was  purged  with  Nj  gas  for  5  to 
10  min.  In  the  case  of  NHi  and  CO,  these  gases  were  then  con¬ 
tinuously  bled  at  a  constant  flow  rate  into  the  chamber  during 
testing.  For  the  liquid  environments,  an  Nj  overpressure  was  main¬ 
tained  to  prevent  the  pickup  of  moisture  from  the  environment.  The 
effect  of  the  water  always  present  in  any  environment  on  the 
interpretation  of  the  crack-growth  data  is  discussed  in  the  follow¬ 
ing  section. 

The  measured  crack  velocities  are  plotted  vs  stress  intensity 
factor  (£i)  in  Fig.  2  for  vitreous  silica  tested  in  liquid  water  and  N, 
gas  of  1  to  2%  relative  humidity.  Note  that  the  functional  form 
V-V0  cxpibK,)  was  used  throughout,  since  this  expression  has  a 
better  theoretical  basis  than  the  commonly  used  power  law.  How¬ 
ever,  the  conclusions  regarding  the  data  would  not  be  altered 
through  use  of  the  latter.  The  important  features  in  this  plot  are: 

(1)  Crack  growth  in  the  low- velocity  regime  in  Nj  gas  is  en¬ 
hanced  (curve  is  shifted  to  lower  K,  and  V  is  a  measurable  function 
of  hi)  by  the  small  amount  of  water  present,  as  was  shown  pre¬ 
viously  by  Wiederhom1  for  soda-lime  glass.  The  relative  amount 
of  enhancement  is  related  to  the  chemical  activity  of  water  (i.  e. 
relative  humidity)  in  the  environment.1 

(2)  A  plateau  in  crack  velocity  is  observed  in  the  moist  N2.  This 
plateau  has  been  linked  to  crack  growth  limited  by  the  rate  of  water 
transport  to  the  tip  of  the  moving  crack.1''  At  sufficiently  high  V, 
the  plateau  behavior  ends  and  rapid  fracture  is  observed. 

(3)  In  N2  gas  above  the  plateau  velocity,  the  V-*,  curve  was  too 
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Fig.  2.  Effect  of  ante*  on  ikm  fracture  in  vitreous  silica 
(curves  represent  best  fit  to  experimental  data). 


Stress  Intensity  K,  (MPa>m'/i) 

Bpk  K,  vs  V  diagraajfar  vtosoasriliw  at^wemigqgiiMre. 
rate  of  slow  crack  growth. 


sleep  to  obtain  measurements,  in  agreement  with  the  spontaneous 
failure  of  vitreous  silica  observed  by  Wiederhom  et  at*  in  a 
vacuum  of  0.133  Fs. 

(4)  Liquid  water  enhances  crack  growth  over  the  entire  range  of 
crack  velocities  (10‘*  to  10'*  m/s).  (At  velocities  >10  1  m/s, 
cavitation  will  occur, N  and  no  effect  of  the  water  will  be  observed. ) 

Figure  3  shows  the  crack-growth  curves  for  ammonia,  hy¬ 
drazine  ,  and  fotmamide  along  with  die  curves  for  N,  gas  and  water 
from  Fig.  2.  The  data  points  for  water  and  ammonia  are  shown  to 
indicate  the  amount  of  scatter  associated  with  the  experimental 
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Stress  Intensity  K|  (MPa*m,J) 

Flf.4.  K,  vs  V  diagram  for  vitreous  silica  at  room  temperature, 
showing  that  carbon  monoxide,  nitrobenzene,  and  acetonitrile  have 
little,  it  any,  effect  on  rate  of  slow  crack  growth. 


measurements;  all  other  data  are  represented  by  best-fit  lines.  The 
most  important  feature  of  this  plot  is  the  absence  of  any  plateau  in 
the  curves  for  ammonia,  hydrazine,  or  formamide.  The  plateau 
occurs  because  crack-growth  rates  controlled  by  small  amounts  of 
the  water  in  a  liquid  or  gas  eventually  become  limited  by  the  rate 
at  which  the  water  can  diffuse  to  the  crack  tip.  For  a  water- 
saturated  gas  this  plateau  would  occur  at  10  *4  m/s  (Ref.  7);  it 
would  be  significantly  lower  in  a  liquid  because  of  viscosity  con¬ 
tributions.1  For  a  relative  humidity  of  0. 17%  in  N2  gas,  the  plateau 
would  occur  at  10  7  m/s. 7  This  water  content  is  below  that  known 
to  occur  in  any  of  the  above  environments  (Table  I).  The  absence 
of  a  plateau  in  velocity  leads  to  the  conclusion  that  the  reaction  of 
ammonia,  hydrazine,  and  formamide  at  the  crack  tip  controls  crack 
velocities  in  preference  to  water  present  in  these  environments. 
The  similarity  in  the  slopes  of  the  V-K,  curves  in  these  environ¬ 
ments  to  that  in  water  suggests  that  stress-activated  processes  are 
also  similar.  It  is  tempting  to  attribute  significance  to  the  overlap 
of  the  data  for  water  and  ammonia.  However,  one  must  realize  that 
while  the  ammonia  is  present  as  a  gas  at  essentially  atmospheric 
pressure,  and  so  has  an  activity  of  1,  the  activity  of  the  water 
referred  to  the  same  state  is  only  ~ 30/760. 

As  Fig.  4  indicates,  neither  carbon  monoxide,  acetonitrile,  nor 
nitrobenzene  controlled  crack  growth  in  preference  to  the  small 
quantities  of  water  present.  The  V-K,  diagram  for  each  displays  a 
plateau  at  the  approximate  velocity  expected  for  the  amount  of 
water  present.  It  should  be  noted  that  the  crack-growth  curves  in 
region  III  vary  for  each  environment,  possibly  reflecting  electro¬ 
static  effects  of  each  species  on  the  firacture  process,  as  described 
by  Wiederhom  et  at. 1 

Finally,  Fig.  3  compares  crack-growth  rates  in  normal  water  and 
deu tented  water  at  room  temperature.  It  can  be  seen  that  the  slopes 
of  the  V-K i  curves  are  parallel  but  the  curve  for  demented  water 
is  displaced  to  lower  velocities. 

IV.  Discussion 

(!)  Implications  of  Model  for  Silica 

Figure  3  shows  that  water,  ammonia,  hydrazine,  and  formamide 
are  effective  in  promoting  stress-corrosion  crack  growth  in  vitreous 
silica.  Each  of  these  environmental  species  is  similar  in  that  its 
molecular  structure  contains  die  combination  of  lone-pair  electron 
orbitals  opposite  a  proton,  with  a  separation  distance  of  H).  I  run. 
This  structural  similarity  is  important  since  it  indicates  that  each  of 
the  environments  that  enhanced  strew  corrosion  is  capable  of  par- 
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Fig.  5.  Comparison  of  effects  of  HjO  and  DjO  on  slow  crack 
growth  in  vitreous  silica  at  room  temperature. 


ticipating  in  the  concerted  reaction  described  in  Fig.  I .  The  sche¬ 
matic  drawings  of  Fig.  6  show  the  molecular  structure  of  each  of 
the  effective  stress -corrosion  agents.  The  dashed  line  represents  the 
bond  which  is  expected  to  cleave  during  the  crack-tip  reaction.  In 
each  case,  one  of  the  resulting  fragments  is  expected  to  be  a  proton 
which  will  form  a  silanol  surface  group.  However,  each  stress- 
corrosion  environment  presents  a  different  remaining  fragment 
which  is  expected  to  form  a  surface  compound  by  bonding  with  a 
silicon  ob  the  glass  surface.  The  difference  in  surface  compounds 
for  silica  fractured  in  different  stress-corrosion  environments 
presents  an  interesting  method  for  potentially  verifying  the  pro¬ 
posed  model  for  stress-corrosion  chemistry.  We  are  currently  pur¬ 
suing  surface  analysis  techniques  which  will  allow  identification  of 
these  surface  compounds. 

So  far  we  have  shown  that  environmental  species  possessing 
lone-pair  orbitals  opposite  proton  donor  sites  promote  stress  corro¬ 
sion  in  vitreous  silica.  What  is  equally  important  is  the  evidence 
that  species  devoid  of  this  structural  similarity  do  not  participate 
strongly  in  the  stress -corrosion  process.  As  demonstrated  in  Fig.  4. 
carbon  monoxide,  acetonitrile,  and  nitrobenzene  show  only  en¬ 
hancement  of  low-velocity  crack  growth,  which  can  be  directly 
related  to  the  low-level  concentration  of  water  in  each.  Although 
each  of  these  environments  contains  lone-pair  orbitals,  none  has  a 
proton  donor  site  directly  across  and,  except  for  carbon  monoxide, 
the  lone  pair  is  pointing  in  toward  the  center  of  the  molecule. 

Since  reference  has  been  made  to  correlations  between  dielectric 
constant"  and  dipole  moment12  of  the  environment  and  slow  crack- 
growth  rates,  we  examined  our  K,  vs  V  data  for  such  relations. 
Table  I  lists  the  environments  studied  in  this  work  along  with  their 
respective  dipole  moment  and  dielectric  constant.  No  systematic 
correlations  were  found  between  the  values  from  Ihble  I  and  the 
slow  crack-growth  results  of  Figs.  3  and  4.  Nitrobenzene  (dipole 
moment  1.41  x  10' 20  Om)  did  not  control  slow  crack-growth 
rates,  whereas  ammonia  (4.9x  10'”)  Cm  was  a  strong  stress- 

■  corrosion  agent.  Formamide  (e'=  109)  was  less  effective  as  a 
stress-corrosion  agent  than  ammonia  (e’(/)* 26)  (e'(g)=l).  It  is 
apparent  from  these  comparisons  and  others  which  can  be  made 
from  the  present  data  that  studies  which  report  effects  of  high 
dipole  moment  or  high  dielectric  constant  on  strength  may  in  fact 
be  reporting  the  effect  of  water  impurities  in  those  environments. 

Figure  3  shows  that  crack-growth  rates  measured  in  deuterated 
water  are  lower  than  those  measured  in  normal  water.  To  interpret 
this  hydrogen  isotope  effect,  it  will  be  necessary  to  discuss  the 
stress-corrosion  model  in  the  context  of  reaction-rate  theory.  This 
discussion  can  begin  by  examining  a  schematic  of  the  reaction 
coordinate  for  crack-tip  bond  rupture  in  a  silica  specimen  which  is 
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Fig.  6.  Representations  of  water  and  other 
nonaqueous  species  having  tone-pair  orbitals 
(electron  donor  sites)  directly  across  from 
protons  (electron  acceptor  sites).  Dashed  line 
indicates  bond  expected  to  cleave  during  crack- 
tip  interaction. 


Reaction  Coordinate 

Fig.  7.  Representation  of  reaction  coordinate 
for  bond  rupture  in  silica;  solid  line  indicates 
rupture  in  vacuum  and  dashed  line  rupture  with 
water  present. 


strained  below  the  critical  stress  intensity  for  catastrophic  failure 
(see  Fig.  7).  The  main  points  of  interest  on  the  reaction  coordinate 
are  the  energy  levels  of  the  reactants  (bridging  bond),  the  activated 
complex  (that  complex  which  exists  at  the  time  of  bond  rupture), 
and  the  products  (those  species  left  on  the  fracture-generated  sur¬ 
faces).  In  vacuum  (solid  line),  the  entire  energy  of  the  bridging 
bond  must  be  overcome  before  bond  rupture  can  occur.  (This 
energy  is  supplied  in  pan  by  mechanical  strain  as  well  as  by 
thermal  contributions).  When  water  is  present  (dashed  line),  the 
height  of  the  energy  barrier  to  bond  rupture  is  reduced  because 
charge  transfer  (step  2  of  crack-tip  interaction)  provides  an  alter¬ 
nate,  lower-energy,  reaction  path.  The  effect  of  a  reduced  energy 
harrier  to  reaction  can  be  seen  in  Eq.  (1). 

k=(RT/Nh)  expf-UG'/RT)  (1) 

where  *  is  the  reaction-rate  constant,  N  Avogadro’s  number,  h 
Planck's  constant,  T  temperature,  R  the  gas  constant,  and  AG*  the 
height  of  the  energy  barrier  between  the  reactants  and  the  activated 
complex.  Since  the  rate  constant  for  a  chemical  reaction  (in  this 
case,  bond  rupture)  is  reciprocally  dependent  on  the  height  of  the 
energy  barrier,  mechanisms  that  provide  a  lower-energy-activated 
complex  will  increase  crack-growth  rates.  It  should  be  noted  that 
factors  such  aS  the  change  in  the  energy  of  the  bridging  bond  due 
to  adsorbed  water  (step  I  of  crack-tip  interaction)  and  the  energy 
of  the  reaction  products  (step  3  of  crack-tip  interaction)  can  affect 
the  energy  barrier  to  fracture  and,  thus,  the  fracture  rate.  It  will  be 
important  in  subsequent  studies  to  evaluate  the  relative  importance 
of  all  three  steps  of  the  bond-rupture  process  on  the  fracture  rate. 
At  present,  the  observed  isotope  effect  indicates  that  the  concerted 
reaction  (simultaneous  electron  and  proton  transfer)  is  the  rate- 
limiting  step  to  the  chemical  bond  rapture  process.  The  isotopic 
effect  arises  from  the  difference  in  ground-state  vibrational  ener¬ 
gies  for  deuterium  and  protium  which,  in  turn,  results  in  different 
activation  barriers  for  the  two  isotopes.  When  a  proton  exchange 
is  involved  in  the  rate-limiting  step,  a  primary  (strong)  isotope 
effect  is  predicted  whereby  the  deuterated  reaction  is  expected  to 
proceed  more  slowly.  Since  we  do  see  a  decrease  in  crack  velocity 
for  deuterated  water,  it  is  reasonable  to  conclude  that  proton  trans¬ 
fer  is  involved  in  the  rate-limiting  step  to  fracture. 

The  preceding  discussion  has  addressed  reaction-rate  theory  in 
an  elementary  manner.  This  topic  is  treated  more  rigorously  by 
solvation  theory,  which  parallels  closely  the  reaction  involved  in 
stress-corrosion  crack  growth.  Although  solvation  chemistry  and 
stress-corrosion  chemistry  are  not  identical,  both  involve  reactions 
which  break  down  bridging  bonds  in  the  solid.  (In  solvation,  there 
is  a  need  to  contain  the  fragments  of  bond  rapture  in  solution. 
However,  this  step  is  not  necessary  in  stress-corrosion  crack 
growth  since  broken  bonds  are  separated  by  mechanical  strain.) 
Reichardt's11  treatment  of  solvation  theory  presents  a  detailed  dis¬ 


cussion  of  possible  interactions  between  solvent  molecules  and 
chemical  reactions  occurring  in  the  solute.  His  treatment  shows 
bow  various  interactions  can  influence  rate  constants  for  the  chemi¬ 
cal  interactions.  Application  of  solvation  theory  to  stress-corrosion 
effects  can  be  very  helpful  in  understanding  the  phenomenon  of 
slow  fracture.  However,  it  should  be  noted  that  phenomenological 
studies  of  chemical  kinetics  alone  do  not  provide  information  con¬ 
cerning  details  of  the  chemical  reaction  mechanism  and  may  not  be 
extremely  useful  for  predicting  fracture  rates  in  other  environments 
or  for  other  brittle  materials. 

(2)  Predictions  for  Other  Materials 

We  have  shown  that  the  mechanistic  model  for  stress  corrosion 
in  silica  can  be  used  to  predict,  qualitatively,  the  effectiveness  of 
various  nonaqueous  environments.  Recent  work  has  shown  that  the 
mechanistic  approach  is  useful  not  only  for  predicting  stress- 
corrosion  behavior  of  other  oxide  materials  but  can  be  extended  to 
ionic  as  well  as  purely  covalent  solids.  Michalske  et  at. 14  showed 
that  single-crystal  sapphire  (AljOj)  shows  slow  fracture  behavior 
similar  to  silica  when  exposed  to  the  nonaqueous  environments 
discussed  in  this  paper.  This  result  was  prcdirtod  since  the  chem¬ 
istry  of  the  Si-O  bond  and  the  Al-O  bond  are  quite  similar. 
However,  the  same  study  showed  that  single-crystal  magnesium 
fluoride  (MgFi)  exhibits  quite  a  different  susceptibility  when  ex¬ 
posed  to  the  same  environments.  Specifically,  ammonia,  which 
was  very  effective  in  promoting  stress  corrosion  in  the  oxide  mate¬ 
rials,  was  not  effective  for  magnesium  fluoride,  whereas  ace¬ 
tonitrile,  which  wu  not  effective  for  the  oxides,  greatly  enhanced 
fracture  rates  in  magnesium  fluoride.  These  results  were  inter¬ 
preted  in  terms  of  the  change  in  bood  character  in  going  from  oxide 
to  fluoride  materials.  In  fluorides,  the  bridging  bond  is  mostly  ionic 
in  character  and  for  this  reason  electrostatic  interactions  with  the 
environmental  molecules  are  expected  to  dominate  stress-corrosion 
behavior.  This  type  of  interaction  is  quite  different  from  the  chemi¬ 
cal  interaction  proposed  for  the  ox'ide  system.  Experimental  in¬ 
vestigations  of  stress  corrosion  in  silicon”  suggest  that  aqueous 
environments  have  little,  if  any,  effect  on  fracture  rate.  This  result 
would  also  be  predicted  from  the  mechanistic  model  since  the 
silicon-silicon  bond  is  completely  covalent  and  would  not  exhibit 
.  a  preference  for  either  end  of  the  water  molecule  (lone-pair  orbital 
or  proton)  and,  thus,  would  not  easily  participate  in  the  concerted 
reaction  hypothesized  for  stress  corrosion  in  oxides. 

Finally,  a  mechanistic  model  for  the  stress-corrosion  process  has 
the  potential  to  interface  with  current  atomistic  models  concerning 
the  physics  of  fracture.  Theoretical  work  by  Fuller  et  at. "  showed 
how  any  arbitrary  chemical  interaction  can  be  placed  in  the  context 
of  a  crack  tip  in  a  discrete  atomic  lattice.  It  will  be  important  in  the 
Aiture  to  add  experimental  parameters  obtained  in  conjunction  with 
mechanistic  models  for  crack  growth  with  theoretical  models  for 
crack-tip  physics  in  order  to  make  quantitative  predictions  of  stress 
corrosion  effects. 
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V.  Summary 

Crack  velocities  were  measured  as  a  function  of  for  vitreous 
silica  in  water  and  several  nonaqueous  environments.  Results  of 
these  studies  showed  that  environments  with  a  molecular  structure 
containing  a  lone-pair  electron  orbital  across  from  a  proton  donor 
site  have  the  greatest  effect  on  crack  growth.  A  qualitative  model 
is  proposed  to  demonstrate  how  molecules  of  this  structural  type 
may  interact  with  Si-0  bonds  at  the  crack  tip.  This  model  enables 
one  to  predict  the  effect  of  a  given  environment  on  crack-growth 
rates  in  silica.  This  molecular  model  can  also  be  used  to  predict 
stress-corrosion  effects  in  other  brittle  materials. 
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The  influence  of  transport-limited  kinetic  crack  growth  oa  the 
fatigue  properties  of  soda-lime  glam  was  examined.  Dynamic 
fatigue  data  were  taken  on  apsclnsaaa  with  controlled  inden¬ 
tation  flaws  and  were  compared  with  the  predicted  respome 
from  measured  crack  velocity  characteristics.  Heptane  was 
used  as  the  operational  test  environment  because  of  its  pro¬ 
nounced  crack  velocity  plateau;  control  tests  in  water  served  to 
establish  a  baseline  reference  for  comparing  the  results.  Free- 
tographic  observations  using  a  stress  wave  marker  technique 
showed  a  complex  growth  history  for  flaws  broken  in  heptane 
compared  to  that  for  flaw*  broken  in  water.  The  magnitude  at 
the  predicted  region  II  influence  is  too  small  to  be  detected  in 
the  dynamic  fatigue  results,  even  allowing  for  the  relatively 
Ugh  degree  of  data  reproducibility.  The  impHcatioas  of  this 
conclusion  for  lifetime  predictions  are  discussed. 

I.  Introduction 

The  long-term  strength  of  brittle  solids  is  governed  by  such 
factors  as  flaw  size  and  shape  and  the  susceptibility  to  chemi¬ 
cally  enhanced  slow  crack  growth.  Implicit  in  the  fracture  me¬ 
chanics  treatments  of  time-dependent  strength  (fatigue)  properties 
are  certain  assumptions  concerning  these  factors,  e.g.  that  the 
flaws  respond  in  essentially  the  same  way  as  macroscopic  cracks 
and  that  the  crack  velocity  can  be  expressed  as  some  simple 
function  (usually  power  law)  of  the  stress  intensity  factor.  The 


Received  October  14.  I9S2;  reviled  copy  received  Februery  II,  19(3;  approved 
Febnuiy  21,  1983. 

Supported  by  the  U  S.  Office  of  Naval  Research.  Metallurgy  end  Ceramics  Pro¬ 
gram.  and  by  the  Deutsche  Fonchuoiigemeinschsft 
'Member,  the  American  Ceramic  Society. 

•Now  with  the  Department  of  Physics.  Chulahmgkon  University.  Bangkok, 
Thailand. 

’On  leave  from  Fraunhofer  InetittM  fttr  Wertaaoffmechanik.  Freiburg,  Federal 
Republic  of  Germany. 


success  of  this  approach  has  led.  in  conjunction  with  statistical 
accountability  of  flaw  populations,  to  useful  engineering  design 
schemes.' 1 

This  paper  examines  the  degree  to  which  macroscopic  crack 
velocity  data  can  be  used  to  determine  the  kinetics  of  failure  in 
strength  testing.  In  particular,  a  critical  look  is  taken  at  the  potential 
complications  in  lifetime  predictions  when  the  crack  velocity  func¬ 
tion  shows  multiregion  behavior.’  Generally,  three  such  regions  are 
identifiable:  region  I  (low  velocities),  controlled  by  rate  of  reaction 
between  environmental  species  and  crack-tip  bonds;  region  II  (in¬ 
termediate  velocities),  controlled  by  rate  of  transport  of  environ¬ 
mental  species  to  the  tip;  region  III  (high  velocities),  controlled  by 
electrostatic  environment-bond  interactions.4  It  is  customary  to 
recognize  the  existence  only  of  region  I  in  the  formal  derivation  of 
fatigue  equations,  on  the  grounds  that  it  is  the  domain  of  slowest 
growth  thin  must  control  the  crack  evolution  to  failure  and  that  in 
the  mote  concentrated  active  environments  (e.g.  water)  the  higher 
regions  are  not  strongly  evident  in  the  crack  velocity  response.  In 
addition,  it  has  been  argued  that  natural  surface  flaws,  by  virtue  of 
their  relatively  small  size  and  their  continuous  accessibility  to  the 
external  environment  at  the  points  of  intersection  with  the  free 
surface,5"7  may  not  be  subject  to  conventional  region  II  effects. 
Considerations  of  this  kind  clearly  raise  questions  about  the  status 
of  macroscopically  determined  crack  velocity  functions  as  a  suit¬ 
able  basis  for  fatigue  strength  analysis. 

The  approach  adopted  here  is  to  tun  dynamic  fatigue  tests  on 
glass  specimens  with  controlled  flaws,  generated  by  Vickers  inden¬ 
tations,*'"  in  an  environment  with  a  pronounced  region  II  crack 
velocity  plateau.  Although  indentation-induced  crack  systems  are 
subject  to  intense  residual  stress  fields,10 11  ”  these  fields  have 
been  well  characterized  and  are  readily  accommodated  into  the 
fracture  mechanics  formalism.*14  The  element  of  control  signifi¬ 
cantly  reduces  the  statistical  aspect,  and  allows  for  monitoring  of 
the  crack  evolution  to  failure,  in  the  present  study,  this  monitoring 
is  accomplished  by  imposing  periodic  stress  markers  onto  the  ffac- 
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ng.i.  Schematic  of  transducer  setup  for  imposing  sonic 
markers  on  surface  of  indentation  flaw  in  flexural  specimen. 


0.3  0.4  0.5  0.6  06 

Stress  Intensity  Factor,  KiMPa  m* 


Fig.  2.  Crack  velocity  data  for  soda-lime  glass  in  heptane 
(6  specimens)  and  water  (2  specimens)  obtained  using 
double-cantilever  technique.  For  heptane,  solid  lines  are 
segmented  least-squares  fits  to  regions  I,  II,  and  HI;  broken 
line  is  extrapolation  of  region  I  to  point  of  intersection 
with  region  III.  For  water,  solid  line  is  prediction  from 
fatigue  data. 


ture  surfaces.  I5~17  The  dynamic  fatigue  data  are  then  compared  to 
predictions  calculated  from  v-K  curves*  obtained  on  double- 
cantilever-beam  specimens. 

II.  Experimental  Procedure 

Soda-lime  glass  was  selected  as  a  model  material  for  all  tests  in 
this  study.  Heptane  was  chosen  as  the  test  medium  because  its  low 
affinity  for  water  (—50  ppm)  leads  to  a  wide  region  0  plateau; 
also,  it  has  been  determined  that  exposure  to  air  does  not  lead  to 
an  increase  in  water  concentration.  Control  tests  were  conducted  in 
water  to  provide  a  convenient  reference  baseline  for  data  analysis. 
Inert  strength  tests  for  parameter  calibration  were  conducted  in  dry 
nitrogen  gas.6'7 

(I)  Macroscopic  Crock  Velocity  Measurements 

Crack  velocities  were  measured  on  specimens  (75  by  12.5  by 
I  mm)  cut  from  glass  microscope  slides,  using  the  double¬ 
cantilever  configuration  with  applied-moment  loading. 16  A  groove 
-0.5  mm  deep  was  cut  down  the  length  of  one  face  to  guide  the 
crack.  The  specimens  were  annealed,  precracked,  and  immersed  in 
the  test  fluid.  An  optical  microscope  was  used  to  monitor  the  crack 
growth.  Stress  intensity  factors  were  calculated  from  the  applied 
load  and  specimen  dimensions." 

*Tha  mtaciipt  it  dropped  than  Kt,  on  the  undantandina  Urn  w«  are  alwaya  dealing 
wkh  node  I  cnclu 


(2)  Failure  Tests  on  Controlled-Fiaw  Specimens 

Controlled  flaws  for  strength  testing  were  introduced  by  indent¬ 
ing  with  a  Vickers  pyramid.  Raws  of  different  dimensions  were 
produced  by  varying  the  contact  load,  P,  so  that  any  size  effect  in 
the  region  □  kinetics  might  be  investigated. 1 1  All  indentations  were 
made  in  air  at  a  fixed  load  duration  of  10  s,  and  were  left  to  stand 
for  —30  min  before  strength  testing  to  allow  any  relaxation  in  the 
residual  stress  field  (e.g.  due  to  lateral  crack  growth')  to  equalize 
for  all  specimens. 

The  failure  tests  were  conducted  in  four-point  flexure,  using  a 
crosshead  machine  to  deliver  the  bending  loads.  The  rods  were 
carefully  oriented  so  that  one  set  of  radial  cracks  emanating  from 
the  indentation  comers  was  normal  to  the  maximum  tensile  stress. 
For  water,  a  droplet  was  simply  placed  on  the  indentation  site 
immediately  prior  to  testing;  in  the  case  of  heptane  a  liquid  bath 
was  necessary  because  of  rapid  evaporation.  In  the  latter  case  the 
specimens  were  dried  for  several  minutes  in  hot  air  before  im¬ 
mersion.  Beam  theory  was  used  to  compute  the  stress  at  failure,  ot. 
Routine  microscopic  examination  of  die  broken  parts  was  carried 
out  to  confirm  that  failure  originated  at  the  indentation  flaw. 

(A)  Stress  Wave  Fractograpky:  Stress  wave  fractographic 
observations'1'17  were  made  on  annealed  glass  specimens  (ISO  by 
8  by  2  mm)  containing  5  N  indentations.  A  transducer  mounted  at 
the  end  of  the  bend  bar  generates  transverse  stress  waves  normal 
to  the  plane  of  the  indentation  crack  which  is  to  lead  to  failure 
(Fig.  1).  These  waves  periodically  modify  the  direction  of  the 
maximum  tensile  stress  generated  by  the  external  load,  without 
significantly  affecting  the  driving  force  on  the  crack.  The  per¬ 
turbations  leave  optically  detectable  time  markers  on  the  fracture 
surface.  Ifcken  in  conjunction  with  the  transducer  frequency,  the 
markers  constitute  a  pictorial  record  of  the  velocity  history  of 
the  crack  system.  The  major  departure  from  previous  studies  using 
this  technique  is  in  the  lower  frequencies  attainable,  <0.1  Hz 
(cf.  the  MHz  region  conventionally  used);  under  such  operating 
frequencies  optically  resolvable  markers  could  be  produced  at 
crack  velocities  <1  /urn'1,  corresponding  to  the  heptane  data 
plateau  region. 

(B)  Dynamic  Fatigue  Curves:  Fatigue  tests  were  carried  out 
at  constant  stressing  rates,  6-..  on  annealed  glass  rods  5  mm  in 
diam.  (inner  span  20  mm,  outer  span  60  mm).  An  effort  was  made 
to  cover  as  wide  a  stressing-rate  range  as  possible  (e.g.  by  incorpo¬ 
rating  a  piezoelectric  load  cell  into  the  system  for  measuring  flex¬ 
ural  forces  at  fast  rates)7  to  maximize  the  prospects  of  detecting  any 
significant  shifts  in  the  comparative  heptane  and  water  fatigue 
curves.  Inert  strengths  were  measured  at  the  fastest  loading  rates  in 
flowing  dry  N2  gas. 

m.  Results  and  Discussion 

Let  us  now  show  the  correspondence  between  the  three  sets  of 
observations,  i.e.  macroscopic  crack  growth,  fractograpky,  and 
dynamic  fatigue. 

The  macroscopic  crack-growth  data  in  Fig.  2  show  clear  evi¬ 
dence  for  multiregion  velocity  behavior  in  heptane.  These  data 
agree  with  those  previously  reported  for  heptane  and  other  al¬ 
kanes."  For  a  macroscopic  crack,  the  level  of  the  plateau  is  a 
function  of  the  concentration  of  water  and  the  viscosity  of  the 
fluid.4  However,  Quackenbush  and  Frechette10  and  Richter16 
showed  that  both  the  crack  velocity  at  which  a  plateau  is  observed 
and  the  crack-front  shape  are  functions  of  specimen  thickness.  The 
work  of  these  authors  suggests  that  the  heptane  v-K  curve  in  Fig.  2 
is  not  unique.  Whereas  regions  I  and  III  represent  reaction  rates 
intrinsic  to  the  material-environment  system,  region  II  is  de¬ 
pendent  on  crack  size  and  shape. 

Now  consider  the  flaw-growth  patterns  in  Fig.  3  obtained  by 
stress  wave  fractography  on  cracks  grown  under  a  constant  bending 
load.  In  both  Figs.  3(A)  and  3(B),  corresponding,  respectively, 
to  heptane  and  water,  two  frequencies  have  been  superimposed  to 
allow  coverage  of  a  wide  range  of  velocities.  First  note  that  the 
markers  for  die  specimen  tested  in  water  remain  closely  elliptical 
in  profile,  and  increase  smoothly  in  spacing,  over  the  entire  range 
of  crack  growth,  corresponding  to  velocities  from  2  pm*s~'  to 
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100  mm-s  '1,  indicative  of  a  single  region  of  crack  propagation. 
The  flaw  growth  in  heptane  is  more  complex.  The  spacing  of  the 
markers  in  the  low-frequency  domain  is  nearly  constant,  and  corre¬ 
sponds  to  a  crack  velocity  of  3  fim-s  near  the  region  11  plateau. 
As  crack  growth  proceeds  within  this  region,  the  crack  front  be¬ 
comes  increasingly  distorted  from  its  initially  near-elliptical  profile 
below  the  tensile  surface.  Once  into  the  high-frequency  domain, 
corresponding  to  velocities  >  1  m-s' the  front  reverts  to  the  ellip¬ 
tical  geometry.  This  difference  in  behavior  between  water  and 
heptane  environments  was  reproducible  over  a  number  of  tests. 

Based  on  the  preceding  fractographic  observations,  it  would  be 
expected  that  the  fatigue  strengths  of  specimens  tested  under  water 
and  heptane  might  differ,  especially  at  stressing  rates  for  which  the 
crack  spends  the  largest  portion  of  its  growth  in  region  11. 

Accordingly,  let  us  examine  the  possible  correspondence  be¬ 
tween  v-K  and  dynamic  fatigue  results  in  terms  of  indentation 
fracture  mechanics.  All  mathematical  details  involved  in  establish¬ 
ing  this  correspondence  are  relegated  to  the  Appendix.  We  simply 
note  at  this  point  that  analytical  solutions  of  the  fatigue  equations 
for  flaws  with  residual  stresses  are  obtainable  only  for  single¬ 
region  crack  velocity  functions  of  power-law  form,  and  that  these 
solutions  are  themselves  of  power-law  form.14  Figure  4  shows  the 
dynamic  fatigue  curves  predicted  from  the  crack-growth  data  in 
heptane,  with  and  without  a  region  II  plateau  included  in  the  calcu¬ 
lations.  as  well  as  the  experimental  dynamic  fatigue  data  in  water 
and  heptane.  Note  that  the  data  have  been  normalized  for  inden¬ 
tation  load,  P.  so  that  all  data  for  each  environmental  condition 
might  be  reduced  to  a  universal  curve. 

The  sequence  of  operations  to  obtain  the  predicted  and  measured 
dynamic  fatigue  curves  in  Fig.  4  was  as  follows: 

(i)  On  the  assumption  that  the  results  for  tests  in  water  may  be 
described  by  single-region  behavior  over  the  data  range  covered, 
a  linear  least-squares  fit  was  appropriately  made  to  the  loga¬ 
rithmically  plotted  dynamic  fatigue  data"  in  Fig.  4.  The  slope  and 
intercept  of  this  fitted  line,  in  conjunction  with  inert  strength 
parameters,  gives,  via  a  suitable  set  of  transformation  equations, 
corresponding  values  for  exponent  and  coefficient  in  the  glass/ 
water  crack  velocity  equation. 

(ii)  A  linear  representation  of  the  result  from  step  (i)  was  made 
on  the  logarithmic  crack-velocity/stress-intensity-factor  diagram. 
Fig.  2.  The  plotted  line  is  seen  to  pass  through  the  experimental 
points  for  a  water  environment,  within  the  scatter  over  the  data 
range  covered.*  Thus  self-consistency  between  results  from  fatigue 
and  velocity  tests  is  established  for  this  single-region  system. 

(iii)  Linear  least-squares  fits  were  made  to  each  of  the  three 
clearly  defined  crack  velocity  regions  in  Fig.  2.  Appropriate  ex¬ 
ponents  and  coefficients  were  evaluated  for  each  fitted  segment. 

(iv)  Using  these  calibrated  velocity  parameters,  together  with 
the  same  inert-strength  parameters  referred  to  in  (i),  numerical 
solutions'2  were  obtained  for  the  dynamic  fatigue  equation  (Appen¬ 
dix)  for  heptane.4  Two  such  solutions  are  plotted  in  Fig.  4,  one 
with  and  the  other  without  region  II  included  in  the  crack  velocity 
function;  in  the  latter  case  region  1  is  taken  to  operate  up  to  the 
extrapolated  intersection  with  the  region  III  curve  (Fig.  2). 

Even  with  the  relatively  high  degree  of  reproducibility  achieved 
by  using  indentation-induced  flaws,  the  magnitude  of  any  region  II 
influence  seems  to  be  too  small  to  be  unequivocally  distinguished 
in  the  dynamic  fatigue  results  for  heptane.  The  fact  that  the  data 
taken  at  different  loads  fall  onto  the  same  curve  for  each  environ¬ 
ment  indicates  that  flaw  size  alone  is  not  a  primary  factor.  Similar 
conclusions  were  reported  by  Chandan  etal..’  although  their  data 
were  less  detailed  than  those  presented  here. 

IV.  Concluding  Remarks 

The  foregoing  observations  suggest  that  indentation  flaws  do  not 
show  a  one-to-one  correspondence  with  macroscopic  cracks  in 
their  multiregion  crack  velocity  response.  This  breakdown  in  cor¬ 
respondence  appears  to  be  confined  to  region  II;  we  recall  that  for 
the  tests  in  water,  where  region  I  effectively  controls  the  kinetics 
over  the  data  range  considered,  mutual  consistency  is  obtained 
between  crack  velocity  and  dynamic  fatigue  results.  Regions  I  and 


Fig.  3.  Fracture  surfaces  of  soda-lime  glass  from  indentations  at  contact 
load  of  5  N  and  subsequent  bending  stress  of  35  MPa:  (4 )  heptane,  marker 
frequencies  60  mHz  and  104  kHz;  (fl)  water,  1.1  Hz  and  1.9  kHz.  (Super¬ 
posed  frequencies  in  each  case  simply  allow  covenge  of  extended  velocity 
range;  arrows  indicate  points  at  which  markers  correspond  to  transition 
from  lower  to  higher  frequency.)  Compare  distorted  marker  pattern  in 
heptane  with  relatively  symmetrical  pattern  in  water. 


200h  P  -  2N  o 
I  1  N  o 


10-«  10-2  10°  10*  104  10» 
Stress  Rate  Parameter,  StP/MPa  S’1  N 


Fig.  4.  Normalized  dynamic  fatigue  curves  for  soda-lime  glass  in  heptane 
and  water.  Each  data  point  is  mean  and  standard  deviation,  in  logarithmic 
coordinates,  of  6  lo  15  specimens  (error  bars  omitted  for  water,  for  clarity). 
Contact  loads  P  used  for  heptane  are  indicated  (loads  used  to  obtain  control 
water  data  covers  range  0.05  to  10  N  (Ref  1 1 ),  but  are  not  differentiated 
here).  Shaded  band  denotes  inert-strength  level.  Solid  line  through  water 
data  is  least-squares  fit.  Solid  line  through  heptane  data  is  prediction  from 
crack  velocity  curve  with  region  II  included;  broken  line  is  corresponding 
prediction  with  region  II  excluded. 


'In  drawing  a  straight  line  through  the  water  data  in  Fig  2  we  intend  no  s 
about  the  fundamental  relation  between  crack  velocity  and  stress  intensity  factor,  ours 
is  strictly  an  empirical  fitting  procedure. 

'inclusion  of  a  possible  threshold  stress  intensity  lector  at  0  3  MPa-m’*3  hi  Fig.  2 
did  not  significantly  affect  the  calculated  fatigue  curves 
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III  are  intrinsic  to  the  material/environment  and  would  not  be 
expected  to  depend  in  any  way  on  crack  geometry.  For  region  II, 
however,  with  its  origins  deeply  rooted  in  transport  mechanisms, 
geometrical  factors  must  play  some  role.  Such  effects  have  been 
noted  in  fracture  mechanics  specimens,  where  the  crack  velocity 
plateau  becomes  a  function  of  specimen  width.16'20  Our  results 
suggest  that  crack  shape  is  an  even  more  important  factor  than 
crack  sire.  Thus,  in  the  present  instance  with  axisymmetric  contact 
loading  the  crack  system  always  has  a  surface  trace  in  direct  con¬ 
tact  with  the  environment,  so  transport  effects  are  apparent  only  in 
the  subsurface  propagation;  region  II  effects  are  then  manifested  as 
a  constraint  on  the  surface  expansion  by  the  more  slowly  moving 
inner  crack  portions. 

The  heptane  results  in  Fig.  2  seem  to  indicate  that  higher  regions 
in  the  crack  velocity  function  are  unlikely  to  influence  fatigue 
properties  strongly,  particularly  in  longer-term  tests,  where  the 
crack  spends  nearly  all  its  growth  time  in  region  I.  This  is  not 
altogether  unexpected,  bearing  in  mind  the  smoothing  effect  of 
integrating  over  the  crack  velocity  function  in  the  failure  mechan¬ 
ics.  Nevertheless,  care  should  be  taken  in  extrapolating  fatigue 
data  to  long  lifetime  domains  which  lie  well  outside  the  data  range. 
Fortunately,  the  conventional  linear  extrapolation  procedure  (i.e. 
without  due  regard  to  the  curvature  in  the  fatigue  plot  caused  by 
higher  regions)  will  tend  to  underestimate  lifetimes,  in  line  with  the 
requirements  of  conservative  design. 

Finally,  we  should  point  out  that  the  insensitivity  to  region  II 
effects  discussed  here  may  not  extend  to  all  strength-test  proce¬ 
dures.  In  proof  testing,  the  position  and  slope  of  region  II  have 
been  shown  to  significantly  affect  subsequent  strength  distribu¬ 
tions.  The  difference  between  proof  testing  and  dynamic  fatigue  is 
that,  in  the  latter,  the  failure  stress  is  controlled  more  by  the  lower 
end  of  the  crack  velocity  spectrum  through  which  the  flaw  grows, 
whereas  in  the  former  it  is  only  the  higher  end  which  controls. 
A  detailed  description  of  proof  testing  in  this  context  will  be 
given  elsewhere.2 


APPENDIX 


In  this  Appendix  a  summary  is  given  of  the  derivation  of  dy¬ 
namic  fatigue  equations  from  the  crack  velocity  function  for  flaws 
with  residual  contact  stress.  Reference  is  made  to  earlier  papers  for 
more  detailed  formulations.6  "  M 
The  analysis  begins  with  the  assumption  that  the  crack  velocity 
v  relates  to  the  stress  intensity  factor  K  via  a  power  law, 


v  =  v„ (K/KS 


(A-l) 


where  n  and  v0  are  empirical  exponent  and  coefficient,  re¬ 
spectively,  for  a  given  material/environment  system  and  Kc  is  the 
material  toughness.  For  cracks  of  characteristic  dimension  c  pro¬ 
duced  at  a  peak  indentation  load  P  and  subjected  to  subsequent 
tensile  stress  cr.,  the  stress  intensity  factor  has  the  form 


K=xP/c'‘2+<1kt„c' 


(A— 2) 


where  ^  and  0  are  dimensionless  constants;  the  first  term  in 
Eq.  (A— 2)  represents  the  residual  contact  field  and  the  second 
represents  the  applied  field.  In  dynamic  fatigue  testing  the  stressing 
rate  is  held  constant  as  a  function  of  time  t,  i.e. 


<r„=<7„r  (<t„= constant) 

Combining  the  above  three  equations  gives 
dc/dl  =  V'nl  xP/K,  c 1,2  +  liiaQC  l,2t/Kr  r 


(A— 3) 


(A— 4) 


which  serves  as  a  master  differential  equation.  Solution  of  this 
equation  involves  computation  of  the  time-to-failure,  tf,  i.e.  the 
time  for  the  crack  to  grow  from  its  original  size  to  a  critical 
instability  configuration  ( K  =KC.  dK/dc> 0),  thereby  defining  the 
failure  strength.  af = ir.t,. 


Solutions  of  analytical  form  are  obtainable  from  Eq.  (A-4)  only 
for  single-region  crack  velocity  functions.  FOr  the  case  where  P  is 
used  as  a  test  variable,  these  solutions  are  of  power-law  form11'14 


tTfP'^ihfaP)' 


(A-5) 


with  exponent  and  coefficient  which  relate  to  the  crack  velocity 
parameters  in  Eq.  (A-l)  via  the  "transformation  equations”14 


n=4n  73-2/3 

v0=(2nn')',2(<jmP',Y(cjP2'>)/X^ 


(A-6 a) 
(A-6 b) 


The  quantities  labeled  with  subscript  m  in  Eq.  (A-6 b)  refer  to  the 
conditions  under  equilibrium  crack  growth,  i.e.  with  K~KC, 
dK/dc  =0  in  Eq.  (A-2),  whence 


,,3=  3AC?'V44'3^(^P  )*'3 
c„/PM=(4*P/Xcfs 


(A-7 a) 
(A-7 b) 


These  instability  conditions  are  measurable  as  the  strength  and 
corresponding  crack  size  in  inert  environments.6’7  Thus  from 
Eq.  (A-S)  we  see  that  a  plot  of  log {afP l'1)  vs  \og(,d.P )  should  be 
universally  linear  for  all  contact  loads. 

For  multiregion  crack  velocity  functions  no  such  analytical  solu¬ 
tions  are  available.  It  is  then  necessary  to  integrate  Eq.  (A-4) 
numerically  over  the  crack  velocity  range.6 
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ABSTRACT 


This  paper  presents  a  review  of  our  current  understanding  of 
environmentally  Induced  slow  crack  growth  in  glasses,  single  crystals  and 
polycrystalline  ceramics.  It  Is  shown  that  the  rate  of  crack  growth  is 
controlled  by  the  chemical  activity  of  the  active  species  in  the  environment 
as  well  as  by  the  stress  Intensity  at  the  crack  tip.  A  recently  developed 
molecular  model  of  stress  Induced  chemical  reaction  between  vitreous  silica 
and  water  is  described.  The  Implications  of  this  model  for  the  effects  of 
other  chemical  species  on  crack  growth  are  discussed.  Finally,  the 
complications  introduced  by  the  presence  of  grain  boundaries  In 
polycrystalline  ceramics  are  pointed  out. 


EFFECTS  OF  CHEMICAL  ENVIRONMENTS  ON  SLOW 
CRACK  GROWTH  IN  GLASSES  AND  CERAMICS 


S.  W.  Frelman 

Inorganic  Materials  Division 
National  Bureau  of  Standards 
Washington,  DC  20234 

INTRODUCTION 

The  objective  of  this  paper  is  to  provide  a  current  view  of  the 
mechanisms  of  environmentally  influenced,  subcritical,  crack  growth  in  both 
glasses  and  ceramics.  This  review  will  concentrate  on  various  aspects  of  the 
stress  aided  chemical  reaction  which  governs  the  crack  process.  Particular 
attention  will  be  paid  to  recent  models  for  stress  assisted  physical  and 
chemical  processes  which  appear  to  govern  crack  growth. 

Before  discussing  crack  growth  in  glasses  and  ceramics,  it  should  be  said 
that  all  of  these  materials  will  be  considered  to  be  ideally  brittle,  i.e. 
that  there  are  no  zones  of  plastic  deformation  at  propagating  crack  tips. 

This  assumption  has  been  shown  to  be  valid  in  MgO  and  aluminum  oxide  [Lawn  et 
al.,  1980].  Even  for  "softer"  materials  such  as  MgFg,  where  dislocations  are 
known  to  be  generated  around  advancing  crack  tips,  fracture  is  still  governed 
by  bond  breaking  mechanisms. 

Most  of  the  crack  growth  data  reported  herein  was  taken  using  double 


n 


cantilever  beam  specimens  in  which  a  constant  load  or  bending  moment  is 


applied  to  the  specimen  and  crack  extension  measured  optically  as  a  function 

-10  “2 

of  time.  Crack  velocities  from  10  to  10  m/sec  can  be  obtained  in  this 


way.  The  data  are  usually  plotted  as  log  V  versus  Kj  ,  where  Kj  Is  the  stress 
intensity  factor  at  the  crack  tip  and  is  calculated  from  the  load  and  specimen 
dimensions. 

GENERAL  CRACK  GROWTH  BEHAVIOR 

V-Kj  plots  for  glasses,  single  crystals  and  polycrystalline  ceramics  can 
be  quite  complex  in  shape  (Figure  1).  Each  segment  can,  however,  be  described 
by  a  different  rate  controlling  mechanism.  Because  they  are  isotropic  and 
homogeneous,  most  of  the  work  in  determining  these  mechanisms  has  been 
performed  on  silicate  glasses.  These  mechanisms  can  be  synopsized  as  follows: 

Region  I. 

Because  slowly  growing  cracks  spend  most  of  their  growth  period  in  this 
regime,  it  is  of  primary  importance.  Slopes  of  curves  range  from  n  s  11  for 
some  binary  glasses  to  n  >  100  for  some  polycrystalline  ceramics.  Wiederhorn 
[1967]  showed  that  crack  velocity  is  controlled  by  the  rate  of  reaction  of  the 
active  ingredient  in  the  environment  with  the  chemical  bonds  of  the  solid  at 
the  crack  tip.  This  reaction  may  or  may  not  involve  the  formation  of 
corrosion  products. 

* 

For  engineering  purposes  and  ease  of  mathematical  manipulation  crack  growth 
data  is  often  expressed  as  V  =  A  Kj°  where  A  and  n  are  constants  which  depend 
on  material  and  environment.  While  this  expression  is  easily  integrated  to 
calculate  times  to  failure,  its  theoretical  basis  is  much  weaker  than  an 
expression  of  the  form,  V  =  VQ  exp(bKj).  Because  of  its  common  usage, 
however,  this  paper  will  use  n  as  a  measure  of  the  slope. 


While  reductions  in  surface  energy  due  to  the  attachment  of  surface 
groups  to  Si02  can  be  considered  as  the  driving  force  for  the  stress  enhanced 
reaction  at  a  crack  tip,  as  suggested  by  Orowan  [1944]  and  more  recently  by 
Parks  [1983],  the  crack  growth  rate  is  determined  by  the  kinetics  of  the 
chemical  process. 

Reaction  rate  theory  can  be  used  to  derive  a  model  for  stress 

assisted  crack  propagation  in  glasses  and  ceramics.  Such  a  derivation 

for  the  Si02~H20  reaction  (the  specific  molecular  chemistry  is  described  in 

* 

the  next  section)  is  shown  in  summary  as  follows: 

A  general  expression  for  H20  reacting  with  an  Si-0  bond  is  given 
by: 


[Si-O-Si]  +  H20  x - ►  2  [SiOH]  (1) 

where  x  is  the  activated  complex.  As  in  most  reaction  rate  theory  [Laidler, 
1965],  the  decomposition  of  the  activated  complex  to  the  products  is  assumed 
to  very  rapid  compared  to  the  reverse  reaction. 

kT  *  a(H20) 

Rate  =  k  exp  (-AG  /RT)  - j-  (2) 

where  k  is  the  transmission  coefficient  (for  most  reactions  assumed  equal  to 
1,  meaning  every  activated  complex  becomes  a  product);  k  is  the  Boltzman 
constant;  T  is  temperature;  h  is  Plancks  constant;  AG  is  the  Gibbs  free 


v.’. 
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For  a  more  detailed  derivation  see  Wiederhorn  et  al.  [1982]. 
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energy  of  activation,  and  a(H20)  is  the  activity  of  the  water  and  f  is  the 
activity  coefficient  for  the  activated  complex. 


Combining  terms  and  assuming  that  the  crack  velocity  is  directly 
proportional  to  the  reaction  rate,  we  obtain: 

* 

v  =  VQ  a(H20)  exp  (-AG  /RT)  (3) 


There  are  at  least  two  Important  aspects  to  equation  (3).  First,  note 
that  the  reaction  rate,  and  therefore  the  crack  velocity,  Is  proportional  to 
the  activity  rather  than  the  concentration  of  water  in  solution  in  the  liquid 
or  gas.  If  water  vapor  is  assumed  to  behave  as  an  ideal  gas,  then 


a.  =  !i 
1  *0 


(4) 


where  PQ  is  the  vapor  pressure  of  pure  water  at  the  temperature  of  the 
experiment.  In  a  gaseous  environment,  is  the  vapor  pressure  of  water  and 
P./P,  is  the  relative  humidity.  For  water  dissolved  in  non-aqueous  liquids, 

P^  is  the  equilibrium  vapor  pressure  of  water  over  the  solution  and  P^/P,  Is 
the  relative  humidity  of  the  gaseous  atmosphere  over  the  solution.  Given 
these  definitions  and  the  relationship  expressed  in  equation  3,  it  becomes 
clear  why  the  data  of  Wiederhorn  [1967]  and  Freiman  [1974]  obey  a  simple 
V  Z  r.h.  relationship  (Figure  2).  The  fact  that  crack  velocity  is  dependent 
on  relative  humidity  rather  than  concentration  has  Important  Implications  for 
predictions  of  structural  reliability.  It  means  that  liquids  such  as  silicone 
oil  and  fluoronated  hydrocarbons  cannot  be  considered  inert  because  their 
solubility  for  water  is  small.  Second,  this  expression  provides  a  fundamental 


xsa 

*v|Y|^1viV 


basis  for  plotting  crack  velocity  data  as  InV  vs.  Kj  rather  InV  vs.  InKj  as 


noted  earlier. 


The  stress  dependence  of  the  reaction  is  contained  in  the  AG  term.  It 
has  been  shown  (Wiederhorn  et  al.  1980)  that  this  term  can  be  expanded  as 
follows: 


*  t  *  .1/2  *  * 

AG  =  -  TAS  +  AE  -  (/rd)  A/%AV  -  (yV  -yV)/p 


where  AS  ,  AE  and  AV  are  the  activation  entropy,  activation  energy  and 

activation  volume  respectively;  d  is  a  dimensional  parameter  depending  on  the 

* 

structure  of  the  crack  tip.  The  last  term  is  included  following  the  approach 

of  Charles  and  Hillig  [1962]  and  accounts  for  changes  in  the  surface  energy  of 

the  crack  tip  during  the  reaction.  The  effect  of  stress  on  the  crack  tip 

1/2  ^ 

reaction  rate  is  included  in  the  (nd)  /£KjAV  term,  d  must  have  units  of 
2  ~l/2 

(dimension)  so  that  d  Kj  has  units  of  stress.  The  use  of  the  activation 
volume  in  this  context  assumes  that  the  tensile  stress  dependence  of  the 
reaction  rate  can  be  expressed  as  the  negative  of  the  pressure  dependence. 
Equations  3  and  5  can  be  combined  to  yield 


V  =  Vo  exp  [(-E  +  b  Kj)/RT] 


where  all  of  the  non-stress  dependent  terms  are  included  in  E  . 


Wiederhorn  et  al.  (1982)  use  the  crack  tip  radius,  p,  but  for  values  of  p 
approximately  that  of  the  network  spacing  in  the  solid,  this  term  probably  has 
little  meaning. 


The  previous  derivation  of  the  kinetics  of  crack  growth  was  completely 
general.  No  details  of  the  chemical  reaction  between  the  environment  and  the 
glass  were  discussed.  One  important  question  is  why  H^O  is  such  an  effective 
crack  growth  agent  for  all  silicate  glasses,  when  the  corrosion  rate  of  a 
material  such  as  vitreous  is  extremely  low? 


Michalske  and  Freiman  [1983]  have  described  a  specific  chemical  mechanism 
by  which  strained  Si-0  bonds  in  vitreous  silica  react  with  molecules  of  a  gas 
or  liquid.  This  model  for  the  HgO-SiOg  reaction  is  shown  in  Figure  3  as  a 
three  step  process: 

Step  1.  An  H^O  molecule  orients  itself  with  respect  to  an  Si-O-Si  bond 
at  the  crack  tip  such  that  the  lone  electron  pair  orbitals  on 
the  oxygen  of  the  water  molecule  are  aligned  toward  the 
silicon,  and  hydrogen  bonding  occurs  between  the  0  and 

the  hydrogen.  The  strain  (as  much  as  20%)  on  the  bridging  S1-0 
bond  clearly  enhances  the  tendency  to  react  at  this  site. 


Step  2.  Electron  transfer  from  the  0water  to  the  SI  occurs  simultaneous 
with  proton  transfer  to  the  0sjnca-  Two  new  bonds  are  formed, 

Si"°water  and  H'°s11ica- 


Step  3.  Rupture  of  the  weak  hydrogen  bond  between  0  and  the 

transferred  hydrogen  occurs  to  yield  Si-OH  surface  groups  on 
each  fracture  surface. 
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Note  that  the  model  does  not  require  prior  dissociation  of  the  water  molecule, 
nor  must  reaction  products  be  removed  from  the  solid.  Wiederhorn  and  Johnson 
[1973]  observed  definite  effects  of  hydrogen  ion  activity  on  crack  growth  in 
silica  and  other  glasses  suggesting  that  OH  and  H+  ions  may  show  a  different 
stress  dependent  reaction  rate  than  H2O.  These  effects  need  further  study, 
however. 

The  above  model  suggests  that  other  environments  should  also  enhance 
crack  growth  in  silica  if  the  species  has  structural  and  bonding  features 
similar  to  water,  namely,  proton  donor  sites  on  one  part  of  the  molecule,  and 
lone  pair  orbitals  on  another.  It  is  likely  that  there  are  size  limitations 

o 

as  well  since  the  Si-0  bond  distance  is  only  1.63  A.  It  is  interesting  to 
note  the  similarity  of  this  model  to  that  of  Griggs  and  Blacic  [1965]  for  the 
hydrolytic  weakening  of  quartz. 

Figure  4  shows  crack  growth  curves  for  amorphous  silica  in  water, 
ammonia  (NHj),  hydrazine  (NgH^),  formamide  (HCONHg)  and  N2  gas.  As  will  be 
discussed  later,  the  plateau  in  the  N2  gas  curve  occurs  when  crack  growth  is 
controlled  by  the  rate  of  diffusion  to  the  crack  tip  of  the  small  quantity  of 
dissolved  water.  The  absence  of  this  plateau  in  the  curves  for  ammonia, 
hydrazine  and  formamide  leads  to  the  conclusion  that  their  reaction  with  the 
crack  tip  bonds  governs  crack  growth  rates  in  preference  to  the  reaction  of 
the  water  present  in  these  environments.  Based  on  the  plateau  in  the  crack 
growth  curves  obtained  in  environments  such  as  carbon  monoxide  (CO), 
acetonitrile  (CHjCN)  and  nitrobenzene  (CgHgNOg),  these  did  not  control  crack 
growth,  as  would  be  predicted  by  the  above  model,  since  each  of  these 


molecules  does  not  meet  all  of  the  required  bonding  specifications.  Although 
each  contains  lone  pair  orbitals,  none  can  donate  protons.  It  should  also  be 
noted  that  no  direct  correlations  were  found  to  exist  between  dielectric 
constant  or  dipole  moment  alone  and  an  environment's  ability  to  cause  crack 
growth  in  preference  to  water. 

Recent  crack  growth  data  obtained  on  a  commercial  soda-lime-silica 
glass  shows  the  same  trends  as  that  described  for  silica  [Freiman  and 
White,  1982].  There  was  no  obvious  participation  of  the  Na+  ions  in  the 
crack  growth  process.  However,  preliminary  data  suggest  that  when  the  Na 
content  is  increased  significantly  as  in  a  binary  33%  Na  -  67%  Si02  glass, 
variations  in  the  crack  growth  mechanisms  appear. 

Regions  IA  and  IB. 

The  mechanism  of  crack  growth  in  the  steeper  IA  regime  in  a  number  of 
glass  seems  to  be  similar  to  that  in  region  I,  but  there  is  some  evidence  that 
there  are  accompanying  changes  in  crack  tip  geometry.  Work  by  Michalske 
[1982]  has  demonstrated  the  occurrence  of  crack  blunting  and  a  stress 
corrosion  limit  in  soda-lime-silica  glass  in  water,  but  except  for  some  data 
in  borosilicate  glasses  [Wiederhorn  and  Johnson,  1973,  Simmons  and  Freiman, 
1980],  determination  of  a  definite  crack  limit  have  not  been  extended  to  other 
glasses  or  ceramics.  For  instance,  no  evidence  for  a  region  of  increased 
slope  has  been  observed  for  vitreous  SiOg  at  velocities  as  low  as  10  **  m/sec. 

A  Region  IB  plateau  in  the  velocity  range  of  10  to  10  m/sec  has  been 


observed  only  in  binary  alkali-silicate  glasses  tested  in  aqueous  solutions 
[Simmons  and  Freiman,  1981].  The  authors  hypothesize  that,  because  of  the 


great  solubility  of  these  glasses  in  water  at  lower  Kj' s,  cracks  grow  at  a 
rate  approximating  that  of  the  Si02  dissolution  rate.  A  very  small  stress  is 
apparently  sufficient  to  remove  corrosion  products  at  the  crack  tip  and 
accelerate  the  dissolution  process  there,  in  preference  to  the  sides  of  the 
crack.  However,  effects  of  stresses  generated  due  to  a  H+  for  Na+  exchange 
cannot  be  ruled  out. 


Region  II. 

This  plateau  appears  when  crack  growth  occurs  in  an  environment  in  which 
the  minor  constituent  in  the  liquid  or  gas  controls  the  crack  tip  reaction. 

As  the  crack  proceeds,  the  species  in  the  environment  which  reacts  with  the 
chemical  bonds  at  the  crack  tip,  e.g.  H20,  is  depleted  in  the  vicinity  of  the 
tip,  thereby  creating  a  concentration  gradient.  As  noted  earlier,  in  Region  I 
the  rate  is  controlled  by  the  reaction  rate  at  the  tip,  but  as  the  crack 
velocity  increases,  the  size  of  the  H20  depleted  zone,  6,  grows.  At  some 
crack  velocity  the  rate  at  which  the  active  species  can  diffuse  through  this 
zone  to  the  tip  becomes  slower  than  the  reaction  rate,  and  hence  becomes  the 
controlling  step  in  determining  rates  of  growth.  Based  on  the  use  of  the 
Stokes-Einstein  expression  for  the  diffusivity  of  a  molecule  in  a  liquid,  it 
was  shown  that  this  velocity  is  given  by:  [Wiederhorn  et  al.,  1982] 


V  -  0-0275  kT I  _o 

P'ateau~  6  nr  6  Uj 


where  k  is  the  Boltzman  constant,  T  Is  temperature,  r  is  the  size  of  the 
diffusing  species,  CQ  is  the  bulk  concentration  of  this  diffusing  species,  and 
q  is  the  viscosity  of  the  solution.  Experimental  data  for  alcohol-water 


solutions  has  verified  the  applicability  of  Equation  (7)  (Figure  5).  Because 
diffusion  can  also  take  place  in  from  the  faces  of  a  double  cantilever  beam 
specimen  for  instance,  the  Region  II  velocity  is  also  crack  front  length 
dependent  and  can  show  a  small  stress  intensity  dependence,  i.e.  a  slope  =  2 
to  6. 


Region  III. 

Crack  growth  curves  in  Region  III  are  quite  steep,  n  >  100.  Wiederhorn 
[1967]  showed  that  crack  growth  rates  in  Region  III  for  soda-lime  glass  are 
independent  of  water  concentration.  Freiman  [1974],  showed  that  while  crack 
growth  rates  in  Region  III  are  not  affected  by  water  in  the  environment,  in 
agreement  with  Wiederhorn,  there  are  effects  of  the  liquid  in  which  the  water 
was  dissolved,  e.g.  the  chain  length  of  an  alcohol.  Similar  data  was 
subsequently  reported  by  Richter  [1977].  A  model  explaining  Region  III  crack 
growth  based  on  an  electrostrictive  interaction  between  the  solid  and  the 
fluid  at  the  crack  tip  has  recently  been  proposed  by  Wiederhorn  et  al.  [1982]. 
Their  model  is  synopsized  as  follows: 


If  we  assume  that  the  slope  of  the  V-Kj  curve  in  Region  III  (or  in  Region 

I)  is  determined  by  the  stress  dependence  of  a  reaction  rate  at  the  crack  tip, 

* 

then  the  activation  volume,  AV  ,  can  be  taken  as  a  measure  of  the  stress 

4 * 

dependent  process.  From  a  physical  point  of  view  AV  represents  the 
difference  in  volume  between  the  reactants  and  the  activated  complex.  It  is 
assumed  that  the  formation  of  charges  during  the  rupture  of  a  strained  Si-0 
bond  at  crack  tip  has  a  significant  effect  on  AV  ,  suggesting  that 
environments  which  can  alter  the  magnitude  of  these  charges  will  influence  AV 
and  therefore  the  slope  of  the  V-Kj  curve. 


Wiederhorn  et  al.  derived  the  following  expression  for  AV  in  a 
solid-solvent  system  in  which  the  electrostatic  contributions  are  of  primary 
importance: 


N  2  ,  ,  „  . 

AVr,  =  (3B)_i  (e.  +  zXL  -  (e,  +  cX*  ^2 


'ES 


(8) 
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where  N  is  Avogadro's  number,  e  is  the  unit  of  charge,  r  is  the  atomic 

radius,  (contributions  due  to  Si  and  0  are  calculated  separately)  B  is  the 

bulk  modulus  of  the  solid  and  and  t 2  are  the  dielectric  constants  of  the 

liquid  and  solid  respectively.  One  can  see  that  the  second  term  on  the  right 

contains  the  partial  derivative  of  the  solid  dielectric  constant  with  respect 

to  presssure,  i.e.  the  stress  contribution.  A  plot  of  log  (velocity)  versus 

Kj,  should  have  a  slope  proportional  to  AV.  Figure  6  is  a  plot  of  AV^ 

calculated  from  Equation  (8)  compared  to  the  determined  from  the  slope 

* 

of  the  V-Kj  curve.  The  correlation  between  the  calculated  and  measured 
values  suggests  the  general  validity  of  the  above  approach. 


EFFECTS  OF  CHEMICAL  BONDING  AND  MICROSTRUCTURE 

Most  ceramic  materials  are  susceptible  to  crack  growth  in  the  presence  of 
water.  What  is  surprising  is  that  these  materials  include  such  diverse 
chemical  compositions  as  A1203,  MgF2,  ZnSe,  and  graphite,  as  well  as  many 
others.  Is  the  mechanism  of  subcritical  crack  growth  the  same  as  that 


A  value  of  d  of  0.5  nm  was  assumed  for  this  calculation.  While  there  is 


clearly  uncertainty  in  this  number,  it  should  be  correct  within  at  least  an 
order  of  magnitude. 


described  for  Si02?  What  are  the  effects  of  changing  the  ratio  of  covalent  to 
ionic  bonding  in  the  solid?  What  are  the  effects  of  strain  on  the  bond 
character?  Complete  answers  to  these  questions  are  not  known  at  present. 
However,  recent  data  obtained  on  Al^  (sapphire)  [Michalske  et  al.,  1982] 
suggests  that  this  material  behaves  identically  to  amorphous  Si02  in  that  the 
same  environments  that  governed  crack  growth  in  the  latter  also  control  in 
A1203  and  to  the  same  relative  degree.  In  the  more  ionically  bonded  MgF2, 
however,  the  above  authors  found  that  ammonia  no  longer  controls  crack  growth 
in  Region  I,  while  acetonitrile  (CHgCN)  is  an  effective  crack  growth 
environment  even  though  it  has  no  proton  donating  capability.  This  behavior 
suggests  that  individual  atom  solvation  may  play  a  part  in  the  fracture  of 
ionic  solids.  There  may  be  a  pure  electrostatic  bonding  of  the  solvent 
molecules  to  the  Mg  and  F  atoms.  The  seeming  lack  of  any  measurable 
subcritical  crack  growth  in  silicon,  a  completly  covalent  solid,  indicates  the 
needs  for  some  degree  of  charge  separation  during  stressing. 

Very  little  detailed  crack  growth  data  exists  for  polycrystalline 
ceramics.  There  are  several  general  observations  that  can  be  made,  however. 
First,  it  is  apparent  that  grain  boundary  structure  and  chemistry  becomes  a 
significant  factor  in  determing  crack  growth  rates.  One  can  show,  for 
instance,  that  MgO,  which  does  not  appear  to  be  subject  to  subcritical  crack 
growth  in  single  crystal  form,  does  undergo  delayed  failure  when  prepared  in 
polycrystalline  form  [Freiman,  1976],  Fracture  occurs  almost  exclusively 
along  grain  boundaries.  Other  polycrystalline  ceramics  in  which  subcritical 
crack  growth  is  observed  also  exhibit  primarily  intergranular  failure  over  the 
regime  in  which  the  crack  is  growing  slowly,  even  in  cases  where  these  is  no 
obvious  grain  boundary  phase  [Mecholsky  and  Freiman,  1981].  A  transition  from 


an  intergranular  to  a  transgranular  failure  mode  is  frequently  observed  at  the 
critical  flaw  boundary  in  these  materials.  Whether  this  kind  of  behavior  is 
due  to  chemical  effects  or  to  changes  in  the  ability  of  the  crack  front  to 
change  directions  at  different  velocities  is  not  known.  In  addition,  there 
are  effects  of  grain  size  on  the  slope  of  the  V-Kj  curve  which  are  not 
understood. 

SUMMARY 

An  overall  picture  of  the  complex  relationship  between  crack  velocity  and 
stress  intensity  factor  has  been  presented.  It  has  been  shown  that  chemical 
reaction  kinetics  can  be  used  to  explain  the  dependence  of  crack  growth  rates 
on  the  activity  of  a  given  environment.  Mechanisms  of  subcritical  crack 
growth  in  glasses  and  ceramics  have  been  described  in  terms  of  models  for  the 
stress  induced  reactions  of  environments  with  crack  tip  bonds.  It  has  been 
shown  that  the  tendency  of  a  given  substance  to  participate  in  such  a  reaction 
can  be  predicted  from  knowledge  of  the  molecular  chemistry  of  the  environment 
and  the  bonding  in  the  solid.  While  only  limited  crack  growth  data  exist  for 
polycrystalline  ceramics,  it  is  clear  that  processes  in  these  materials  are 
dominated  by  grain  boundary  chemistry  and  structure. 

In  closing,  it  is  perhaps  worthwhile  to  point  out  that  no  mention  was 
made  in  this  paper  of  the  possible  effects  of  zeta  potential  on  crack  growth 
even  though  other  authors  in  this  conference  discuss  it  rather  extensively, 
i.e.  Westwood  [1983],  Dunning  [1983].  It  is  our  observation  that  no  clear 
correlation  has  been  made  between  zeta  or  surface  potential  and  the  rates  of 
subcritical  crack  growth.  There  is  serious  question  whether  zeta  potentials 


measured  on  bulk  solids  or  powders  will  be  similar  in  any  way  to  those  within 
the  small  confines  of  a  crack  tip.  In  addition,  these  measurements  suggest 
what  has  happened  to  a  surface  after  reaction  with  an  environment,  while  we 
are  interested  in  the  kinetics  of  the  interaction  process. 
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Figure  1 


Schematic  of  the  different  regions  of  a  typical  crack  growth  curve 
for  glasses  and  ceramics  tested  in  air  for  instance. 
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Figure  2.  Crack  velocity  in  soda-lime  silica  glass  at  a  Kj  -  0.45  MPaM1^2  as  a 
function  of  the  relative  humidity  i.e.  activity  of  water,  in  various 
straight  chain  alcohols.  (After  Freiman  [1974]). 
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Figure  4.  Crack  growth  curves  for  amorphous  silica  tested  in  environments 
whose  molecular  structure  should  be  favorable  to  cause  stress 
induced  reactions.  (After  Michalske  and  Freiman  [1982]). 
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Figure  6. 


Activation  volume,  AVtntal >  determined  from  the  slope  of  the  V-Kj 
curve  for  soda-lime  silica  glass  tested  in  environments  whose 
dielectric  constant  varies  from  1  (vacuum)  to  78  (H2O  plotted 

ePaiSt  [VlMlfP  ca1culated  from  Ration  8.  (After  Wiederhorn 
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Effects  of  Water  and  other  Dielectrics  on  Crack  Growth 

S.  M.  Wiederhorn,  S.  W.  Freiman,  E.  R.  Fuller,  Jr. 
and  C.  J.  Simmons 
National  Bureau  of  Standards 
Washington,  D.C.  20234 

ABSTRACT 

The  effect  of  water  and  a  variety  of  organic  liquids  on  the  crack 
growth  rate  in  soda  lime  silica  glass  was  investigated.  When  water  is 
present  in  organic  liquids,  it  is  usually  the  principal  agent  that 
promotes  subcritical  crack  growth  in  glass.  In  region  I,  subcritical 
crack  growth  is  controlled  primarily  by  the  chemical  potential  of  the 
water  in  the  liquid;  whereas  in  region  II,  crack  growth  is  controlled 
by  the  concentration  of  water  and  the  viscosity  of  the  solution  formed 
by  the  water  and  the  organic  liquid.  In  region  III,  where  water  does 
not  affect  crack  growth,  the  slope  of  the  crack  growth  curves  can  be 
correlated  with  the  dielectric  constant  of  the  liquid.  It  is  suggested 
that  these  latter  results  can  be  explained  by  electrostatic  interactions 
between  the  environment  and  charges  that  form  during  the  rupture  of  Si-0 
bonds . 


Effects  of  Water  and  other  Dielectric  Fluids  on  Crack  Growth 

S.  M.  Wiederhorn,  S.  W.  Frelman,  E.  R.  Fuller,  Jr. 
and  C.  J.  Simoons 
National  Bureau  of  Standards 
Washington,  D.C.  20234 

1.  Introduction: 

Subcritical  crack  growth  in  glasses  is  believed  to  result  primarily 
from  a  stress-enhanced  chemical  reaction  between  a  chemical  environment, 
usually  water,  and  the  Si-0  bonds  of  the  glass.  The  reaction  occurs 
preferentially  at  crack  tips  where  stresses  are  high.  Prior  investiga¬ 
tions  of  subcritical  crack  growth  in  glass  have  indicated  a  rather  complex 
dependence  of  crack  velocity  on  the  applied  stress  Intensity  factor  and 
on  the  amount  of  water  in  the  environment.  If  crack  velocity  is  plotted 
as  either  an  exponential,  or  power  function  of  applied  stress  intensity 
factor,  a  trlmodal  curve  is  obtained,  Figure  1;  each  of  the  three  regions 
of  the  curve  results  from  a  different  mechanism  of  crack  growth.1'2 

Region  I  crack  growth  behavior  is  characterized  by  a  dependence  of 
crack  velocity  on  stress  intensity  factor,  K^,  and  on  the  partial 
pressure  of  water  in  the  environment.1  If  either  of  these  two  quantities 
are  increased,  the  crack  velocity  is  observed  to  increase.  Crack  growth 
studies  conducted  in  either  nitrogen  gas1  or  in  alcohols3'4  suggest  that 


the  crack  velocity  is  reaction  rate  limited  in  region  1,  and  therefore, 
is  controlled  by  the  chemical  potential  of  water  in  the  environment  at 
any  level  of  applied  stress  intensity.  Hence,  identical  crack  growth 
behavior  should  be  obtained  in  environments  in  which  the  chemical 
potential  of  water  is  the  same. 

Region  II  crack  growth  behavior  is  also  characterized  by  a 
dependence  of  crack  velocity  on  the  amount  of  water  in  the  environment. 

In  contrast  to  region  1  behavior,  however,  the  crack  velocity  in  region 
II  does  not  depend  strongly  on  the  applied  stress  intensity  factor. 
Theoretical  interpretations  of  region  II  crack  growth  behavior  lead  to 
the  conclusion  that  the  crack  velocity  is  controlled  by  the  rate  of 
transport  of  water  from  the  environment  to  the  crack  tip.1*2’5  Although 
it  is  recognized  that  crack  growth  in  region  II  is  limited  by  factors  that 
affect  transport  processes,  experimental  data  to  quantify  this  relation 
have  not  yet  been  obtained. 

Region  III  crack  growth  behavior  is  characterized  by  a  strong 
dependence  of  crack  growth  rate  on  stress  intensity  factor,  and  by  a 
complete  absence  of  any  dependence  of  crack  growth  rate  on  water  in  the 
environment.  Of  the  three  regions  of  crack  growth,  region  III  behavior 
is  the  least  understood.  Region  III  crack  growth  has  been  observed  in 
vacuum  for  a  variety  of  commercial  silicate  glasses,  but  not  for  glasses 
that  have  a  high  concentration  (>90  percent)  of  silica.6  Because  of 
these  observations  it  has  been  suggested  that  region  III  crack  growth  is 
controlled  by  fracture  mechanisms  that  do  not  depend  on  environment. 


Studies  of  region  III  crack  growth  in  normal  alcohols,  however,  conflict 
with  this  view.3*1*  In  normal  alcohols,  it  has  been  shown  by  Freiman3 
that  the  slope  of  the  crack  growth  curves  is  not  as  steep  as  that 
obtained  in  vacuum,  or  in  dry  nitrogen.  Furthermore  Freiman3  and 
Richter4  have  demonstrated  that  the  positions  of  the  crack  growth  curves 
depend  on  the  chain  length  of  the  alcohol.  Hence,  in  some  cases  an 
environmental  dependence  of  region  III  crack  growth  is  apparent. 

This  paper  provides  additional  experimental  data  to  clarify  our 
understanding  of  subcrltical  crack  growth  in  glass.  Crack  growth  data 
obtained  in  alcohols  confirm  the  proposed  role  of  chemical  potential  in 
controlling  the  crack  growth  rate  in  region  I.  In  region  II,  however, 
water  concentration,  rather  than  chemical  potential,  will  be  shown  to 
control  the  rate  of  crack  growth.  A  model  for  diffusion  controlled  crack 
growth  is  developed  to  show  that  fluid  viscosity  is  also  a  parameter  of 
importance.  Finally,  data  is  presented  to  confirm  the  fact  that  crack 
growth  behavior  in  region  III  is  affected  by  the  chemical  composition  of 
the  test  media.  Crack  growth  in  region  III  is  influenced  by  an  electro¬ 
static  interaction  between  the  environment  and  charges  that  develop  at 
the  crack  tip  during  the  fracture  process. 

2.  Experimental  Procedure 

Crack  propagation  studies  were  conducted  on  soda  lime  silica 
glass  in  normal  alcohols,  heptane,  formamide,  trichloroethylene  and 
acetonitrile.  The  alcohols  were  selected  as  test  media  because  the 
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solubility  of  water  in  alcohols  decreases  with  increasing  alcohol  chain 
length,  so  that  it  is  possible  to  produce  solutions  of  water  and  alcohol 
that  have  equal  chemical  potentials  even  though  they  contain  different 
concentrations  of  water.  These  solutions  were  used  to  elucidate  the 
effect  of  water  concentration  on  region  II  crack  growth  behavior.  Where 
possible,  the  partial  pressure  of  water  in  the  solutions  was  determined 
by  direct  analysis  with  an  immersible  hygrometer.  The  amount  of  water 
in  solution,  (i.e.  the  concentration)  was  determined  by  the  Karl-Fisher 
technique.  This  procedure  permitted  us  to  make  independent  comparisons 
of  effects  of  chemical  composition  and  chemical  potential  on  crack  growth 
in  region  II.  A  set  of  ethyl  and  butyl  alcohol  solutions,  used  primarily 
for  the  study  of  region  I  crack  growth,  was  made  by  adding  predetermined 
amounts  of  water  to  the  "pure"  alcohols.  The  chemical  potential  of 
these  solutions  were  then  determined  by  reference  to  phase  equilibrium 
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diagrams  for  these  solutions.  The  viscosities  of  solutions  used 
for  region  II  crack  growth  studies  were  determined  by  direct  analysis 
using  a  Cannon-Fenske  viscosimeter  (Table  1). 

Crack  growth  studies  in  region  I,  II  and  III  were  conducted  on  soda- 
lime-silica  glass  microscope  slides  using  both  the  double  cantilever  beam 
technique  and  the  double  torsion  technique.  9*10  Region  I  studies  were 
conducted  using  double  cantilever  beam  specimens,  7  x  25  x  2  mm,  with  a 
0.5  mm  deep  groove  along  the  midplane  on  both  sides  of  the  specimen  to 
control  the  direction  of  crack  growth.  Double  torsion  specimens, 
(ungrooved)  75  x  25  x  1  mm,  were  used  for  crack  growth  measurements  in 
regions  I,  II  and  III.  In  regions  I  and  II  the  load  relaxation  method 
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was  used  to  collect  crack  growth  data,  whereas  in  region  III,  both  load 
relaxation  and  constant  displacement  techniques  were  used. 

Studies  in  region  II  were  made  in  saturated  alcohols  prepared  by 
mixing  excess  water  with  the  alcohols  and  agitating  the  liquids  several 
times  a  week  for  six  months.  A  chemical  analysis  of  the  vater  in 
solution  and  a  measurement  of  the  partial  pressure  of  the  water  in  the 
alcohol  are  presented  in  Table  1.  With  the  exception  of  the  octyl 
alcohol,  all  of  the  alcohols  were  found  to  be  saturated  vlth  water.  We 
were  somewhat  surprised  by  the  observation  that  the  octyl  alcohol  used 
in  our  study  was  not  fully  saturated,  and  have  no  explanation  for  this 
result.  However,  crack  growth  studies  (region  II)  to  be  discussed 
below  are  consistent  with  this  observation. 

The  double  torsion  technique  was  also  used  to . investigate  crack 
growth  behavior  in  region  III.  To  determine  if  the  position  of  the 
crack  velocity  curve  in  alcohol  differed  significantly  from  that  in  dry 
nitrogen,  crack  propagation  data  were  obtained  by  the  relaxation  method 
on  the  same  specimen,  first  in  dry  nitrogen,  then  in  alcohol.  Three 
crack  velocity  curves  were  first  obtained  in  nitrogen  gas.  Without 
disturbing  the  specimen,  either  methyl,  amyl,  or  decyl  alcohol  were 


added  to  the  test  chamber  and  three  more  crack  velocity  curves  were 
obtained.  This  procedure  enabled  us  to  reduce  systematic  variations 
and  thus  obtain  crack  propagation  curves  that  accurately  reflect  the 


positions  of  the  crack  growth  curves  in  alcohol,  relative  to  those  in 
dry  nitrogen  gas. 


In  a  second  study  of  region  III  crack  growth,  the  constant  displace1 
ment  rate  technique  was  used  to  obtain  crack  growth  data  in  dry  nitrogen 
gas,  methyl  alcohol,  decyl  alcohol,  acetonitrile  and  heptane.  At  least 
10  experimental  runs  were  made  at  each  displacement  rate  in  order  to 
reduce  systematic  error  that  resulted  in  variations  in  measured  K^.  The 
mean  value  and  standard  deviation  of  were  obtained  for  each  displace¬ 
ment  rate  measurement,  and  at  least  three  displacement  rates  were  used 
to  determine  the  position  and  slope  of  each  crack  velocity  curve. 
Displacement  rates  used  in  these  studies  were  selected  so  that  the 
velocities  were  in  region  III.  The  velocities  had  to  be  low  enough, 
however,  to  assure  that  cavitation  of  the  fluid  at  the  crack  tip  did  not 
occur.  11  Some  experimental  results  in  region  III  were  also  obtained  on 
vitreous  silica  (C7940)  in  octyl  alcohol  using  the  applied  moment, 
double  cantilever  beam  technique.  12  Since  silica  glass  normally  does 
not  exhibit  subcrltical  crack  growth  in  region  III,*  the  demonstration 
of  subcrltical  crack  growth  in  octyl  alcohol  is  a  clear  Indication  of 
an  interaction  between  glass  and  alcohol. 

3.  Results 

The  crack  propagation  data  for  soda-lime-silica  glass  in  ethyl  and 
butyl  alcohols  are  shown  in  Figures  2,3.  The  shapes  and  positions  of 
these  curves  are  similar  to  those  obtained  earlier  on  this  glass  tested 
in  nitrogen  gas1  and  in  other  alcohols3’4.  In  regions  I  and  II,  crack 
velocity  is  observed  to  decrease  as  the  amount  of  water  in  the  alcohols 
is  decreased.  In  region  III,  the  data  for  the  butyl  alcohol  appears  to 


approach  a  conanon  curve,  as  it  does  for  nitrogen  gas1.  Hence  the 
mechanisms  for  fracture  of  soda  lime  silica  glass  in  ethyl  and  butyl 
alcohols  are  probably  identical  to  those  observed  in  nitrogen  gas  and  in 
long  chain  alcohols  studied  earlier  by  Freiman3  and  Richter4. 

The  crack  growth  data  from  region  I  of  the  present  paper  (Figures 
2  and  3)  can  be  compared  with  earlier  data  collected  in  nitrogen  gas1 
by  plotting  the  logarithm  of  the  crack  velocity  as  a  function  of  the 
logarithm  of  the  relative  humidity  of  the  water  in  the  alcohol.  To  do 
this,  crack  velocities  were  taken  from  Figures  2  and  3  for  an  applied 

L 

stress  intensity  factor  of  0.563  MPa-m  ,  which  is  the  same  value  of 

« 

used  in  the  earlier  study.  The  relative  humidities  were  determined  from 
the  concentrations  of  water  in  the  alcohols,  using  experimental  phase 
diagrams7*8 .  The  resultant  plot.  Figure  U,  indicates  a  relationship 
between  crack  velocity  and  relative  humidity  that  is  roughly  linear 
over  the  range  of  relative  humidities  studied.  Furthermore ,  the  alcohol 
data  are  found  to  lie  relatively  close  to  the  earlier  data  collected  in 
nitrogen  gas.  Since  the  scatter  of  crack  growth  curves  obtained  by  the 
double  cantilever  beam  technique  can  be  as  much  as  a  factor  of  *v2  along 
the  velocity  axis,  the  difference  between  the  three  sets  of  data  can  be 
attributed  to  experimental  scatter.  A  comparison  of  the 

* 

This  procedure  necessitated  a  linear  extrapolation  of  the  region  I  v-Kj 
data  for  the  butyl  alcohol  (Figure  3)  slightly  beyond  the  range  of  the 
data  in  region  I. 


ethyl  and  butyl  alcohol  data  with  crack  growth  data  collected  by 
Freiman3  on  long  chain  alcohols  was  also  made;  agreement  between  the 
two  sets  of  data  was  similar  to  that  shown  in  Figure  4.  These  results 
confirm  earlier  conclusion  that  crack  growth  in  region  1  is  controlled 
by  the  relative  humidity  (i.e.  partial  pressure)  of  the  environment  in 
contact  with  the  crack  tip.1'3  As  will  be  discussed  in  section  4.1,  such 
behavior  is  consistent  with  reaction  rate  theory,  in  which  the  effect  of 
environment  on  the  crack  growth  rate  depends  on  the  chemical  potential 
or  activity  of  the  water  in  the  alcohol. 

In  contrast  to  the  results  presented  above,  crack  growth  in  region 
11  does  not  appear  to  be  controlled  directly  by  the  relative  humidity 
of  the  water  in  the  alcohol.  This  conclusion  is  supported  by  plotting 
the  region  II  crack  growth  data  presented  in  Figures  2  and  3  as  a 
function  of  relative  humidity  of  water  in  the  alcohol.  Data  plotted  in 
this  manner  are  compared  in  Figure  5  with  similar  data  obtained  by 
Freiman  from  studies  of  crack  growth  in  octyl  alcohol3.  To  assure 
ourselves  that  the  two  sets  of  data  were  consistent,  new  crack  growth 
data  were  collected  in  nominally  saturated  octyl  alcohol.  The  crack 
velocity  in  region  11  for  the  octyl  alcohol  was  found  to  be  in  line  with 
the  earlier  data  by  Freiman.  A  comparison  of  the  data  collected  in  ethyl 
alcohol  and  butyl  alcohol  with  the  data  collected  in  octyl  alcohol  shows 
that  the  ethyl  and  butyl  alcohol  data  lie  far  from  the  curve  obtained 
for  the  octyl  alcohol.  Since  the  separation  of  the  data  is  much  greater 
than  can  be  explained  by  experimental  scatter,  we  conclude  that  factors 
other  than  the  relative  humidity  of  water  control  crack  growth  in 
region  II. 


The  conclusion  stated  above  is  supported  by  the  double  torsion 
studies  (relaxation  method)  conducted  in  saturated  alcohols.  Figure  6. 
With  the  exception  of  the  octyl  alcohol,  which  was  not  fully  saturated, 
the  crack  velocity  in  region  II  is  observed  to  decrease  as  the  chain 
length  of  the  alcohol  is  increased.  This  decrease  in  velocity  correlates 
directly  with  the  concentration  of  water  in  these  alcohols,  and  inversely 
with  the  viscosity  of  the  alcohols.  Since  the  alcohols  used  to  obtain 
Figure  6  were  fully  saturated  (with  the  exception  of  the  octyl  alcohol), 
they  had  the  same  relative  humidity  of  water  as  the  water  with  which 
they  were  in  contact.  The  change  of  crack  velocity  of  approximately  1.5 
orders  of  magnitude  in  region  II,  therefore,  cannot  be  explained  by  a 
variation  in  the  relative  humidity  of  the  water  in  the  alcohol  solutions. 
A  model  will  be  presented  in  Section  4.3  of  this  paper,  which  explains 
region  II  crack  growth  behavior  in  terms  of  water  concentrations  and 
viscosity  of  the  water-alcohol  solution. 

The  data  obtained  in  region  III  by  the  double  torsion  technique 
(relaxation  method)  is  presented  in  Figure  7,  for  methyl,  amyl  and  decyl 
alcohol.  In  agreement  with  other  studies3*4,  our  data  shows  a  small, 
but  systematic  shift  of  the  crack  propagation  curve  to  higher  values  of 
as  the  length  of  the  alcohol  chain  is  increased.  Furthermore, 
regardless  of  alcohol  chain  length,  the  curves  obtained  in  alcohol  had 
slopes  that  were  less  than  those  obtained  in  dry  nitrogen  gas. 

The  most  surprising  aspect  of  the  data  shown  in  Figure  7  has  to  do 
with  the  position  of  the  crack  growth  curve  obtained  in  decyl  alcohol. 
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Three  separate  runs  using  this  alcohol  resulted  in  data  that  were 
similar  to  that  shown  in  Figure  7.  In  all  cases,  the  crack  growth  curve 
in  decyl  alcohol  crossed  the  curve  obtained  in  nitrogen  gas,  suggesting 
that  at  higher  crack  growth  rates,  decyl  alcohol  was  retarding  the  crack 
growth.  This  finding,  and  the  observation  of  a  systematic  shift  and  a 
lower  slope  of  the  crack  growth  data  in  the  alcohols  suggests  that  the 
alcohols  are  involved  in  some  intimate  way  with  the  fracture  process  in 
region  Ill.  To  further  check  the  possibility  of  an  intimate  interaction 
between  organic  fluids  and  the  silica  network  of  the  glass,  experimental 
data  were  collected  on  silica  glass  in  octyl  alcohol.  The  results  of 
this  study  are  shown  in  Figure  8.  If  silica  were  tested  in  nitrogen 
gas,  it  would  not  have  been  possible  to  obtain  crack  growth  data  in 
region  III  since  fracture  occurs  abruptly  in  vitreous  silica  once  the 
crack  has  passed  from  region  II. 6  The  fact  that  region  III  crack  growth 
was  obtained  In  octyl  alcohol  supports  our  contention  of  a  direct  effect 
of  the  alcohol  on  the  fracture  process  of  glass. 

To  substantiate  the  data  shown  in  Figure  7,  data  were  obtained 
using  the  double  torsion  technique,  constant  displacement  rate  method. 
The  data  obtained  in  this  way  (Figure  9)  support  the  data  shown  in 
Figure  7;  the  slope  of  the  curve  obtained  in  dry  nitrogen,  or  in  decyl 
alcohol  is  greater  than  that  obtained  in  methyl  alcohol.  Furthermore, 
the  data  for  decyl  alcohol  and  nitrogen  gas  fall  in  a  significantly 
higher  range  of  Kj  than  the  data  for  methyl  alcohol.  It  is  interesting 
to  note  that  the  data  for  decyl  alcohol  also  falls  in  a  higher  range  of 
Kx  than  that  for  nitrogen  gas,  a  finding  that  is  consistent  with  the 


data  in  Figure  7.  However,  because  of  the  experimental  scatter  in  both 
the  nitrogen  gas  and  decyl  alcohol  curves  this  difference  in  location 
cannot  be  construed  as  being  significant  from  this  figure  alone. 
Finally,  data  obtained  in  acetonitrile  also  had  a  slope  that  was 
significantly  lower  than  that  obtained  in  nitrogen  gas,  whereas  the 
slope  of  the  curve  obtained  in  heptane  was  similar  to  that  of  nitrogen 
gas. 

4.0  Discussion  of  Results 

4.1  Reaction  Rate  Theory 

The  results  presented  in  this  paper  can  be  understood  in  terms  of 
chemical  reaction  rate  theory,  which  has  been  used  recently  by 
Wiederhorn,  Fuller  and  Thomson*  3  to  derive  a  model  for  crack  growth  in 
ceramics  and  glasses.  To  ease  our  discussion  of  the  present  set  of 
data,  a  brief  review  of  the  model  is  presented  in  this  section,  giving 
the  underlying  theory  and  the  equations  that  are  pertinent  to  the 
results.  In  Section  4.4,  the  model  is  extended  to  account  for  electro- 
static  interactions  that  are  believed  to  occur  between  the  environment 
and  charges  that  form  during  rupture  of  Si-0  bonds  at  the  crack  tip. 

Consider  a  chemical  reaction  of  the  general  type: 


y.~  %'  •/  *w'  ^ 
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If  the  reverse  reaction  is  negligible,  the  rate  of  reaction  is  given  by: 


Rate  =  k  [A]a  [b]b  .  .  . 
r 


where  [A],  [B],  .  .  .  are  the  concentrations  of  the  reactants.  The 
order  of  the  reaction  with  regard  to  each  component  is  given  by  the 
exponents  in  equation  2  (i.e.  the  order  with  regard  to  A  is  a),  and  the 
overall  order  of  the  reaction  is  given  by  the  sum  a  +  b  +  .  .  ..  The 
reaction  rate  constant,  k^,  can  be  determined  by  transition  rate  theory 
and  is  given  by  the  following  relation:- 1*1**5 


kr  =  k  (kT/h)  exp  (-AG*  /RT)  (fAa  tb  .  .  ./f*) 


where  k  is  the  transmission  coefficient  and  is  usually  assumed  equal 
to  1;  AG^  is  the  Gibbs  free  energy  of  activation;  and  f^,  f^,  ...» 
f  are  the  activity  coefficients  for  the  reactants  and  for  the  activated 
state.  The  Gibbs  free  energy  of  activation  is  the  standard  Gibbs  free 
energy  change  associated  with  the  activation  process,  and  can  be 
related  to  the  standard  chemical  potentials  of  the  reactants  and  the 
activated  complex . 14*  *  1 5»  1 6 


._•*  *  a  b 

AG  =  y  -  v,  y_ 
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By  combining  equations  2  and  3,  and  by  using  the  definition  of 
chemical  activity  (a^  *  [A]  f^),  the  rate  of  chemical  reaction  can  be 


expressed  in  terms  of  the  chemical  activity  of  the  reactants: 


Rate  «  (ic/f*)  (kT/h)  a^  a^  .  .  .  exp  (-AG^/RT)  (5) 

where  a^,  a0,  .  .  .  are  the  chemical  activities  of  reactants  A,  B,  etc. 

Equation  S  can  also  be  expressed  in  terms  of  the  chemical  potential, 
l>A»  Mg.  •  ...  of  the  reactants: 

Rate  -  (ic/f*)  (kT/h)  exp  (ayA  +  by^  +  .  .  .  -  y^/RT  (6) 

where  the  chemical  potentials  of  the  reactants,  y^.  Mg.  etc.,  are  given 
by  the  following  type  of  relation:  yA  «  y^  +  RT  In  aA«  The  free  energy 
change  in  the  exponent  of  equation  6  includes  effects  of  chemical 
composition  on  reaction  rate.  By  contrast,  compositional  effects  on 
reaction  rate  are  included  in  the  activities  of  equation  5. 


For  purposes  of  this  paper,  the  crack  growth  rate  is  assumed  to  be 

proportional  to  the  rate  of  chemical  reaction  of  the  crack  tip,  and  the 

reaction  is  assumed  to  be  first  order  with  respect  to  the  water  in 
* 

solution.  For  these  conditions,  equations  5  or  6  can  be  written  as: 


The  assumption  of  a  first  order  reaction  may  not  be  valid  for  very  low 
activities  of  water.  In  an  earlier  study1,  it  was  noted  that  the  slight 
curvature  of  the  data  for  nitrogen  shown  in  Figure  4  suggested  a  change 
in  the  order  of  the  chemical  reaction  from  first  order  at  high  relative 
humidities  to  one-half  order  at  low  relative  humidities.  This  observa¬ 
tion  was  later  substantiated  for  long  chain  alcohols  by  Frelaan.3  At 
the  high  humidities  used  in  the  present  paper,  the  data  for  ethyl  and 
butyl  alcohols  can  be  represented  by  a  first  order  reaction  of  the 
water  in  solution  with  the  glass. 


V 


(7) 


vn  a  exp  (m  +  U  -  yj  )/RT 
o  s  so  g  o 


where  a  is  the  activity  of  the  water  in  solution,  and  where  y  ,  y 
s  so  g 

and  y*  refer  respectively  to  the  water  in  solution,  the  reactive 
species  in  the  glass  (presumably  Si-0  bonds)  and  the  activated  complex. 
The  standard  state  for  the  water  is  taken  as  pure  water  in  the  liquid 
state  at  1  atmosphere  pressure  and  a  temperature  of  25°  C,  a  convention 
that  follows  the  practice  used  in  references  7  and  8.  The  frequency 
factor,  kT/h,  the  transmission  coefficient,  k,  the  activity  coefficient, 
f*  and  geometric  factors  relating  the  number  of  bonds  broken  to  the 
crack  velocity  are  all  Incorporated  into  vq  of  equation  7.  It  Is  worth 
noting  that  the  dependence  of  crack  velocity  on  the  amount  of  water  in 
the  environment  is  all  contained  in  the  activity  a^,  whereas  the 
dependence  of  crack  velocity  on  stress  is  all  contained  within  the 
exponent  of  equation  7. 


Although  the  free  energy  of  activation  in  equation  7  contains  three 

terms  (AG*  -  y*  -  y„  -  y _ ) ,  these  terms  are  not  all  equivalent  with 

o  g  so 

regard  to  water  enhanced  crack  growth.  The  chemical  potentials  of  both 

the  glass,  y  ,  and  the  activated  state,  y*  depend  on  the  state  of 
g  o 

stress  at  the  crack  tip,  and  hence,  are  the  prime  determinants  of  the 

slope  of  the  v-K^.  curves  shown  in  Figures  1-3.  As  the  stress  applied 

to  the  crack  tip  is  Increased,  y*  -  y  must  decrease  to  account  for 

o  g 

an  increase  in  the  rate  of  crack  growth.  The  slope  of  the  v-K^.  curve 
depends  on  the  sensitivity  of  y*  -  y^  to  the  applied  stress,  the  greater 
the  sensitivity,  the  steeper  the  slope  of  the  v-Kj  plot.  By  contrast 


the  chemical  potential  of  the  water  does  not  depend  on  the  state  of  stress 
at  the  crack  tip,  and  hence,  has  no  effect  on  the  slope  of  the  v-K^  curve. 
It  is  for  this  reason  that  all  of  the  curves  In  region  I  of  Figures  1-3 
have  approximately  the  same  slope. 

A  quantitative  treatment  of  the  dependence  of  applied  stress  on  AG^ 

has  been  developed  by  modelling  the  crack  tip  as  an  elliptical  slit  in  a 

continuum.13  The  chemical  reaction  is  assumed  to  occur  at  the  tip  of 

the  slit  and  the  individual  bonds  at  the  crack  tip  are  assumed  to  be 

acted  upon  by  stresses  calculated  from  the  Inglia  relation.17  The 

following  expression  is  obtained  for  the  free  energy  of  activation:13 

AG*  -  -  TAS*  +  AE  *-[2  Kj/Oip)1*]  AV*  -  <Y*V*-  yV)/p  (8) 

±  i  i 

where  AS  ,  AE  T  and  AV  are,  respectively,  the  activation  entropy,  the 

activation  energy  and  the  activation  volume.  The  parameter  p  is  the 

radius  of  curvature  of  the  crack  tip,  and  is  assumed  to  be  of  the  order 

* 

of  the  lattice,  or  network  spacing  in  the  solid:  0.5  nm  for  glass. 

The  term  containing  K^.  on  the  right  hand  side  of  the  equation 
is  the  stress  to  which  the  reacting  species  at  the  crack  tip  are 
subjected.  This  portion  of  the  activation  free  energy  determines 
the  stress  dependence  of  crack  growth  rate.  More  specifically, 
if  crack  growth  data  are  fitted  to  the  following  empirical  relation: 

_ 

A  curved  crack  tip  is  assumed  primarily  to  estimate  the  stress  level 
at  the  crack  tip.  The  nature  and  structure  of  real  crack  tips  in 
glass  are,  at  present,  not  fully  understood. 


v  ■  vq  exp  (-E*  +  bKj)/RT,  then  the  slope  of  the  v-K.^  data,  b,  can  be 
used  to  determine  the  activation  volume  for  crack  growth.  Equating  b  to 


the  coefficient  of  In  equation  8,  the  following  relation  is  obtained 
for  AV^: 

AV*=  (b/2)  (up)1*  (9) 

This  expression  will  be  used  later  in  the  paper  to  relate  the 
crack  growth  data  (i.e.  the  slope  b)  to  a  variable  calculated  from  an 
electrostatic  model  of  environmental  Interactions  at  crack  tips. 


The  last  term  on  the  right  hand  side  of  equation  8  is  included  to 
account  for  the  fact  that  surface  curvature  can  modify  the  chemical 
potential  of  the  reacting  species.13  The  parameters  y  and  V  are  the 
surface  tension  and  molar  volume  of  the  reacting  species  (the  strained 
Si-0  bonds  of  glass  for  example)  before  reaction,  while  y  ^  and  V^are 
these  same  quantities  for  the  activated  complex  on  surface  during 
reaction.  The  inclusion  of  this  surface  energy  term  follows  the  approach 
first  taken  by  Charles  and  Hillig18  in  their  dissolution  model  of  static 
fatigue. 


4.2  Region  I  Crack  Growth 


The  crack  growth  data  collected  in  the  present  paper  can  best  be 

compared  with  equation  7  of  the  theory  presented  in  Section  4.1.  Since 

v  and  the  chemical  potentials  y  ,  y  and  y^  are  constants  for  a  given 
o  so  g  o 

temperature  and  state  of  stress,  the  equation  implies  that  at  a  constant 
applied  stress  intensity  factor,  the  crack  velocity  will  depend  only  on 


the  activity,  a  ,  of  water  in  the  alcohol  solution.  If  the  standard 
s 

state  of  the  chemical  potential  of  the  water  in  the  alcohol  is  taken  as 

pure  water,  then  the  chemical  activity  of  the  fully  saturated  alcohols 

used  in  the  present  studies  is  'V'l,  at  saturation,  (i.e.  the  chemical 

potentials  of  the  water  are  almost  the  same  for  all  of  the  alcohols). 

Furthermore ,  if  water  vapor  at  room  temperature  is  assumed  to  behave  as 

a  perfect  gas,  then  the  activity  of  the  water  Is  given  by  the  ratio  of 

the  partial  pressure,  pg,  of  water  over  the  solution  to  that,  P8Q»  over 

pure  water:  a  =  p  /p  .  16  Since  p  /p  is  defined  as  the  relative 
s  s  so  s  so 

humidity,  the  crack  velocity  is  predicted  to  be  proportional  to  the 
relative  humidity  of  the  water,  as  is  observed  experimentally  (Figure  U). 
Hence  we  conclude  that  the  chemical  activity,  or  what  is  equivalent,  the 
chemical  potential  of  water  in  solution  is  the  prime  driving  force  for 
crack  growth  in  region  I. 

The  above  observation  has  rather  broad  implications  with  regard  to 
the  effect  of  water  on  the  lifetime  of  structural  ceramics.  As  long  as 
the  chemical  potential  of  water  is  the  same  in  two  solutions,  the  response 
of  cracks  in  structural  ceramics  to  stress  vill  be  identical  with  regard 
to  region  I  crack  growth.  Since  the  lifetime  of  ceramics  is  controlled 
mainly  by  region  I  (i.e.  slow)  crack  growth,  the  chemical  potential  of 
water  in  solution  is  the  prime  environmental  determinant  of  structural 
reliability.  Although  several  investigations  of  the  effect  of  relative 
humidity  on  component  lifetime  are  consistent  with  this  prediction, 
this  subject  has  not  been  studied  systematically  and  certainly  warrants 
further  investigation. 
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4.3  Region  II  Crack  Growth 


The  mechanism  described  in  Section  4.1  will  control  crack  growth  as 
long  as  the  number  of  water  molecules  reaching  the  crack  tip  maintains 
the  activity  of  water  at  the  crack  tip  constant.  When  the  velocity  is 
fast  enough,  the  water  will  be  used  up  at  a  rate  faster  than  can  be 
supplied  by  diffusion  from  the  environment.  When  this  happens,  the  rate 
of  transport  will  control  the  rate  of  crack  motion. 

Using  a  one  dimensional  model  of  diffusive  transport  of  water  to  a 
crack  tip  it  has  been  shown  that  the  rate  of  crack  growth  in  region  II 
can  given  by:1 

v  -  0.0275  D  XQ/n6  (10) 

where  D  is  the  dlffusivity  of  water  in  the  alcohol,  6  la  the 

diffusive  boundary  layer  thickness  in  front  of  the  crack  tip,  Xq  is  the 

molar  concentration  of  the  water  in  the  solution  and  n  is  the  order  of 

* 

the  chemical  reaction.  For  relative  humidities  >0.1,  the  parameter  n 
has  been  shown  to  be  1  for  water  reactions  in  nitrogen  gas  and  will  be 
assumed  to  be  equal  to  1  for  the  alcohols. 


The  constant  0.0275  is  appropriate  when  all  variables  in  equation  10 
are  expressed  in  SI  units. 


The  diffusivity,  D,  of  water  in  alcohols  has  been  measured  in  only  a 
few  of  the  shorter  chain  length  alcohols.21  Hence  in  order  to  compare 
the  data  with  equation  4,  it  is  necessary  to  use  a  theoretical  expression 
for  the  diffusion  coefficient.  For  this  purpose  we  use  the  Stoke s- 
Einstein  relation:  22 

D  =  (kT)/(6irrn)  (ll) 

where  n  is  the  fluid  viscosity,  r  is  the  radius  of  the  diffusing 
molecule  and  k  is  the  Boltzman  constant.  Substituting  equation  11  Into 
equation  10,  the  following  equation  is  obtained  for  the  crack  velocity 
in  region  II: 

v  =  (0.0275  kT/(6irr6)  )(X  /n)  (12) 

o 

If  the  form  of  equation  12  is  correct,  then  the  crack  velocity  in 
region  II  should  be  proportional  to  Xo/n  (i.e.  the  ratio  of  the  molar 
concentration  of  water  in  the  alcohol  to  the  viscosity  of  the  alcohol). 
The  data  plotted  in  this  manner  (Figure  10)  scatters  about  a  line  of 
slope  1,  as  would  be  expected  from  the  theory.  Much  of  the  scatter  in 
Figure  10  can  be  attributed  to  the  uncertainty  in  selection  of  the 
characteristic  velocity  of  region  II.  As  can  be  seen  from  examination 
of  Figure  6,  the  error  in  estimating  the  crack  velocity  in  region  II  can 
be  as  much  as  a  factor  of  2.  Furthermore,  the  fact  that  crack  velocity 
in  region  II  exhibits  a  small  but  significant  dependence  on  applied 
stress  intensity  factor  contributes  to  the  uncertainty  in  selecting 


characteristic  velocities.  Considering  these  sources  of  scatter  in 
the  experimental  data,  the  correlation  between  the  variables  shown  in 
Figure  10  lends  credence  to  our  model  of  diffusion  controlled  crack 
growth . 

The  data  shown  in  Figure  10  can  be  used  to  obtain  an  estimate  of  the 
diffusive  boundary  layer,  6,  at  the  crack  tip.  The  intercept  on  Figure 
10  is  approximately  2.5  x  10~8  m/s;  substituting  this  value  of  crack 
velocity  into  equation  12  and  assuming  T  =  300  K  and  r  (the  molecular 
radius  of  H2O)  =  0.3  nm  a  value  of  fi  =  0.8  pm  is  obtained.  This  value 
is  approximately  one-third  that  obtained  in  earlier  studies  of  crack 
growth  in  nitrogen  gas.  Since  the  collision  distances  in  liquids 
(^.5  nm)  are  less  than  the  collision  distances  in  gases  (^50  nm  at 
standard  temperature  and  pressure),  the  boundary  layer,  6,  in  a  liquid 
is  expected  to  be  less  than  that  in  a  gas.  The  thickness  of  the  boundary 
is  likely  to  vary  from  one  type  of  liquid  to  another,  depending  on  the 
nature  of  the  interactions  that  occur  between  the  liquid  and  the  water 
that  is  dissolved  in  it. 

Region  III  Crack  Growth 

The  data  presented  in  this  paper  and  in  the  papers  by  Freiman3  and 
Richter4  suggest  that  alcohols  and  other  organic  liquids  affect  crack 
growth  in  Region  III.  The  influence  of  these  environments  is  indicated 
by  small  shifts  in  the  position  of  the  crack  propagation  curves  and  by 
small  decreases  in  the  slopes  of  these  curves  when  compared  to  the 


position  and  the  slope  of  the  curve  obtained  in  nitrogen  gas.  Tvo 
possible  mechanisms  were  considered  to  explain  our  experimental 
observations.  The  first  involves  viscous  drag  of  the  liquid  on  the 
freshly  formed  fracture  surfaces,  which  would  retard  crack  growth, 
causing  small  shifts  in  the  crack  velocity  curves.  If  viscous  drag 
were  an  important  mechanism,  then  the  crack  growth  curve  would  be 
expected  to  shift  to  higher  values  of  as  the  viscosity  of  the  liquid 
increased.  This  type  of  shift  in  the  position  of  the  crack  growth 
curve  is  in  fact  observed  for  the  alcohols.  The  possibility  of  a 
viscous  drag  contribution  to  the  crack  growth  behavior  was,  however, 
rejected  for  two  reasons:  (l)  the  shape  of  the  crack  growth  curves 
obtained  experimentally  differs  from  that  predicted  from  fluid  dynamic 
equations;  (2)  viscous  drag  cannot  explain  the  experimental  crack  growth 
behavior  at  low  velocities,  since  all  of  the  crack  growth  curves  would 
be  expected  to  Join  up  with  the  nitrogen  curve  as  the  viscous  drag 
effects  vanished.  To  support  these  conclusions,  detailed  discussion  of 
viscous  drag  is  given  in  Appendix  A. 

The  second  mechanism  considered  to  explain  our  experimental  findings 

involves  electrostrictive  interactions  between  the  glass  and  the  fluid  at 

the  crack  tip.  If  one  accepts  the  chemical  interaction  model  of 

subcritical  crack  growth  in  ceramic  materials,  then  the  slope  of  the 

crack  propagation  curves  in  regions  I  or  III  is  determined  by  the 

activation  volume,  AV^,  of  the  chemical  reaction  that  occurs  at  the 

crack  tip.13  When  fracture  is  caused  by  the  reaction  of  water  with 
♦ 

glass,  AV  represents  the  difference  in  molar  volume  between  the 


reactants  and  the  chemical  complex  in  the  activated  state.  If  fracture 
occurs  in  the  absence  of  water,  then  represents  the  change  in  volume 
that  occurs  as  the  strained  Si-0  bonds  move  into  an  activated  state. 
Referred  to  as  a  fission  reaction  in  the  chemical  literature,  this  is 
the  type  of  reaction  of  primary  interest  in  fracture  processes. 

When  a  fission  reaction  occurs  in  a  liquid  environment,  it  has 
been  shown  that  the  activation  volume  normally  consists  of  two  parts: 
one  part  relates  to  the  stretching  of  the  bond  that  is  being  broken; 
the  second  part  relates  to  volume  changes  in  the  fluid  surrounding  the 
reacting  molecule.15*23”25  The  second  part  is  usually  attributable  to 
electric  charge  formation  or  destruction  during  the  chemical  reaction, 
which  results  in  a  localized  electrostriction,  and  hence  a  change  in 
volume  in  the  fluid  surrounding  the  reacting  species.  This  change  in 
volume  must  be  included  as  part  of  the  activation  volume  because  the 
volume  change  is  concomitant  with  the  fiBBion  reaction.  For  normal 
chemical  reactions,  charge  formation  or  destruction  results  in  rather 
large  contributions  to  the  activation  volume.  In  fact,  the  electrostatic 
portion  of  the  activation  volume  is  usually  3  to  10  times  that  due  to 
changes  in  bond  length  alone.  Hence,  if  charges  are  either  formed  or 
destroyed  during  the  rupture  of  a  silicon-oxygen  bond  in  glass,  electro¬ 
static  effects  in  the  glass  during  crack  growth  will  have  to  be 


considered 


In  order  to  model  the  effect  of  environment  on  a  fission  reaction 
in  glass,  a  procedure  similar  to  that  used  by  Wiederhom,  Fuller  and 
Thomson13  is  employed.  The  procedure  models  the  crack  tip  as  an 
ellipse  and  uses  the  stress  at  the  major  axis  of  the  ellipse,  calculated 
from  the  Inglis  relation,17  as  the  driving  force  for  the  chemical 
reaction.  The  crack  tip  stress  then  replaces  pressure  as  the  thermo¬ 
dynamic  parameter  usually  used  to  discuss  chemical  reactions.  The  crack 
tip  stress  is  of  course  a  negative  pressure.  The  effect  of  charge 
formation  is  then  treated  by  elementary  solvation  theory:16'22  the 
charged  atoms  or  molecules  formed  during  the  reaction  being  approximated  by 
charged  conducting  spheres;  the  medium  surrounding  the  reacting  species 
being  treated  as  a  continuum  with  a  dielectric  constant,  e.  For  the 
fission  of  a  silicon-oxygen  bond  during  crack  growth  in  glass,  positive 
and  negative  charges  are  assumed  to  form  on  the  silicon  and  oxygen 
atoms  respectively  as  they  are  separated.  This  assumption  that  a  full 
unit  of  charge  develops  during  the  activated  state  is  based  on  the  work 
of  Akmed  et  al.,24  who  reported  charge  formation  on  silicon  and  oxygen 
atoms  during  the  crushing  of  silica  glass.  Suggestions  of  charge 
formation  during  fracture  were  also  made  in  earlier  glass  literature.27'28 

The  procedure  used  in  the  present  paper  to  estimate  the  electro¬ 
static  component  of  the  activation  volume  differs  somewhat  from  that 
used  by  Wiederhom,  Fuller  and  Thomson,13  in  that  the  reacting  species 
(the  oxygen  and  the  silicon)  are  assumed  to  lie  at  a  stressed  surface 
of  the  glass.  This  assumption  is  necessary  in  order  to  evaluate  the 
effect  of  environment  on  the  reacting  species.  Furthermore,  we  assume 


that  the  electrostatic  component  of  the  activation  volume  (or  activation 
free  energy)  is  entirely  the  result  of  charging  the  oxygen  and  silicon 
atoms  during  the  fission  process.  Finally,  we  assume  that  the  reaction 
is  sufficiently  slow  that  the  static  component  of  the  dielectric 
constants  can  be  used  to  estimate  the  electrostrictive  processes  at  the 
crack  tip,  and  that  both  the  glass  and  the  environment  at  the  crack  tip 
can  be  treated  as  a  continuum.  These  last  two  assumptions  are  usually 
made  to  treat  the  kinetics  of  chemical  reactions  in  liquids.15*23-21* 

A  schematic  diagram  of  the  model  used  in  the  present  paper  to 
account  for  the  electrostatic  component  of  the  activation  volume  is 
given  in  Figure  11.  Charges  are  assumed  to  form  at  a  flat  Interface 
during  the  reaction.  Each  Si  and  0  atom  in  the  activated  state  is 
assumed  to  react  independently  with  the  glass  and  the  surrounding 
environment  during  the  fission  reaction.  This  assumption  is  often 
made  in  evaluating  the  electrostatic  component  of  the  activation 
volume  and  tends  to  neglect  dipole  or  quadrupole  effects  on  reaction 
rate.  The  configuration  shown  in  Figure  11  was  chosen  because 
electrostatic  effects  of  the  dielectric  fluid  on  the  rupture  process 
are  easily  estimated  from  this  model.  The  model  does  not  take  into 
account  effects  that  result  from  curvature  of  the  crack  near  the 
crack  tip  or  stearic  hinderance  due  to  limits  in  the  dimensions  of 
the  crack  near  its  tip.  It  is  assumed  that  the  change  in  free  energy, 

AG  ,  when  a  charge  forms  at  the  interface  can  be  estimated  by  the  energy 
required  to  charge  conducting  spheres  of  radius,  r^  and  r^,  equal  to  the 
radius  of  the  silicon  and  oxygen  atoms  (0.04  and  0.14  nm)  respectively. 
Each  atom  is  assumed  to  develop  a  full  unit  of  charge  as  rupture  occurs, 


60 


a  positive  charge  for  the  silicon  atom  and  a  negative  charge  for  the 

oxygen  atom.  The  free  energy  of  formation  of  each  charged  atom,  AGeg, 

« 

is  given  by: 

AG  *  (N  e2/r)/(e.+e~)  (13) 

es  a  l  c 

per  mole  of  charge  formed  for  each  ion.  is  Avogadro's  constant,  and 
and  are  the  dielectric  constants  of  the  glass  and  the  liquid 
respectively. 

The  electrostatic  component  of  the  activation  volume  for  the  reac¬ 
tion  is  determined  from  equation  13»  using  the  relation,  AVgg  ■  3AGgg/«)P, 
where  P  is  the  pressure  that  constrains  the  chemical  reaction.  When  the 
reaction  occurs  at  the  surface  of  the  glass,  this  pressure  is  equal  to' 
the  negative  of  the  crack  tip  stress.  Hence  for  a  chemical  reaction 
constrained  by  a  pressure,  P,  the  activation  volume  A V  is  given  by: 

AVes  =  (Na  ®2/r)  t-(ei+e2)-:L  31n  r/3P  "  (e1+e2)“23e1/3P]  (1*0 

for  each  ion.  The  derivation  of  equation  (lU)  is  based  on  the  assumption 
that  the  pressure  is  applied  only  to  the  glass,  hence  the  partial 
derivative  is  taken  only  with  regard  to  e^.  Equation  lU  can  be  written 
in  terms  of  the  bulk  modulus  B: 

Aves  *  (Na  e2/r)[(3B)"1(e1+e2)"1  -  («1+*2>"2**1/* P]  {l5) 

_ 

This  equation  is  derived  in  Appendix  B. 


From  equation  15.  we  note  that  the  electrostatic  contribution  to  the 
activation  volume  depends  on  four  parameters:  the  size  of  the 
electrostatic  charge  formed  during  the  reaction;  the  value  of  the  bulk 
modulus;  the  value  of  the  dielectric  constants  (e^  and  e2)  and  the 
partial  derivative  of  the  dielectric  constant  with  regard  to  the  applied 
pressure.  As  already  mentioned,  equation  15  neglects  any  dipole,  or 
quadrupole  contributions  to  AV  .  The  partial  derivative  of  e1  with 
regard  to  P  is  usually  negative  for  solids;  however,  some  positive  values 
have  been  measured.29 When  3e^/3P  is  negative,  then  the  activation 
volume  is  positive.  However  when  3e^/3P  is  positive,  the  activation 
volume,  AV  ,  cm  be  positive  or  negative  depending  on  the  size  of  the 
two  terms  in  the  square  brackets.  Finally,  the  dependence  of  AV  on 
environment  depends  primarily  on  the  dielectric  constant,  e2,  of  the 
environment. 

Crack  growth  data  can  be  compared  with  the  electrostrictive  model 
for  crack  growth  by  plotting  AV*,  the  total  activation  volume  (equation 
9),  as  a  function  of  AV  ,  the  electrostatic  portion  of  the  activation 
volume.  It  is  worth  noting  that  AV*  determined  from  equation  9  contains 
both  the  stretching  and  the  electrostatic  components  of  the  activation 
volume.  For  crack  growth  studies  conducted  on  soda-lime-silica  glass  in 
different  environments,  AV  is  expected  to  depend  primarily  on  e0,  the 
dielectric  constant  of  the  fluid  at  the  crack  tip.  For  fluids  that  have 
a  high  dielectric  constant  (i.e.  water,  formamide)  AV  is  expected  to 
be  small,  and  as  a  consequence,  the  value  obtained  for  AV* will  be  due 
primarily  to  stretching  of  the  Si-0  bonds  by  the  applied  forces.  For 


environments  that  have  low  dielectric  constants  (i.e.  alkanes,  vacuum) 

AV  will  he  large  and  AV  ^  will  contain  components  due  both  to  electro- 

static  interactions  and  to  stretching.  Figure  12  presents  a  plot  of  AV 

* 

as  a  function  of  AVga.  Despite  the  experimental  scatter,  a  correlation 

between  AV  *  and  AV  is  apparent  from  the  figure.  The  figure 
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demonstrates  a  definite  dependence  of  crack  growth  data  on  the 
dielectric  constant  of  the  test  media,  and  as  such  lends  support  to 
the  idea  that  electrostatic  interactions  between  the  environment  and 
the  highly  strained  bonds  at  the  crack  tip  influence  crack  growth.  The 
scatter  of  the  data  in  Figure  12,  illustrates  the  difficulty  of 
obtaining  accurate  values  of  the  slope  of  the  v-Kj  data  in  region  III, 
especially  when  the  slopes  are  high,  as  they  are  in  alkanes,  nitrogen 
gas,  and  vacuum.  Further  confirmation  of  the  theory  will  require  the 


development  of  more  accurate  methods  of  determining  the  Blope  of  crack 
growth  curves.  These  techniques  should  then  be  applied  to  crack  growth 


The  authors  could  find  no  values  of  3c^/3p  in  the  literature  for 
soda-lime-silica  glass.  The  value  of  Se^/aP  (-9*72  x  10”11  Pa-1) 
used  for  Figure  12  was  obtained  by  Colwell90  on  an  alkali  containing 
glass  made  at  the  National  Bureau  of  Standards.  The  fluid  dielectric 
constants  e g  and  the  slopes,  b,  of  the  crack  growth  curves  used  for 
Figure  12  are  given  in  Table  2. 


studies  in  low  dielectric  constant  liquids,  <8,  where  electrostatic 
effects  are  expected  to  be  greatest.  Furthermore,  a  variety  of  glasses 
and  single  phase  ceramics  should  be  studied  to  evaluate  the  effect  of 
3e1/3P  on  crack  growth  data. 

5.  Summary 

This  paper  presents  the  results  of  a  series  of  studies  to  determine 
the  effect  of  water  and  a  variety  of  organic  liquids  on  crack  growth  in 
soda-lime-silica  glass.  Throughout  the  paper  the  crack  growth  process 
is  viewed  as  a  chemical  reaction  and  the  data  are  interpreted  in  terms 
of  chemical  reaction  rate  theory.  In  the  presence  of  a  chemically 
active  environment,  rupture  of  the  Si-0  bondB  of  the  glass  network  1b 
viewed  as  a  fission  reaction  in  which  Si-0  bonds,  strained  by  the  crack 
tip  stresses,  are  cleaved  by  the  chemical  environment.  In  the  absence 
of  a  chemical  reaction  between  the  glass  and  the  environment  the  fission 
reaction  occurs  as  a  result  of  a  stress  induced  cleavage  of  the  Si-0 
bonds  of  the  glass  network  structure.  The  three  regions  of  crack  growth, 
observed  in  crank  growth  studies  are  interpreted  in  terms  of  reaction 
rate  theory. 

In  region  I,  crack  growth  in  organic  liquids  is  the  result  of  a 
chemical  reaction  between  the  active  component  in  the  environment 
(usually  H2O)  and  the  glass;  the  rate  of  growth  is  controlled  by  the 
chemical  potential  of  the  active  component.  This  finding  confirms  those 
of  earlier  studies  and  is  consistent  with  the  chemical  reaction  rate 


theory  of  crack  growth.  In  region  II  the  rate  of  crack  growth  is 
determined  by  diffusion  of  the  reactant  to  the  crack  tip.  We  show 
that  crack  growth  in  region  II  is  not  controlled  by  the  chemical  potential 
of  the  reactant,  but  by  its  concentration.  Using  the  Stokes-Einstein 
relation  for  the  diffusivity,  a  linear  relation  is  obtained  between 
the  crack  velocity  and  the  ratio  of  the  concentration  of  water  to 
the  viscosity  of  the  organic  fluid.  Finally  in  our  discussion  of 
crack  growth  in  region  III,  we  Introduce  the  idea  that  electrostrictive 
effects  at  crack  tips  can  influence  the  crack  growth  behavior  of  solids. 
The  concept  of  an  electrostatic  component  to  the  activation  volume  is 
commonly  used  by  physical  chemists  to  explain  the  effects  of  solvents 
on  rates  of  chemical  reactions.  In  this  paper,  a  model  was  developed 
to  quantify  the  effect  of  dielectric  media  on  crack  growth  rate  in  glass. 
The  model  follows  classical  chemical  lines  and  has  the  virtue  of  providing 
a  relation  between  the  slope  of  v-K^  data  and  the  dielectric  constants  of 
the  glass  and  the  test  medium.  A  plot  of  the  total  activation  volume 
obtained  from  the  slope  of  crack  growth  data  versus  the  electrostatic 
component  of  the  activation  volume  calculated  from  theory,  lends  support 
to  the  model  and  suggests  that  electrostrictive  effects  influence  the 
growth  of  cracks  in  glass.  Despite  this  agreement,  the  approximate 
nature  of  the  model  should  be  recognized.  Further  studies  may  show 
that  some  variables  neglected  in  the  model  are  in  fact  important  to  our 
understanding  of  region  III  crack  growth.  Clearly,  additional  data  is 
needed  to  fully  evaluate  the  theory  presented  in  this  paper. 
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APPENDIX  A 


The  Effect  of  Viscous  Plow  on  Crack  Motion 


In  an  earlier  section  of  this  paper,  the  idea  that  the  position  of 
the  v-K-j.  curves  in  alcohol  could  be  explained  by  viscosity  effects  was 
rejected.  In  this  section  we  develop  equations  that  relate  the  crack 
motion  to  the  viscosity  of  the  fluid  at  the  crack  tip  and  show  why  the 
idea  was  rejected. 


Assume  that  the  crack  can  be  represented  by  a  parabola  that 
intersects  an  external  surface,  Figure  1A.  The  parabola  is  selected 
because  the  crack  has  the  shape  of  a  parabola  near  its  crack  tip  in 
the  linear  elastic  approximation.  The  crack  length  is  equal  to  c,  the 
y-coordinate  is  perpendicular  to  the  plane  of  the  crack  and  the  x 
coordinate  lies  in  the  plane  of  the  crack.  As  the  crack  moves  along, 
it  sucks  fluid  into  the  crack  opening,  which  is  exposed  to  cm  infinite 
reservoir  of  fluid  at  an  ambient  pressure,  P  .  In  order  for  the  fluid 

fit 

to  flow  to  the  crack  tip,  the  fluid  pressure  decreases  from  the  crack 
opening  to  the  crack  tip.  Therefore  the  fluid  within  the  crack  is  at  a 
negative  pressure  relative  to  the  ambient  pressure,  P  .  The  negative 
pressure  loads  the  surface  of  the  crack  in  such  a  way  as  to  tend  close 
the  crack.  Thus,  the  closing  force  on  the  crack  surface  produces  a 
negative  stress  intensity  factor,  which  opposes  the  positive,  applied 
stress  intensity  factor,  K  ,  that  is  driving  the  crack  to  propagate. 


To  calculate  the  decrease  in  pressure  from  the  crack  opening  to  the 
crack  tip,  and  to  estimate  the  negative  stress  intensity  factor,  we 
assume  that  the  coordinate  system  is  fixed  to  the  crack  tip.  Then  the 
crack  surface  appears  to  recede  from  the  crack  tip  at  a  velocity  vq, 
which  is  also  the  velocity  of  the  crack.  The  half  width  of  the  crack, 
y  ,  is  determined  by  the  standard  (plane  strain)  fracture  mechanics 
relation  for  the  displacement  in  the  vicinity  of  a  crack  that  is  subject 
to  a  stress  intensity  factor  (see,  for  example.  Ref.  34) 

yQ2  =  [2KI(l-v)/G]2(x/2ff)  (1A) 

where  G  is  the  shear  modulus  of  the  solid  and  v  is  Poisson's  ratio. 

The  stress  intensity  factor, . K^,  therefore  describes  the  stresses  and 

displacements  at  the  crack  tip  and  has  to  be  equal  to  the  sum  of  K&  and 

K  :  Kt  *  K  +  K  .  This  approach  to  the  crack  tip  stresses  is  similar 
p  I  a  p 

to  the  approach  used  by  Bareriblatt35  in  his  discussion  of  crack  tip 
forces. 

At  any  distance,  x,  from  the  crack  tip,  the  pressure  drop  across  a  • 
section,  dx,  of  the  crack  is  given  by  the  Poiseille  flow  equation  for  a 
parallel  sided  channel : 

vQ  =  (l/3n)  (dP/dx)  yQ2  (2A) 

where  n  is  the  viscosity  of  the  fluid  in  the  crack.  As  noted  by 
Batchelor,36  this  approximation  is  valid  provided  the  slope  of  the 


walls  of  the  channel  change  slowly  with  distance  along  the  channel. 
Substituting  equation  LA.  into  equation  2A  and  integrating  the  pressure 
drop  from  a  point  within  the  crack  to  the  crack  mouth,  the  following 
equation  is  obtained  for  the  pressure,  P,  at  any  point  within  the 
channel : 

P  =  Pa  +  6irvQn  [G/2Kj(l-v)]2  In  (x/c)  (3A) 

Note  that  as  the  position  x  approaches  the  crack  tip  (i.e.  x  ■*  o),  large 
negative  pressures  occur  within  the  fluid. 

The  negative  stress  intensity  factor ,  K^,  that  results  from  the 
reduction  of  pressure  within  the  crack  can  be  calculated  from  a  standard 
equation  for  a  two  dimensional  crack  vithin  a  solid  that  Is  subject  to 
an  applied  stress  o(x):  3I* 

K  -  2  (c/ir)1/2  f  o(x’)  dx'/[c2-(x')2]1/2  (UA) 

p  *o 

where  x'  =  c-x  measures  the  distance  from  the  crack  mouth. 

The  stress  o(x')  at  any  point  on  the  crack  surface  is  determined 
from  the  difference  between  the  pressure  within  the  crack  and  the 
external  pressure:  o(x')  *  P  -  P&.  Note  from  3A  that  o(x' )  is  a 
negative  quantity,  as  it  should  be  for  a  closing  force  on  the  crack  tip. 
Substituting  the  equation  for  o(x')  into  equation  hk  and  integrating,36 
the  following  equation  is  obtained  for  Kp: 


Kp  =  -1.5  (irln2  +  4C)  VQn  [G/IC^l-v]  (ire )±,c 


=  -8.76  v  n  [G/K  (1-v)]2  (irc)l/2 
o  1 

where  C  =  0.915966  is  Catalan's  constant. 


(5A) 


For  a  constant  crack  length,  c,  Kp  depends  only  on  the  crack 
velocity,  vq,  and  the  crack  tip  stress  intensity  factor,  K^;  since  Kj 
is  a  slowly  varying  function  of  vq  for  ceramic  materials,  Kp  depends 
primarily  oa  V 

Normally  when  a  v-K^.  curve  is  determined  experimentally,  it  is  K& 

that  is  measured  and  presented  as  a  plotting  parameter.  As  long  as  Kp 

is  small,  the  v-K^  curve  and  the  v-K&  curve  will  be  indistinguishable. 

However  when  Kp  becomes  large  relative  to  Kft,  the  experimental  curve 

will  diverge  from  the  "true"  v-K^  curve.  This  occurs  when  part  of  the 

applied  load  is  transferred  from  the  specimen  to  the  fluid  within  the 

crack.  The  range  of  velocities  over  which  this  occurs  can  be  calculated 

by  evaluating  equation  5A  and  comparing  it  with  crack  tip  stress 

intensity  factor.  As  shown  in  Table  1A,  where  P-P  and  K  are  calculated 

a  p 

as  a  function  of  velocity  for  water,  the  transition  range  over  which 
these  viscous  effects  are  important  occur  over  a  relatively  narrow  range 
of  crack  velocities.  The  data  in  Table  1A  cure  plotted  in  Figure  2A  to 
give  a  predicted  v-K  curve.  We  note  that  the  curve  bends  over  rather 

a 
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sharply ,  over  a  velocity  range  of  about  two  orders  of  magnitude. 
Behavior  of  this  type  has  been  observed  for  water  by  Schdnert  et  al.2 
and  by  Michalske  and  Frechette  for  water  and  for  methyl  alcohol.12 
Since  the  shape  of  the  curves  observed  in  the  present  paper  are  not 
similar  to  those  obtained  for  the  alcohols  in  region  III,  we  conclude 
that  viscous  effects  of  the  crack  tip  fluid  played  little  if  any  role 
in  the  present  experiments,  thus  Justifying  our  statements  in  section 


APPENDIX  B 


Si 


Electrostatic  Energy  of  a  Charged  Sphere 
on  a  Dielectric  Interface 


The  model  used  in  section  4.3  to  describe  the  region  III  crack  growth 


behavior  relates  the  change  in  free  energy  during  the  crack-tip  fission 


reaction  to  the  electrostatic  energy  required  to  charge  a  conducting  sphere 


embedded  on  a  dielectric  interface  (see  Figure  12).  In  this  appendix  we 


derive  an  expression  for  this  electrostatic  energy. 


In  a  linear  dielectric  medium,  the  electrostatic  energy  of  a  systc 


of  charges  is  related  to  the  electric  field  E  that  they  produce  by  the 


relation:3® 


t  /•  m‘ 


where  e  is  the  dielectric  constant  and  the  volume  integral  is  over  ell 


space.  Accordingly,  the  problem  reduces  to  finding  the  electric  field  E, 


or  its  scalar  potential  t,  since  then  E  ■  -  V*,  and  integrating  the  results 


according  to  Eq.  (IB). 


The  electrostatic  potential  *  is  determined  as  a  solution  of  Poisson's 


equation: 


V2*  »  -  4lfv 


s 


where  x  is  the  charge  density.  The  boundary  conditions  for  the  present 
problem  require:  (l)  that  the  surface  of  the  conducting  sphere  is  at  a 
constant  potential;  (2)  that  the  displacement  field,  eE,  normal  to  the 
sphere  is  equal  to  4n  times  the  surface  charge  density;  (3)  that  the 
surface  charge  density  integrated  over  the  surface  of  the  sphere  is  Q; 
and  (h)  that  the  tangential  components  of  E  is  continuous  across  all 
boundaries.  Rather  than  solve  Eq.  (2B)  directly,  we  note  that  a  solution 
of  Poisson's  equation  is  unique  when  the  potential  is  specified  on  all 
closed  conducting  surfaces  (Dirchlet  boundary  conditions),  as  in  the 
present  problem.  Thus,  if  we  obtain  a  solution  by  any  means,  we  have 
the  unique  solution. 

A  solution  to  the  present  problem  is  obtained  from  the  electrostatic 
potential  of  a  point  charge  Q  embedded  on  a  planar  interface  betveen  to 
two  semi-infinite  dielectrics.  This  solution:38 


has  spherical  equipotential  surfaces  and  hence  gives  the  proper  form  of 
the  solution  for  distance  greater  than  the  radius  r  of  the  conducting 
sphere  (i.e.,  R>r).  Inside  the  sphere  the  potential  is  constant : 


sa 


t 


$ 

c 


The  electric  field  vector  is  given  by  the  negative  of  the  gradient  of  the 


potential ,  and  only  has  a  radial  component  since  the  potential  is 


spherically  symmetric: 


t  for  R<r 


[eg  +  e^)Rz  ,  for  r<R 


It  is  easy  to  verify  that  this  solution  satisfies  the  remaining  boundary 


conditions. 


The  electrostatic  energy  required  to  charge  the  conducting  sphere  is 


obtained  by  integrating  Eq.  (5B)  in  Eq.  (IB): 


68  "dielectric  2  *2 


<  43R  *  *£/electrlc  *1<  ^ 


,e2  +  el,r 


In  section  k.3  the  charge  Q  is  taken  to  be  a  full  unit  of  electronic 


charge,  Q  *  e. 


Table  1:  Viscosity,  Water  Concentration  and  Relative  Humidity  of  Alcohols 
used  in  Region  II  Crack  Growth  Study. 


Alcohol 

Viscosity 

(Pa*s) 

Water  Concentration 
(mol/liter) 

Percent  Relation 
Humidity 

Butyl 

3.13  x  10”3 

9.2 

100 

Amyl 

3A6  x  10"3 

k.OQ 

100 

Hexyl 

3.95  x  10"3 

2.0k 

100 

Heptyl 

5.16  x  10”3 

2.52 

100. 

Octyl 

6.31  x  10“3 

1.13 

50 

Decyl 

10.99  x  10"3 

1.68 

100 

Ethyl  (1%) 

-3* 

1.1  x  10  3 

0.555 

7 

Butyl  (0.2%) 

-3* 

2.95  x  10  J 

o.m 

5 

Butyl  (151) 

2.95  x  10"3* 

0.555 

21 

Butyl  (55*) 

2.95  x  10“3# 

2.78 

70 

# 

The  viscosities  of  the  partially  saturated  ethyl  andQbutyl  alcohols 
were  assumed  to  be  equal  to  the  pure  alcohols  at  20  °C. 


Table  2:  A  comparison  of  the  slopes,  b,  of  crack  growth  data  with  the 

dielectric  constant,  e.,  of  the  test  medium.  All  crack  growth 
measurements  were  made  on  soda- lime— silica  glass;  slopes  were 
determined  from  an  empirical  fit  to  the  equation:  v  »  vQexp 
bKj/RT. 


Test 

Dielectric 

Slope,  b 

Test 

Medium 

Constant, 

MKS  units 

Method 

Formamide 

109 

0.168 

D.C.B. 

Water 

78 

0.111 

D.C.B.,  ref.  31 

Acetonitrile 

39 

0.232 

Double  Torsion 

Methanol 

33 

0.182 

Double  Torsion 

0.169 

Double  Torsion 

Hexanol 

13.3 

0.223 

D.C.B. ,  ref.  33 

Heptanol 

12.2 

0.283 

D.C.B. ,  ref.  33 

Octanol 

10.3 

0.228 

D.C.B.,  ref.  33 

Decanol 

8.0 

0.240 

D.C.B. ,  ref.  33 

0.252 

Double  Torsion 

0.391 

Double  Torsion 

Trichloroethylene 

3.4 

0.463 

D.C.B. 

Alkanes 

2 

0.488 

D.C.B.,  ref.  32 

Heptane 

2 

0.361 

Double  Torsion 

0.772 

Double  Torsion 

0.724 

D.C.B.,  ref.  32 

Nitrogen 

1 

0.636 

Double  Torsion 

0.542 

Double  Torsion 

0.823 

D.C.B.,  ref.  1 

Vacuum 

1 

0.88 

D.C.B.,  ref.  6 

Table  1A.  Effect  of  viscous  forces  on  crack  notion  calculated  from 

equations  3A  and  5A.  The  pressure  Is  calculated  for  a  1  ran 
crack  at  a  distance  of  1  in  from  the  crack  tip.  (G  »  28  GPa, 


V  ■ 

0.2k  and  n  »  1  nPa*s) 

Crack  Velocity 

(m/s) 

(mPa»n’) 

P‘Pa 

(nPa) 

K 

p  «s 
(mPaun*) 

1.00  x  10_1 

0.570 

-13.6 

-0.205 

5.72  x  10~2 

0.560 

-  8.1 

-0.122 

2.63  x  10"2 

0.51*3 

-  3.9 

-0.059 

1.20  x  10”2 

0.530 

-  1.9 

-0.028 

5.25  x  10“3 

0.505 

-  0.91 

-O.OlU 

2.U0  x  10~3 

0.U90 

-  O.kk 

-0.007 

1.00  x  10“3 

0.U70 

-0.20 

-0.003 
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Figure  2.  Crack  propagation  in  ethyl  alcohol  containing  water.  The 
concentration  of  water  is  expressed  In  volume  percent. 
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PERCENT  RELATIVE  HUMIDITY 


Figure  5.  Crack  growth  rate  in  region  II  as  a  function  of  relative 
humidity.  The  data  was  taken  from  figures  2  and  3  using 
the  plateau  values  of  the  velocity  and  relative  humidities 
calculated  from  the  concentration  using  thermodynamic  data 
reported  in  references  7  and  8. 


Figure  6.  Creek  growth  date  on  ■oda-llne-illica  glees  collected  tor 
the  double  torsion  technique  (relexetion  method)  in 
alcohols.  With  the  exception  of  the  octyl  alcohol,  all 
alcohols  were  fully  saturated.  The  octyl  alcohol  was  *50 
percent  saturated. 
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Figure  8.  Region  III  crack  growth  data  collected  on  silica  glass 

(C79U0)  by  the  double  cantilever  bean  technique  (constant 
monent  method). 
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Figure  12.  Plot  of  the  total  activation  volume  (eq.  9)  aa  a  function 
of  the  electrostatic  component  of  the  activation  volume 
(eq.  15).  The  value  of  8s,/3P  (-9*72  x  10"11  Pa"1)  use  taken 
from  data  by  Colvell30  collected  on  an  alkali  containing  glass 
made  at  NBS.  Double  cantilever  beam  data  are  indicated  by 
circles,  double  torsion  data  by  triangles.  Data  taken  from 
the  literature  are  identified  by  reference  numbers  on  the 
figure.  ACH  is  an  abreviation  for  acetonitrile;  PHD  is  an 
abbreviation  for  formamlde;  TCE  is  an  abbreviation  for 
trlchlorethylene.  The  other  chemical  symbols  are  standard 
chemical  designations. 
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Figure  2A:  Calculated  effect  of  fluid  viscosity  on  v-K_  curve:  soda- 
liae-silica  glass  in  water.  The  K_  portion1 of  the  curve  Is 
linear  extrapolation  of  data  collected  at  velocities  less 
than  10"  3  a/s. 
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ABSTRACT 

The  effect  of  subcritical  crack  growth  on  proof  testing  Is  examined. 
Crack  velocity  curves  obtained  by  fracture  mechanics  techniques  are  used  to 
predict  theoretical  strength  distributions  for  specimens  that  survive  proof 
testing.  These  theoretical  distributions  are  then  compared  with  experimental 
distributions  obtained  on  soda-lime  silica  glass  slides.  The  comparison 
reveals  a  strong  sensitivity  of  the  proof  test  results  to  the  exact  position 
and  shape  of  the  crack  growth  curve.  Minor  changes  in  the  crack  growth  curve 
result  in  major  shifts  in  position  and  shape  of  the  strength  distribution 
curves  after  proof  testing.  The  importance  of  crack  geometry  and  specimen 
configuration  to  crack  growth  behavior,  and  hence  to  the  strength 
distribution,  is  emphasized. 


1.  INTRODUCTION 


In  order  to  Improve  the  reliability  of  ceramic  materials  In  structural 

applications,  proof  testing  Is  often  used  as  a  method  of  controlling  the 

strength  distribution.  Weak  components  are  broken  by  the  proof  stress  leaving 

a  population  of  components  with  a  minimum  value  of  the  strength,  which  In  the 

absence  of  subcrltlcal  crack  growth,  Is  equal  to  the  stress  used  during  the 

proof  test.  Thus,  proof  testing  truncates  the  strength  distribution,  and 

guarantees  that  all  components  will  be  at  least  as  strong  as  the  minimum 

value.  Because  of  Its  simplicity,  proof  testing  has  been  used  successfully  as 

a  basis  of  design  in  many  structural  applications.*”^ 

Potential  problems  in  the  application  of  proof  testing  occur  when 

subcritical  crack  growth  accompanies  proof  testing.  An  example  of  this  Is 

seen  In  a  recent  study  on  soda-lime  silica  glass. tested  In  air  or  water.1*  It 

was  shown  that,  for  certain  test  conditions,  the  proof  test  did  not  truncate 

the  strength  distribution.  Furthermore,  test  results  did  not  fully  agree  with 

12 

the  theory  that  was  used  to  explain  them.  The  study  suggested  the  need  for 

additional  research  on  the  effect  of  multiregion  crack  growth  on  proof 

testing.  It  is  because  of  this  need  that  the  present  study  was  conducted. 

13 

The  study  is  similar  to  one  recently  completed  by  Ritter  et  al.  ,  and  the 
results  are  complementary  to  those  obtained  by  those  authors. 

2.  EXPERIMENTAL  PROCEDURE 

Heptane  was  chosen  as  a  test  medium  for  this  study  because  of  its  low 
affinity  for  water.  At  normal  temperatures  and  pressures  heptane  contains 
only  ~  50  ppm  water,  which  corresponds  to  a  partial  pressure  of  water  equal  to 
30%  relative  humidity.  The  net  effect  of  the  low  concentration  of  water  is  to 
suppress  region  II  crack  growth  behavior  to  a  velocity  of  ~  10  ®  m/s,  but  to 
leave  the  curve  for  region  I  crack  growth  virtually  unchanged.  The  range  of 
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the  crack  growth  curve  In  region  II,  0.5  MPa-m  to  0.7  MPa-nn  ,  Is  greater 
than  that  obtained  In  most  other  environments,  and  therefore,  is  expected  to 
enhance  the  effect  of  multi  region  crack  growth  on  proof  testing. 

Proof-test  specimens  were  unannealed  microscope  slides  with  dimensions  of 
75  x  25  x  1  mm.  Test  conditions  and  the  number  of  specimens  used  in  each  test 
are  given  in  Table  1.  All  specimens  were  loaded  to  the  proof  load  in 
four-point  bending  at  a  rate  of  33  MPa/s.  As  the  equipment  was  automated  to 
reduce  the  load  as  soon  as  the  proof  test  load  was  reached,  the  hold  time  at 
the  proof  test  load  was  less  than  0.1  s  and  is  assumed  to  be  zero  in  the 
present  study.  Based  on  crack  growth  data  reported  in  reference  14,  and 
preliminary  strength  measurements  In  heptane,  the  unloading  rates  were 
selected  to  yield  the  maximum  range  of  predicted  effects  of  proof  testing  on 
the  strength  distribution  curve  after  proof  testing.  Following  the  proof 
test,  specimens  were  loaded  to  failure  at  a  loading  rate  of  33  MPa/s  without 
changing  their  position  in  the  test  fixture.  Thus,  the  stress  distribution 
within  the  test  specimen  was  the  same  when  the  strength  was  measured  as  it  was 
during  the  proof  test. 

Crack  velocity  data  were  obtained  using  two  fracture  mechanics 

configurations.  Applied-moment,  double  cantilever  beam  specimens 

(75  x  25  x  1  mm)  were  used  to  collect  crack  velocity  data  over  the  velocity 

range  10  ^  to  10~^  m/s;  crack  velocities  were  measured  by  direct  observation 

15 

of  the  crack  tip  as,,a  function  of  time  .  Single-edge,  notched  tensile 

specimens  (130  x  25  x  2  mm)16  were  used  to  collect  data  over  the  velocity 

-9  -9  -7 

range  10  to  100  m/s.  In  the  velocity  range  10  to  10  m/s,  velocities 

were  measured  by  monitoring  the  crack  tip  visually.  Higher  velocity 

measurements  were  made  by  the  use  of  an  acoustic  pulse  technique  described  in 

detail  elsewhere.^  A  sonic  pulse  was  used  to  mark  the  fracture  surface  so 


that  the  position  of  the  crack  front  was  established  as  a  function  of  time. 
Velocities  were  then  determined  after  completion  of  the  experiment  by 
microscopic  examination  of  the  crack  surface. 

To  compare  the  crack  velocity  data  with  the  strength  data  a  computer  code 

was  designed  to  calculate  the  breaking  stress  of  the  specimens  after  the  proof 

test.  The  computer  code  is  a  modified  .version  of  the  one  discussed  in 
reference  12.  Input  data  for  the  code  are:  (1)  the  Weibull  parameters  SQ  and 
m  obtained  in  the  test  environment  (heptane);  (2)  a  v-Kj  plot  expressed  as  a 
piecewise  power-law  fit  to  the  v-Kj  data;  (3)  the  fracture  mechanics 
parameters  KIc*  KIscc  and  Y  (Y  Is  defined  by  the  equation  Kj  =  Yo/a  where  a  Is 

the  characteristic  crack  size,  o  the  applied  stress  and  Kj  the  stress 

A 

intensity  factor)  ;  and  (4)  the  proof  test  parameters  (stressing  rate  during 
loading  and  unloading,  proof  stress,  and  the  hold  time  which  is  assumed  to  be 
zero  in  the  present  experiment). 

The  computer  code  first  calculates  the  initial  distribution  of  strength 

from  the  initial  distribution  obtained  in  heptane  and  the  input  data 
** 

enumerated  above.  This  procedure  provides  a  common  basis  for  comparing 

*  1/2 
The  value  of  Kjc  used  in  the  present  study  was  0.75  MPa-m  .  The  reason 

for  selecting  this  value  is  discussed  in  the  Appendix. 

** 

The  initial  values  of  m  and  S,.  determined  in  this  study  are  listed  in 

o 

Table  la  of  the  Appendix.  This  calculation  is  needed  because  subcritical 
crack  growth  occurs  during  the  strength  tests.  The  initial  distribution, 
calculated  from  the  computer  code,  is  an  Idealized  distribution  that  would 
be  determined  if  strength  measurements  could  be  made  without  the  occurrence 
of  subcritical  crack  growth. 
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crack  growth  data  and  proof  test  data,  and  was  used  by  Ritter  et  al.  in 
their  study.  After  the  initial  distribution  has  been  determined,  the  computer 
code  determines  the  effect  of  the  proof  test  cycle  on  this  distribution,  by 
first  calculating  the  degree  of  strength  degradation  (i.e.  the  amount  of  crack 
growth)  during  the  proof  test  cycle,  and  then  by  calculating  the  breaking 
stress  in  heptane  of  the  specimens  that  survived  the  proof  test  cycle.  The 
results  of  this  analysis  can  then  be  compared  directly  with  experimental 
strength  distribution  data  obtained  in  heptane.  The  main  output  of  the 
computer  code  is  a  graphical  representation  in  which  Wei bull  probability 
coordinates  are  used  to  display  the  strength  distribution  after  proof  testing. 
On  command,  the  computer  code  was  capable  of  plotting  the  initial  strength 
distribution  in  heptane,  the  calculated  initial  distribution,  and  the 
distribution  after  proof  testing. 

The  value  of  Y  (1.5)  used  in  the  present  paper  is  representative  of  the 

maximum  stress  intensity  factor  at  a  semi-elliptical  surface  crack  with  a 

18 

minor  to  major  axis  ratio  of  ”0.6.  This  ratio  was  observed  in  the  present 
study  for  indentation  cracks  in  glass  surfaces  (see  fig.  7,  below).  We  assume 
that  Y  remains  constant  as  the  surface  crack  increases  in  size.  This 
assumption  permits  us  to  use  a  one-dimensional  model  to  determine  the  strength 
distribution  after  proof  testing.  For  real  cracks,  the  shape  of  the  crack, 
and  consequently,  the  value  of  Y  probably  changes  during  crack  growth. 

Although  this  aspecjt  of  crack  growth  is  Important,  the  more  simple  treatment 
adopted  here  should  permit  us  to  identify  the  principal  factors  that  influence 
strength  degradation  during  proof  testing. 


3.  RESULTS 


Crack  Growth  Data 

The  results  of  the  crack  velocity  measurements  are  given  in  figure  la. 

Approximately  100  data  points  were  collected  by  each  of  the  two  techniques. 

Data  for  the  single-edge,  notched  (SEN)  specimens  were  spread  evenly  over  the 

velocity  range  10  to  10  m/s;  the  double  cantilever  beam  (DCB)  data  were 

-10  -3  * 

spread  evenly  over  the  velocity  range  10  to  10  m/s  .  Thus,  the  data  from 
the  single  edge  notched  specimens  are  concentrated  in  a  higher  range  of  crack 
velocities  than  the  data  from  the  double  cantilever  beam  specimens.  Only 
5  data  points  were  obtained  on  single  edge  notched  specimens  at  velocities  of 
<  10~7  m/s. 

Because  the  two  sets  of  velocity  data  overlap  over  much  of  the  data 

range,  differences  between  the  data  are  shown  more  clearly  by  fitting  each  set 

** 

of  data  piecewise  with  straight  lines,  figure  lb  .  When  this  is  done,  the 

crack  growth  data  obtained  by  the  single-edge,  notch  technique  appear  to  lie 

at  slightly  higher  values  of  Kj  than  the  data  obtained  by  the  double 

cantilever  beam  technique.  A  small  lateral  translation  of  the  two  sets  of 

1/2 

data  by  approximately  0.05  MPa-m  brings  them  into  close  coincidence,  and 
for  practical  purposes,  makes  the  two  sets  of  data  indistinguishable  over  most 
of  the  data  range. 

-  j  . 

The  abbreviation  SEN  and  DCB  are  changed  to  SE(T)  and  DB(Mx)  in  the  E616-82 
Standard  Terminology  Relating  to  Fracture  Testing.  A  discussion  of  the 

terminology  will  be  found  in  the  1983  Annual  Book  of  ASTM  Standards. 

**  ... 

A  summary  of  the  piecewise  linear  fit  to  the  crack  growth  data  is  given  in 


Table  2a  of  the  Appendix. 


Ey: 


The  difference  between  the  two  sets  of  date  In  region  II  possibly  results 
from  the  difference  In  thickness  of  the  specimens  used  In  the  two  techniques. 
The  greater  thickness  of  the  SEN  specimens  Increases  the  diffusion  path  for 
water  from  the  bulk  to  the  crack  tip  environment  and  suppresses  the  region  II 
data  to  lower  crack  velocities.  This  type  of  data  has  been  reported  earlier 
by  Richter.^ 

-4 

At  velocities  that  exceed  ~  10  m/s,  the  slope  of  the  v-Kj  curve  is 

observed  to  decrease,  almost  as  if  there  were  a  second  region  II  behavior.  At 
1/2 

Kj  ~  0.93  MPa-m  the  crack  velocity  is  observed  to  increase  precipitously  by 

about  one  to  two  orders  of  magnitude.  This  sudden  increase  in  velocity  has 
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whereas  the 

10  *  m/s  is 
20 


been  attributed  to  cavitation  of  the  fluid  near  the  crack  tip,' 
decrease  In  slope  of  the  v-Kj  curve  at  velocities  that  exceed 
probably  due  to  viscous  drag  of  the  fluid  on  the*fracture  surfaces.' 

Strength  Data 

The  initial  strength  distribution  obtained  by  breaking  628  microscopic 
slides  in  heptane  Is  shown  in  figure  2.  Expressed  as  a  two  parameter  Weibull 
distribution,  the  parameters  for  the  regression  line  are  SQ  =  112  MPa  for  the 
intercept,  and  m  =  7  for  the  slope.  Despite  the  large  number  of  specimens 
tested,  the  curve  still  exhibits  experimental  scatter  about  the  regression 
line.  The  applicability  of  the  strength  distribution  to  the  specimens 
subjected  to  the  proof  test  can  be  checked  qualitatively  by  comparing  the 
number  of  specimens^roken  during  the  loading  portion  of  the  proof  test  cycle 
with  those  predicted  from  figure  2.  This  procedure  is  permissible  since  the 
loading  rate  was  the  same  for  all  4-point  bend  tests  in  this  study.  As  can  be 
seen  from  Table  2,  the  measured  and  predicted  number  of  failures  differ  by 
1  percent  for  the  medium  and  slow  unloading  rates.  The  difference  of 
7  percent  (34  versus  41  percent  failures)  for  the  rapid  unloading  rate  can  be 
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accounted  for  by  a  four  percent  error  in  either  the  proof  test  level,  or  SQ 
for  the  proof  tested  specimens.  In  view  of  this  small  difference  in  strength, 
we  feel  justified  in  using  the  Weibull  parameters  determined  from  figure  2  to 
describe  the  underlying  strength  distributions. 

Results  of  the  proof  test  study  (figure  3)  indicate  the  occurrence  of  a 

large  number  of  failures  below  the  proof  stress  for  all  three  unloading  rates. 

For  the  medium  unloading  rate  approximately  42  percent  of  the  specimens  failed 

below  the  proof  stress  when  they  were  reloaded  to  failure.  For  the  rapid  and 

slow  unloading  rates  16  percent  of  the  survivors  failed  below  the  proof  stress 

upon  reloading  to  failure.  This  decrease  in  strength  Is  consistent  with 

theoretical  expectations,  since  subcritlcal  crack  growth  is  expected  to  weaken 

all  specimens  that  survive  the  proof  test.  Specimens  that  are  weakened  most 

by  the  proof  test  are  those  with  strengths  that  lie  closest  to  the  proof  test 

stress  prior  to  testing.  The  weakest  specimen  after  testing  had  strengths  of 

74,  61  and  50  percent  of  the  proof-stress  for  the  rapid,  medium  and  slow 

unloading  rates  respectively.  This  finding  agrees  with  that  reported  by 
13 

Ritter  et  al.  who  noted  that  the  strength  degradation  of  specimens  just 
passing  the  proof  test  is  greater  as  the  unloading  rate  is  decreased. 

The  strength  distribution  curves  shown  in  figure  3  give  little  indication 
that  the  proof  test  truncated  the  distribution  curves.  Only  the  fastest 
unloading  rate  gave  evidence  of  truncation.  However  for  this  curve,  a  large 
number  of  specimens^Jbroke  at  strengths  that  were  less  than  the  proof-stress. 
Strength  distribution  curves  for  the  medium  and  slow  unloading  rates  differ 
only  slightly  from  the  initial  strength  distribution.  Because  of  the  large 
number  of  specimens  tested,  however,  we  believe  that  the  small  changes  in 
shape  and  position  of  strength  distribution  curves  shown  in  figure  3  are 
statistically  significant.  As  will  be  shown,  these  changes  in  shape  and 


position  can  ba  explained  by  the  occurrence  of  subcrltlcal  crack  growth  during 
the  proof-test. 

4.  COMPARISON  OF  CRACK  VELOCITY  OATA  WITH  PROOF  TEST  DATA 

Strength  Distribution  Curves 

A  comparison  of  the  strength  distribution  curves  with  the  curves 
calculated  from  the  v-Kj  data  are  given  In  figure  4.  The  straight  lines  in 
this  figure  represent  the  initial  distribution  of  specimen  strengths  obtained 
in  heptane.  The  data  points  were  selected  from  figure  3  to  Indicate  the 
location  of  the  data;  not  all  of  the  data  points  are  given.  The  curves  In 
figure  4  are  the  theoretical  estimates  of  the  strength  distribution  after 
proof  testing,  determined  from  the  crack  velocity  data  in  figure  lb. 

For  the  rapid  unloading- rate,  the  two  theoretical  curves  predict  sharply 

truncated,  strength  distributions.  Considering  this  prediction,  few  strengths 

should  have  been  measured  with  values  of  less  than  "  100  MPa.  However,  in 

contrast  to  this  expectation,  the  experimental  curve  was  not  sharply 

truncated;  approximately  16  percent  of  the  specimens  broke  at  stresses  less 

than  100  MPa.  Although  the  experimental  data  fall  along  the  theoretical 

curves  at  strength  levels  greater  than  100  MPa,  there  is  a  low-strength  "tail" 

in  the  data  that  is  not  explained  by  either  theoretical  curve.  This  "tail" 
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was  also  observed  by  Ritter  et  al.  in  their  recent  study  of  soda-lime  silica 
glass  proof  tested  in  heptane. 

The  strength  distribution  for  the  rapid  unloading  rate  (figure  4)  differs 
considerably  from  the  distribution  obtained  in  an  earlier  study  In  which  the 
proof  testing  was . conducted  in  dry  nitrogen  gas.11  In  the  earlier  study,  the 
experimental  data  closely  approximated  the  theoretical  curve.  The  difference 
in  results  between  studies  is  probably  not  attributable  to  water  In  the 
environment.  In  both  studies  the  unloading  rate  was  sufficiently  fast  that 


the  strength  degradation  during  the  proof  test  should  have  been  controlled  by 
the  "water  free"  portions  of  the  crack  growth  curves.  It  is  possible  that  the 
difference  in  results  is  due  to  physical  rather  than  chemical  aspects  of  the 
two  environments.  Viscosity,  for  example,  is  much  greater  for  liquids  than 
for  gases;  therefore  viscous  drag  by  the  liquid  on  the  crack  surfaces  could 
retard  the  motion  of  the  crack  and  permit  some  specimens  to  survive  which 
would  otherwise  have  failed  during  the  proof  test.  Clearly,  additional 
studies  are  needed  to  confirm  this  possibility,  and  more  generally,  to  explore 
differences  in  the  effect  of  liquids  and  gases  on  the  strength  of  materials. 
From  a  practical  point  of  view,  however,  results  from  the  present  study  in 
heptane  and  the  earlier  one  in  nitrogen  gas  suggest  that  proof  testing  in 
liquids  may  not  be  a  wise  procedure,  because  the  strength  distribution  curve 
may  not  be  sharply  truncated  when  the  test  is  performed  in  a  liquid. 

In  the  case  of  the  medium  and  slow  unloading  rates  (figure  4)  the 
theoretical  curves  depart  significantly  from  the  experimental  data:  the 
departure  is  much  greater  than  would  be  expected  from  the  scatter  of  the 
experimental  data.  Furthermore,  the  two  theoretical  curves  in  these  figures 
show  a  divergence  that  is  surprisingly  large  considering  the  small  difference 
between  the  crack  velocity  curves  in  figure  1.  These  differences  between  the 
experimental  data  and  the  theoretical  curve  on  the  one  hand  and  the  two 
theoretical  curves  on  the  other,  suggest  that  the  strength  distribution  curves 
after  proof  testing  are  extremely  sensitive  to  the  exact  location  of  the 
curves  that  describe  crack  motion.  As  can  be  seen  by  comparing  figure  1  with 
the  medium  and  slow  unloading  rate  curves  of  figure  4,  small  changes  in  the 
position  of  the  v-Kj  curve  result  in  very  large  differences  in  the  predicted 
position  of  the  strength  distribution  curve  after  proof  testing.  This  extreme 
sensitivity  of  the  strength  distribution  to  the  crack  velocity  also  explains 


the  differences  between  the  experimental  and  theoretical  strength  distribution 
curves  in  figure  4. 

Maximum  Crack  Velocity  During  Proof  Testing 

The  sensitivity  of  the  strength  distribution  curves  to  the  exact  shape 
and  position  of  the  crack  velocity  curves  can  be  understood  in  terms  of  the 
maximum  velocity  achieved  by  the  critical  flaw  in  the  course  of  the  proof 
test.  During  a  proof  test,  strength  loss  occurs  at  all  levels  of  applied 
stress  (above  the  fatigue  limit)  as  a  consequence  of  subcritical  crack  growth. 
The  strength  loss  is  most  severe  where  the  crack  velocity  is  highest.  Thus, 
the  high  velocity  regime  of  the  v-Kj  diagram  is  most  Important  for  determining 
the  distribution  of  strengths  after  proof  testing.  In  fact,  because  most 
crack  growth  occurs  at  the  maximum  crack  velocity,  strength  degradation  should 
be  related  to  the  maximum  velocity  achieved  in  the  course  of  the  proof  test: 
the  higher  the  velocity  the  greater  the  amount  of  strength  degradation. 

To  examine  the  relation  between  the  maximum  crack  velocity  during  proof 

testing  and  the  strength  distribution  curve,  the  maximum  crack  velocity  was 

calculated  for  each  of  the  points  used  to  plot  the  theoretical  curves  shown  in 

figure  4.  The  maximum  velocity  was  determined  and  then  printed  next  to  each 

point  on  the  theoretical  curve  to  Indicate  the  maximum  velocity  during  the 

proof  test.  An  example  of  the  type  of  data  generated  is  shown  in  figure  5  for 

the  slow  unloading  rate  (results  for  the  fast  and  medium  unloading  rates  are 

similar  to  figure  5..  and  are  not  discussed  here).  As  can  be  seen  from  figure 

5  the  range  of  maximum  velocities  was  relatively  narrow  for  both  sets  of 

curves,  covering  approximately  two  orders  of  magnitude  (10  ®  to  10  ®  m/s)  for 
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the  SEN  specimens  and  one  order  of  magnitude  (10  to  10  m/s)  for  the  DCB 
specimens.  Regardless  of  the  unloading  rate,  or  the  v-Kj  data,  the  narrow 
velocity  range  was  common  to  each  of  the  strength  distribution  curves  obtained 


in  the  present  study.  Thus  only  a  small  portion  of  the  v-Kj  plot  is  effective 
in  establishing  the  position  and  shape  of  the  strength  degradation  curve  after 
proof  testing.  This  point  is  illustrated  in  figure  6  which  shows  the  range  of 
crack  velocities  that  was  effective  in  establishing  the  strength  degradation 
curves  shown  in  figure  5.  Moving  the  v-Kj  plot  by  a  small  amount  dramatically 
changes  the  portion  of  the  v-Kj  curve  that  determines  the  strength 
distribution  after  proof  testing.  For  example,  in  figure  5  the  shape  of  the 
DCB  curve  is  determined  primarily  by  region  I  crack  growth  behavior.  Lowering 
the  region  II  portion  of  the  v-Kj  curve  slightly,  as  has  been  done  for  the  SEN 
specimen,  permits  a  large  fraction  of  the  specimens  to  achieve  velocities  that 
are  characteristic  of  region  II  crack  growth  behavior.  Because  the  slope  of 
the  strength  distribution  curve  depends  on  the  slope  of  the  crack  growth 
curve  ,  this  change  in  crack  growth  behavior  dramatically  alters  the  shape  of 
the  strength  distribution  curve  after  proof  testfng. 

The  general  shape  and  appearance  of  the  strength  distribution  curves  in 

figure  5  are  a  direct  consequence  of  the  fact  that  these  curves  were 

calculated  from  the  crack  velocity  data  in  figure  1.  Both  the  crack  growth 

and  strength  curves  are  multi-regioned.  Because  of  the  importance  of  the 

maximum  velocity  as  a  determinant  for  strength,  each  segment  of  the  strength 

distribution  curve  is  related  to  a  segment  of  the  crack  growth  curve. 

However,  the  two  curves  are  inverted  relative  to  one  another:  the  high 

velocity  region  of'the  crack  growth  curve  maps  onto  the  low  probability  region 

of  the  strength  distribution  curve  and  vice  versa.  Furthermore,  the  slopes  of 

corresponding  segments  of  the  two  curves  are  also  related:  the  slope,  m,  of 

the  strength  distribution  curve  is  given  by  n-2  where  n  is  the  slope  of  the 
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crack  velocity  curve.  This  one-to-one  relation  between  the  two  curves  plays 
an  important  role  in  determining  the  sensitivity  of  the  strength  distribution 


curve  to  the  crack  velocity  curve,  especially  when  only  a  narrow  range  of 
crack  velocity  data  Is  relevant  to  the  strength  distribution  curve,  as  is  the 
case  in  the  present  study* 

5.  CORRELATION  BETWEEN  THEORY  AND  EXPERIMENT 

In  view  of  the  above  discussion,  the  lack  of  correlation  between  the 
experimental  and  theoretical  strength  distributions  can  be  understood  in  terms 
of  the  high  sensitivity  of  the  strength  distribution  curve  to  the  exact  shape 
and  position  of  the  crack  growth  curve.  Any  physical  or  chemical  process  that 
disturbs  the  position  of  this  v-Kj  curve,  or  prevents  an  accurate 
determination  of  the  curve  can  affect  the  position  and  shape  of  the  strength 
distribution  curve,  and  thus  result  in  a  lack  of  agreement  between  the 


theoretical  and  experimental  strength  distribution  curves.  Some  of  the  causes 
of  poor  agreement  between  the  theoretical  and  experimental  distribution  curves 
will  now  be  discussed.  As  will  be  seen,  the  causes  are  of  two  types: 
experimental  errors  that  prevent  an  accurate  determination  of  the  "true"  v-Kj 
curve;  and  physical  or  chemical  effects  that  alter  the  behavior  of  small 
surface  cracks  in  such  a  way  that  their  response  to  applied  stress  is 
significantly  different  from  that  obtained  on  large  macroscopic  cracks.  Of 
the  two  causes,  the  latter  are  more  fundamental  to  our  understanding  of  crack 
growth  and  strength  degradation.  However,  in  practice  both  can  be  equally 
important  for  the  determination  of  accurate  strength  distribution  curves. 

Exper i mental  Uncertai nty 

The  experimental  uncertainty  to  be  associated  with  the  crack  growth  data 
is  illustrated  in.  figure  1,  by  the  small  difference  in  position  of  the  two 
v-Kj  curves  and  by  the  random  scatter  of  the  data  of  each  curve.  The 
difference  between  the  two  curves  (figure  la)  represents  a  systematic 
uncertainty  in  measurement  that  has  its  origin  in  the  technique  used  to  make 
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the  measurement.  This  type  of  error  can  be  reduced  only  by  more  thorough 
theoretical  and  experimental  investigations  of  the  two  crack  velocity 
techniques  to  trace  the  origin  of  the  measurement  difference.  The  random 
error  in  figure  1  is  given  by  the  scatter  of  data  about  their  mean  value  for 
any  given  value  of  Kj  or  v.  As  can  be  seen  from  this  figure  the  scatter  in  Kj 
ranges  from  3  to  10  percent  of  Kj,  whereas  the  error  in  v  ranges  from  ~  0.2 
to  ~  1  order  of  magnitude  in  crack  velocity,  depending  on  the  value  of  v  and 
Kj  and  the  technique  used  to  determine  the  v-Kj  plot.  These  errors  are 
typical  of  those  that  have  been  reported  previously  in  the  literature  for 
crack  velocity  studies  on  glass. 

With  regard  to  the  present  experiment,  the  experimental  scatter  of  the 
data  in  figure  1  is  sufficient  to  preclude  a  determination  of  the  “true"  v-Kj 
curve  with  sufficient  accuracy  to  predict  the  position  and  shape  of  the 
strength  distribution  curve  after  testing.  As  a  consequence  of  this  finding, 
it  is  suggested  that  more  accurate  measurement  techniques  are  needed  to 
determine  crack  growth  behavior.  Without  more  accurate  methods  of 
measurement,  one  should  be  wary  of  using  crack  growth  data  to  predict  proof 
testing  results.  However,  even  with  better  measurement  techniques  there  may 
be  more  basic  reasons  that  strength  distribution  curves  cannot  be  determined 
accurately  from  crack  growth  data.  These  reasons  relate  to  the  equivalence  of 
crack  growth  data  from  large  and  small  cracks. 

Large  versus  Small  Cracks 

Although  it  is  usual  to  assume  that  the  growth  of  microscopic  size  cracks 
in  ceramic  materials  is  described  by  crack  growth  data  obtained  by 
macroscopic,  fracture-mechanics  techniques,  there  are  good  physical  reasons 
why  this  assumption  may  not  be  valid.  In  other  words,  small  cracks  that 
control  the  strength  of  glass  may  not  behave  in  the  same  way  as  large  cracks 


108 


VA-.  A  j  v. 


J 


used  in  fracture  mechanics  studies.  Three  factors  that  may  contribute  to  the 
difference  in  crack  growth  behavior  are:  (1)  residual  stresses  due  to  plastic 
deformation  at  the  origin  of  surface  cracks;  (2)  blunting  of  small  cracks  due 
to  stress  corrosion;  and  (3)  differences  in  the  v-Kj  curve  resulting  from 
geometric  differences  between  large  and  small  cracks. 

Residual  Stresses 

The  effect  of  residual  stresses  on  the  strength  of  glass  was  first 

21  22 

discussed  by  Marshall  and  Lawn  *  .  The  effect  results  from  the  fact  that 

surface  cracks  in  glass  are  usually  formed  by  mechanical  damage  during  the 
grinding  and  polishing  of  glass,  or  by  mechanical  contact  with  hard  substances 
after  the  glass  has  been  produced.  Residual  stress  fields  always  form  at 
points  of  mechanical  contact  as  a  consequence  of  micro-plastic  deformation  at 
contact  asperities.  Unless  the  residual  stresses  are  relieved  by  annealing, 
they  contribute  to  the  stress  field  that  forces  cracks  to  grow  during  a  proof 
test.  Therefore,  these  stresses  should  be  taken  into  account  to  determine  the 
strength  distribution  after  proof  testing.  As  the  crack  grows,  the  effect  of 
residual  stresses  gradually  becomes  less  Important;  the  effect  is  least  for 
the  weakest  specimens  that  survive  the  proof  test.  Hence,  the  effect  of 
residual  stresses  on  the  strength  distribution  is  expected  to  be  most 
pronounced  for  the  high  probability  portion  of  the  stress  distribution.  In 
this  paper  it  is  assumed  that  residual  stresses  played  little  role  in  the 
strength  degradation  process.  This  assumption  is  supported  by  the  fact  that 
the  proof  test  results  in  figure  3  are  very  similar  to  those  obtained  by 
Ritter  et  al.  on  annealed  specimens  .  However,  to  fully  understand  the  role 
of  residual  stresses  on  proof  testing,  additional  studies  are  clearly  needed. 
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Crack  Blunting 

By  applying  transition  rate  theory  to  ceramic  materials,  Charles  and 

Hillig  suggested  the  existence  of  a  static  fatigue  limit  caused  by  crack  tip 

blunting  In  materials  such  as  glass.  At  stresses  above  the  fatigue  limit, 

surface  flaws  sharpen  and  grow,  whereas  at  stresses  below  the  fatigue  limit, 

crack  tip  blunting  and  crack  arrest  occurs.  The  Charles-HIllig  theory  Is 

consistent  with  all  known  experimental  data  on  the  strength  of  glass. 

Recently,  Michalske,  using  fracture  mechanics  techniques,  has  demonstrated  a 

fatigue  limit  In  soda-lime  silica  glass  at  a  value  of  Kj  “  0.25  MPa-m1^.** 

The  fatigue  limit  was  Identified  by  a  hysteresis  In  the  crack  velocity  curve, 

and  by  mlcrostructural  features  that  were  left  on  the  face  of  a  crack  that  was 

repropagated  after  having  been  blunted  by  stress  corrosion.  Resharpening  of 

cracks  In  soda-lime  silica  glass  occurs  almost  Instantaneously  when  the 

1/2 

stress- Intensity  factor  exceeds  0.425  MPa-m  .  If  the  stress  Intensity  for 
resharpening  falls  within  the  range  of  Kj  that  Is  effective  In  determining  the 
strength  distribution  after  proof  testing,  then  crack  blunting  can  have  a 
significant  effect  on  the  shape  of  the  strength  distribution  curve  after  proof 
testing.  Our  studies  on  this  point  are  quite  limited;  however,  they  do  show 
that  crack  blunting  is  more  Important  for  proof  tests  that  were  conducted 
using  slow  unloading  rates  rather  than  fast  ones.  In  addition  our  studies 
show  that  the  strongest  specimens  are  not  affected  at  all  by  the  proof  test, 
whereas  the  weakest  ones  suffer  the  same  degree  of  strength  degradation  as 
would  occur  if  there  were  no  static  fatigue  limit.  Additional  studies  to 
clarify  the  role  of  crack  blunting  on  proof  testing  should  include  the 
establishment  of  resharpening  kinetics  and  the  incorporation  of  these  kinetics 
into  the  strength  treatment. 
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difference  In  crack  growth  behavior  are:  (1)  residual  stresses  due  to  plastic 
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geometric  differences  between  large  and  small  cracks. 

Residual  Stresses 

The  effect  of  residual  stresses  on  the  strength  of  glass  was  first 
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discussed  by  Marshall  and  Lawn  *  .  The  effect  results  from  the  fact  that 

surface  cracks  in  glass  are  usually  formed  by  mechanical  damage  during  the 
grinding  and  polishing  of  glass,  or  by  mechanical  contact  with  hard  substances 
after  the  glass  has  been  produced.  Residual  stress  fields  always  form  at 
points  of  mechanical  contact  as  a  consequence  of  micro-plastic  deformation  at 
contact  asperities.  Unless  the  residual  stresses  are  relieved  by  annealing, 
they  contribute  to  the  stress  field  that  forces  cracks  to  grow  during  a  proof 
test.  Therefore,  these  stresses  should  be  taken  into  account  to  determine  the 
strength  distribution  after  proof  testing.  As  the  crack  grows,  the  effect  of 
residual  stresses  gradually  becomes  less  Important;  the  effect  Is  least  for 
the  weakest  specimens  that  survive  the  proof  test.  Hence,  the  effect  of 
residual  stresses  on  the  strength  distribution  is  expected  to  be  most 
pronounced  for  the  high  probability  portion  of  the  stress  distribution.  In 
this  paper  it  Is  assumed  that  residual  stresses  played  little  role  In  the 
strength  degradation  process.  This  assumption  is  supported  by  the  fact  that 
the  proof  test  results  in  figure  3  are  very  similar  to  those  obtained  by 
Ritter  et  al.  on  annealed  specimens  .  However,  to  fully  understand  the  role 
of  residual  stresses  on  proof  testing,  additional  studies  are  clearly  needed. 
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Crack  Blunting 

By  applying  transition  rate  theory  to  ceramic  materials,  Charles  and 
23 

Hillig  suggested  the  existence  of  a  static  fatigue  limit  caused  by  crack  tip 

blunting  in  materials  such  as  glass.  At  stresses  above  the  fatigue  limit, 

surface  flaws  sharpen  and  grow,  whereas  at  stresses  below  the  fatigue  limit, 

crack  tip  blunting  and  crack  arrest  occurs.  The  Chari es-Hi 11 ig  theory  is 

consistent  with  all  known  experimental  data  on  the  strength  of  glass. 

Recently,  Michalske,  using  fracture  mechanics  techniques,  has  demonstrated  a 

fatigue  limit  in  soda-lime  silica  glass  at  a  value  of  Kj  ~  0.25  MPa-m^.^ 

The  fatigue  limit  was  identified  by  a  hysteresis  in  the  crack  velocity  curve, 

and  by  microstructural  features  that  were  left  on  the  face  of  a  crack  that  was 

repropagated  after  having  been  blunted  by  stress  corrosion.  Resharpening  of 

cracks  in  soda-lime  silica  glass  occurs  almost  instantaneously  when  the 
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stress- intensity  factor  exceeds  0.425  MPa-m  .  If  the  stress  intensity  for 
resharpening  falls  within  the  range  of  Kj  that  is  effective  in  determining  the 
strength  distribution  after  proof  testing,  then  crack  blunting  can  have  a 
significant  effect  on  the  shape  of  the  strength  distribution  curve  after  proof 
testing.  Our  studies  on  this  point  are  quite  limited;  however,  they  do  show 
that  crack  blunting  is  more  important  for  proof  tests  that  were  conducted 
using  slow  unloading  rates  rather  than  fast  ones.  In  addition  our  studies 
show  that  the  strongest  specimens  are  not  affected  at  all  by  the  proof  test, 
whereas  the  weakest  ones  suffer  the  same  degree  of  strength  degradation  as 
would  occur  if  there  were  no  static  fatigue  limit.  Additional  studies  to 
clarify  the  role  of  crack  blunting  on  proof  testing  should  include  the 
establishment  of  resharpening  kinetics  and  the  incorporation  of  these  kinetics 
into  the  strength  treatment. 


Crack  Geometry 

Effects  of  environment  on  both  the  strength  and  the  crack  growth  rate 

have  been  studied  by  a  number  of  Investigators.  Studies  on  soda-1 ime-si llca 

glass  Indicate  relatively  good  agreement  between  strength  measurements  and 

crack  growth  measurements,  when  measurements  are  conducted  either  In  water 

where  region  I  crack  growth  is  dominant,  or  in  an  inert  environment  such  as 

dry  nitrogen,  where  region  III  crack  growth  is  dominant.  For  both  techniques 

the  stress  corrosion  susceptibility,  n,  was  found  to  be  similar  regardless  of 

* 

the  technique  used  to  make  the  measurement.  Therefore  it  is  unlikely  that 

the  behavior  of  large  cracks  differs  significantly  from  small  ones  when  only 

2S 

region  I  or  region  III  crack  growth  is  involved  in  degrading  the  strength. 

The  same  conclusion,  however,  may  not  be  valid  when  region  II  crack  growth 
occurs. 

Crack  growth  in  region  II  has  been  studied  by  a  number  of  investigators. 
26  •  27 

Wiederhorn  and  Schonert  et  al.  ,  were  the  first  to  attribute  region  II 
crack  growth  to  an  environment  limited  transport  of  water  to  the  crack  tip. 
Schonert  et  al.  demonstrated  that  the  crack  front  changed  its  shape  as  the 
crack  moved  from  region  I  to  region  II,  and  again  when  it  moved  from  region  II 
to  region  III.  When  studies  are  conducted  with  through  cracks  in  flat 
specimens  (edge-notched  or  double  cantilever  beam  specimens),  the  portion  of 
the  crack  front  in  the  center  of  the  specimen  always  leads  the  portion  of  the 


* 

When  residual  stresses  effect  the  crack  growth  behavior,  a  correction  must  be 

made  for  the  effect  of  these  stresses  on  the  stress  corrosion  susceptibility, 

n.  However,  even  in  this  case,  the  value  of  n  for  crack  growth  can  be 

22  23 

related  to  that  obtained  from  fatigue  data.  * 


crack  front  that  intersects  the  specimen  surface  in  region  I  or  III.  In 

region  II,  however,  the  reverse  is  true:  the  portion  of  the  crack  front  near 

the  surface  always  leads  that  in  the  center. 

28  29 

Quackenbush  and  Frechette  and  Richter  confirmed  the  findings  of 

Schonert  et  al.  and  attributed  crack  growth  in  region  II  to  a  combined  effect 

of  region  I  and  region  III  crack  growth  behavior.  Since  water  is  freely 

available  at  the  specimen  surface,  crack  motion  is  enhanced  by  water;  whereas, 

away  from  the  specimen  surface  the  crack  tip  is  starved  for  water  because  of 

the  longer  diffusion  path,  and  crack  motion  is  retarded.  In  specimens  with 

through  cracks,  the  cracks  move  as  a  unit  In  region  II,  the  center  portion  of 

the  crack  being  pulled  along  by  portions  near  the  surface.  As  the  applied 

stress-intensity  factor  in  region  II  is  increased,  the  crack  velocity  is 

observed  to  increase,  but  not  as  rapidly  as  in 'regions  I  or  III.  Furthermore 

as  Kj  is  increased,  more  of  the  crack  perimeter  behaves  as  if  it  were  depleted 

of  water,  and  finally  when  the  outer  portion  of  the  crack  is  completely 

depleted  of  water,  it  too  enters  region  III  and  the  crack  perimeter  regains 

its  normal  shape.  Thus,  region  II  represents  a  transition  in  crack  growth 

behavior  between  regions  I  and  III;  depending  on  the  availability  of  water  at 

the  crack  tip,  different  portions  of  the  crack  front  behave  somewhat 

* 

independently  as  the  crack  moves. 
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In  an  earlier  paper,  Varner  and  Frechette  suggest  that  crack  arrest  occurs 

periodically  in  the  central  portion  of  the  crack  in  region  II  (for  studies 

conducted  in  nitrogen  gas).  A  later  study  by  Richter  (conducted  in  air) 

and  by  Quackenbush  and  Frechette  (conducted  in  decane)  gave  no  indication 

of  crack  arrest  in  region  II.  The  results  of  the  present  study  on 

Indentation  crack  also  give  no  indication  of  crack  arrest  in  region  II. 


The  effects  just  described  also  occur  for  surface  cracks  that  are 
responsible  for  the  failure  of  strength  specimens.  The  effect  is  shown  in 
figure  7  which  illustrates  successive  positions  of  a  surface  crack  that  was 
propagated  to  failure  in  normal  heptane.  The  crack  was  formed  by  the  use  of  a 
diamond  pyramidal  indenter;  the  velocity  was  measured  by  the  use  of  a  sonic 
pulse  generator.  The  crack  position  is  indicated  by  ripple  marks  on  the 
fracture  surface,  each  set  of  marks  indicating  successive  positions  of  the 
crack  as  it  grows  larger.  In  figure  7a  the  ripple  marks  closest  to  the 
indentation  (half-way  between  the  two  arrows)  represent  region  I  crack  growth, 
whereas  those  furthest  from  the  indentation  represent  region  III  crack  growth. 
The  arrows  mark  the  transition  between  region  II  and  region  III  crack  growth. 
As  can  be  seen  from  figure  7b,  the  shape  of  the  crack  in  region  II  differs 
from  the  normal  semi -elliptical  shape  that  is  characteristic  of  region  I  or 
region  III  crack  motion.  As  in  the  case  of  the  through  cracks  discussed 
above,  the  change  in  shape  can  be  attributed  primarily  to  water  depletion  at 
portions  of  the  crack  that  are  distant  from  the  surface.  Because  of  water 
depletion,  these  portions  of  the  crack  move  slower  than  those  near  surface 
where  the  water  is  plentiful.  Thus,  the  shape  of  the  surface  crack  is 
modified  by  the  water  content  of  the  environment  near  the  crack  perimeter. 
Because  of  the  difference  in  geometry  between  small  surface  cracks  and  the 
large  cracks  that  are  typical  of  fracture  mechanics  specimens,  physical 

* 

constraints  on  the  j;rack  front  are  different  for  the  two  types  of  cracks  . 
Therefore,  crack  growth  in  region  II  will  probably  not  be  the  same  for  the 
two  types  of  cracks  even  though  crack  growth  in  region  II  represents  a 

* 

For  example,  Kj  for  an  elliptical  crack  can  vary  by  as  much  as  50  percent 

along  the  crack  perimeter. 


transition  in  behavior  between  regions  I  and  III  for  both  crack  types.  The 
differences  between  the  theoretical  and  experimental  strength  distribution 
curves  in  figure  4  may  be  the  result  of  these  geometric  considerations. 
Substantiation  of  this  belief,  however,  will  require  the  solution  of  three 
dimensional  elasticity/diffusion  boundary  value  problems  to  determine  exactly 
how  crack  growth  occurs  for  surface  cracks.  This  task  is  beyond  the  scope  of 
the  present  paper. 

6.  SUMMARY 

This  paper  presents  a  critical  analysis  of  the  applicability  of  crack 
growth  data  to  proof  testing  as  a  means  of  improving  the  reliability  of 
structural  ceramic  materials.  In  particular,  the  paper  deals  with  the 
question  of  predicting  strength  distribution  curves  from  crack  growth  data 
obtained  by  standard  fracture-mechanics  techniques.  To  deal  with  this 
question,  both  strength  data  and  crack  growth  data  were  obtained  on 
soda- lime- silica  glass  specimens  that  were  tested  in  liquid  heptane.  Crack 
growth  data  were  obtained  by  the  applied-moment,  double-cantilever  beam  (DCB) 
technique,  and  by  the  single-edge,  notch  tension  (SEN)  technique.  Although 
there  were  systematic  differences  between  the  two  sets  of  crack  velocity  data, 
a  partial  overlap  of  the  data  was  observed  over  most  of  the  data  range. 
Considering  that  two  different  test  techniques  were  used  to  collect  the  data, 
the  two  sets  of  crack  growth  data  were  considered  to  be  consistent  with  the 
degree  of  scatter  tjjat  has  been  reported  for  crack  growth  data. 

Proof-test  experiments  were  conducted  on  soda-lime  silica  glass 
microscope  slides,  that  were  subjected  to  4-point  bending.  Strength 
distribution  curves  after  proof  testing  were  obtained  using  three  different 
unloading  rates  during  the  proof  test  cycle.  Truncation  of  the  strength  data 
was  obtained  only  for  specimens  subjected  to  a  rapid  (165  MPa/s)  rate  of 


unloading.  Specimens  that  were  subjected  to  either  slow  (0.33  MPa/s),  or 
medium  (3.3  MPa/s)  rates  of  unloading  gave  curves  that  hardly  differed  from 
the  initial  strength  distribution  curve.  Furthermore,  the  three  strength 
distribution  curves  obtained  from  the  proof-tested  population  of  glass  slides 
could  not  be  predicted  from  the  crack  growth  data  obtained  in  the  present 
experiment. 

In  the  course  of  our  discussion  and  analysis,  we  concluded  that  the 
position  and  shape  of  the  strength  distribution  curve  after  proof  testing  was 
highly  sensitive  to  the  exact  position  and  shape  of  the  crack  growth  data.  We 
were  able  to  show  that  slight  shifts  in  the  position,  or  shape  of  the  v-Kj 
plot  resulted  in  substantial  changes  in  the  position  and  shape  of  the 
predicted  strength  distribution  curve.  As  a  consequence,  any  physical  or 
chemical  process  that  disturbs  or  changes  the 'position  of  the  crack  growth 
curve,  or  prevents  an  accurate  determination  of  the  crack  growth  curve  will 
cause  the  strength  distribution  curve  determined  from  crack  growth  data  to 
differ  considerably  from  the  one  determined  experimentally.  Specific  causes 
of  the  differences  are  (1)  the  experimental  uncertainty  inherent  in  obtaining 
v-Kj  data,  and  (2)  physical  and  chemical  effects  that  alter  the  behavior  of 
small  surface  cracks  so  that  their  response  to  applied  stress  differs 
significantly  from  the  response  of  large  macroscopic  cracks  typical  of 
fracture  mechanics  specimens. 

From  a  practical  point  of  view  our  strength  data  suggests  that  proof 
testing  in  an  "inert"  liquid  may  not  be  a  wise  procedure  because  only  a  modest 
degree  of  truncation  is  observed  (only  at  high  unloading  rates)  and  because 
after  proof- testing  many  specimens  have  strengths  that  lie  below  the  proof 
test  level.  Based  on  these  observations  and  results  from  an  earlier  study, 
which  included  proof  test  studies  in  dry  nitrogen  gas,  we  suggest  that  proof 


testing  be  conducted  in  dry,  gaseous  environments,  and  that  rapid  rates  of 
unloading  be  used  to  assure  structural  reliability.  When  this  procedure  is 
followed,  satisfactory  truncation  of  the  strength  distribution  can  be 
achieved. 
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TABLE  1.  TEST  CONDITIONS  FOR  SODA- LIME  SILICA  GLASS  PROOF  TESTED  IN  HEPTANE 


LOADING 

UNLOADING 

PROOF 

NUMBER 

PROOF 

RATE 

RATE 

LOAD 

TESTED 

TEST 

(MPa/s) 

(MPa/s) 

(MPa) 

SURVIVORS 

33 

STRENGTH 

TESTS 

628 

33 

165 

103 

673 

417 

33 

3.3 

100.5 

630 

179 

33 

0.33 

90.5 

497 

143 

* 

All  strength  measurements  were  made  by  4-polnt  bending  using  an  Inner  span 
of  19  mm  and  an  outer  span  of  64  mm. 


UNLOADING 

RATE 

(MPa/s) 


TABLE  2.  TEST  OF  WEIBULL  DISTRIBUTION 


Sq  =  112  MPa,  M=7 


PROOF  FRACTION  BROKEN 
LOAD  DURING  LOADING 
(MPa) 


MEASURED  PREDICTED 


APPENDIX  A 


The  Selection  of  Kjc 

In  developing  the  computer  code  used  in  the  present  study,  we  found  it 

necessary  to  establish  a  criterion  for  failure.  A  standard,  albeit  somewhat 

idealized,  failure  criterion  is  to  assume  that  failure  occurs  when  the 

31 

strength,  S,  is  equal  to  the  applied  stress,  a.  In  fracture  mechanics 
terms,  this  criterion  is  equivalent  to  assuming  KI  *  KIc  and  d  Kj/da  >  0  for 
the  most  critical  crack. 

In  more  practical  terms,  however,  failure  by  crack  growth  is  determined 
by  a  crack  velocity  criterion:  failure  occurs  when  the  crack  is  moving  so 
rapidly  that  the  specimen  cannot  be  unloaded  fast  enough  to  avoid  failure. 

The  critical  velocity  for  failure  depends  on  the  rate  of  unloading:  as  the 
unloading  rate  is  increased,  the  critical  velocity  also  is  increased.  This 
failure  criterion  is  exactly  applicable  to  proof-testing,  since  the  number  of 
specimens  that  pass  the  proof  test  Is  determined  by  the  rate  of  unloading. 

When  the  v-Kj  curve  is  steep  as  it  is  for  soda- lime  silica  glass  in 
nitrogen  gas,  the  two  criteria  are  almost  identical.  The  value  of  Kjc  is 
insensitive  to  the  crack  velocity,  and  the  value  selected  (within  region  III) 
for  Kjc  does  not  affect  the  truncation  that  results  from  proof  testing. 
However,  when  the  v-,Kj  curve  is  complicated  by  physical  effects,  such  as 
viscous  drag  or  cavitation,  as  it  is  for  soda-lime  silica  glass  tested  In 
heptane,  then  one.  must  be  concerned  with  the  equivalence  of  the  two  failure 
criteria.  In  particular,  it  must  be  determined  whether  KIc  determined  by 
"conventional"  methods  can  be  used  to  establish  a  criterion  for  crack 


stability  for  a  proof  test. 


In  the  present  study  we  used  a  value  of  KIC  =  0.  75  MPaVm  as  obtained  on 
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DCS  specimens  that  were  loaded  to  failure  in  dry  nitrogen  gas.  To  test  the 

validity  of  this  convenitonal  value  of  Kjc,  the  v-Kj  curve  for  the  SEN 

specimen,  figure  lb,  was  used  to  determine  a  strength  distribution  curve  for 

three  values  of  "Kjc":  0.75,  0.90  and  1.10  MPa-m^.  For  each  of  the  proof 

test  cycles  used  in  the  present  study,  the  strength  distribution  after  proof 

testing  and  the  number  of  specimens  that  broke  during  the  proof  test  were 

found  to  be  independent  of  the  value  of  "Kjc"  used  in  the  program.  The  reason 

for  this  apparent  insensitivity  of  the  strength  distribution  curve  to  Kjc  is 

the  fact  that  the  distribution  curve  is  determined  by  a  relatively  narrow 

range  of  velocities  on  the  v-Kj  plot.  As  long  as  "Kjc"  lies  at  a  value  that 

is  above  this  range  of  velocities,  the  strength  curve  will  not  be  sensitive  to 

its  value.  The  main  effect  of  changing  "Kjc*'  is- to  shift  the  position  of  the 

* 

initial  strength  distribution  .  Increasing  the  value  of  "Kjc"  increases  the 
value  of  SQ  and  decreases  the  value  of  m  for  the  initial  distribution.  Thus, 
"Kjc"  plays  the  role  of  a  scaling  factor  in  a  proof  test.  In  view  of  this 
conclusion,  the  fracture  mechanics  value  of  Kjc  (i.e.  0.75  MPa-m  )  was  used 
for  the  present  study. 

*The  initial  distribution  referred  to  here  is  the  distribution  calculated  from 
the  input  data  discussed  in  Section  2. 


TABLE  la.  INPUT  DATA:  WEIBULL  PARAMETERS 


Condition  SQ  ■ 

Experimental  112  7 

Inert:  calculated  from 
v-Kj  data 

SEN  123.1  6.55 

DCB  130.1  6.28 


■>v; 


TABLE  2a.  CRACK  GROWTH  DATA 

0.75  MPa-m172 
0.2  MPa-ro1/2 
1.5 

Slopes  and  Intercepts  for  v-Kj  curves;  values  based  on  Pa  as  units  for  stress 


Condition 

Region 

Slope,  n 

In  A 

SEN 

1 

16.59 

-  233.13 

2 

5.07 

-  81.90 

3 

245.66 

-3315.52 

4 

56.00 

-  762.88 

5 

13.02 

-  180.91 

6 

6.97 

-  98.23 

7 

382.50 

-5255.06 

8 

16.03 

-  218.66 

9 

6.60  - 

-  88.35 

DCB 

1 

17.10 

-  237.95 

2 

4.09 

-  67.85 

3 

125.83 

-1700.19 

Note:  The  number  of  significant  figures  given  in  the  above  parameters  are 

Included  for  purposes  of  calculation  only,  and  do  not  indicate  accuracy 
of  the  various  parameters. 


0.5  0.7  04  1.1 


Omm  Int amity  Factor,  K-  MPm1'1 


04  0.4  0.5  0.6  0.7  04  0414  14 

Stmt  Intensity  Factor,  K|,  MPa-m* 


Crack  velocity  data,  soda-lime  silica  glass  In  heptane,  (a)  Crack 
velocity  data  collected  by  the  applied  moment,  double  cantilever 
techniques  and  the  single-edge,  notch  technique,  (b)  Piecewise, 
straight-line  fit  of  the  crack  growth  data  shown  In  (a).  Curve  marked 
SEN  represents  the  single  edge,  notch  data;  curved  marked  DCB  represents 
the  double  cantilever  data. 
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100  140  180 


60  100  140  180 


Strength,  MPa 

Strength  distribution  after  proof  testing  in  heptane:  (a)  rapid  (165 
MPa/s)  unloading  rate;  (b)  medium  (3.3  MPa/s)  unloading  rate;  (c)  slow 
(0.33  MPa/s)  unloading  rate. 
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0.99 
=  0.90 
|  0.60 
|  0.30 
%  0.10 


60  80  100  120  140  160 

Strength,  MPa 


Maximum  velocities  (m/s)  achieved  during  proof  testing  printed  next  to 
calculated  strength  distribution  curves.  The  maximum  crack  velocities 
and  the  calculated  strength  distribution  curves  were  obtained  from  the 

v-Kj  data  shown  in  figure  1:  slow  (0.33  MPa/s)  unloading  rate. 

* 


1q2  k  Slow  Unloading  Rata 


10° 

1 

>  10"a 


1  10-* 


2  10-* 


10-10 


0.3  0.4  0J5  0.6  0.7  0.8  0.9 1.0  1.2 


Strata  (ntanslty  Factor,  Kb  MPa-mH 


Range  of  crack  velocities  that  determine  the  strength  distribution  curves 
for  the  slow  (0.33  MPa/s)  unloading  rate. 


** 

Successive  positions  of  a  surface  crack  that  was  propagated  to  failure  in 
heptane:  (a)  arrows  mark  the  transition  between  region  II  and  region  III 
crack  growth;  (b)  higher  magnification  of  (a)  showing  the  shape  of  the 
crack  fronts  in  region  II. 


To  be  published  In  Ferroelectrlcs,  1984. 


FRACTURE  OF  FERROELECTRIC  CERAMICS 


R.F.  COOK,  S.W.  F REIMAN,  B.R.  LAWN 

Inorganic  Materials  Division,  National  Bureau  of  Standards, 
Washington.  DC  20234 

R.C.  POHANKA 

Office  of  Naval  Research,  Code  431,  Arlington,  Virginia 

Abstract  This  paper  surveys  the  temperature,  microstructural 
and  environmental  variations  of  the  fracture  properties  of 
ferroelectric  ceramics.  Earlier  work  shows  that  fracture 
toughness  decreases  on  heating  through  the  Curie  temperature. 
There  is  also  anomalous  behavior  in  the  strength  at  small 
crack  sizes,  indicative  of  a  grain  size  effect.  Further,  the 
strength  properties  are  known  to  be  adversely  affected  by  the 
presence  of  water  in  the  atmosphere.  Data  from  recent 
indentation  studies  on  barium  titanate  are  used  to 
investigate  these  phenomena. 


INTRODUCTION 


Ferroelectric  ceramics  such  as  barium  titanate  (BaTiO^)  and  lead 
zirconate  titanate  (PZT)  are  susceptible  to  brittle  failure  from 
small  processing  and  handling  flaws.  To  characterize  the  strength 
properties  of  brittle  materials  we  need  to  understand  the  nature 
of  all  forces  acting  on  the  flaws.  We  shall  address  these  issues 
in  terms  of  data  from  indentation  fracture  studies. 


INDENTATION  TESTING 


Reproducibility  in  strength  data  may  be  optimized  by  introducing 
controlled  indentation  flaws  Into  prospective  test  specimens.* 
The  driving  force  on  such  flaws  comes  from  two  main  sources,  the 
externally  applied  stress  field  and  a  local  residual  field  about 
the  elastic/plastic  impression.  Under  equilibrium  fracture 
conditions  the  strength  is  given  by 
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AK  4/3/P1/3 
c 


(1) 


where  A  is  a  dimensionless  constant,  K  is  the  material  toughness 
and  P  is  the  original  indentation  load.  Departures  from  the 
predicted  variation  of  o  on  P  in  Eq.  (1)  could  be  effected  either 
by  a  third,  unaccounted  driving  force  on  the  indentation  flaw,  or 
by  a  systematic  dependence  of  K^  on  crack  size. 

Accordingly,  controlled  flaw  tests  were  run  on  BaTiO^. 
Indentations  were  made  on  annealed  bar  specimens  with  a  Vickers 
diamond  pyramid,  in  air,  at  prescribed  loads.  Failure  was 
produced  in  four-point  flexure  in  a  controlled  environment.  At 
low  loads,  failure  tended  to  initiate  from  natural  rather  than 
introduced  flaws,  thereby  imposing  a  lower  limit  on  the  data 
range . 


TEMPERATURE  EFFECTS 


Strength  variations  were  followed  as  a  function  of  temperature  for 

a  BaTiO,  material  of  nominal  grain  size  7  pm.  The  indentation 

load  for  these  tests  was  fixed  at  30  N,  so  that  the  cracks 

produced  were  well  in  excess  of  the  grain  size.  The  flaws  were 

introduced  at  room  temperature,  and  the  strength  tests  were  run  in 

a  heated  oil  bath.  Figure  1  shows  the  results.  There  is  a 

monotonic  decrease  in  strength  to  the  vicinity  of  the  Curie 

temperature,  T  ,  then  an  apparent  levelling  out.  The  falloff 

between  23°C  and  I50°C  corresponds  to  a  reduction  of  in 

(see  Eq.  1),  consistent  with  trends  noted  previously  in 

2  3 

monocrystalline  BaTiO^.  ’  Specimens  tested  after  cycling  through 
Tfi  showed  the  same  strengths,  within  experimental  scatter,  as 
those  heated  directly  from  room  temperature. 

In  view  of  this  reversibility,  it  can  be  concluded  that  the 
results  in  Fig.  1  reflect  intrinsic  thermal  effects  in  the 
toughness  parameter.  The  reported  appearance  of  domain  wall 
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FIGURE  I  Strength  ea  function  of  temperature. 

markings  on  fracture  surfaces  produced  below  Tfi  suggests  that 

crack/ twin  interactions  could  account  for  the  trends  observed 

2 

here  ;  the  falloff  of  K  with  temperature  would  then  be 

c 

attributable  to  the  thermally  activated  annihilation  of  domain 

4 

walls  as  the  transition  is  approached. 

5.6 

Analogous  thermal  effects  have  been  reported  for  PZT. 
MICROSTRUCTURAL  EFFECTS 

The  Influence  of  microstructure  on  strength  properties  was 

investigated  by  systematically  reducing  the  controlled-f law  size. 

via  the  Indentation  load.  Data  are  plotted  in  Fig.  2  for  the  7  urn 

(Fig.  I)  and  a  1  ym  BaTlO^.  Several  points  may  be  noted: 

(1)  At  I50*C  the  data  fit  the  dependence  predicted  by  Eq.  (1)  over 

the  entire  load  range,  for  both  materials.  Moreover,  the 

representative  straight  lines  ars  Indistinguishable,  indicating 

that  K  is  the  same  in  both  cases.  This  grain-sise  invariance  of 
c 

toughness  above  the  Curie  point  is  consistent  with  data  on  other 
cubic  ceramics. ^ 
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FIGURE  2  Strengch  as  function  of  load. 

(ii)  At  25°C  the  data  are  well  behaved  over  only  part  of  the  load 

range.  In  this  well-behaved  region  there  is  overlap  with  the  1$0“C 

data  for  the  I  pm  material,  whereas  for  the  7  pm  material  the 

temperature  dependence  noted  in  the  strength  data  of  Fig.  1  is 

apparent.  These  results  imply  a  grain-size  dependence  in  K£  below 

the  Curie  point.  A  systematic  study  of  this  dependence  in  BaTiO^ 

materials  over  a  grain-size  range  1-150  pm  shows  that  K  can  vary 

2  C 

by  as  much  as  a  factor  of  two. 

(iii)  At  low  loads  the  25°C  data  for  both  materials  drop  off 
below  the  theoretical  line.  As  alluded  earlier,  such  a  departure 
could  be  due  either  to  the  presence  of  an  additional  driving  force 
on  the  cracks  or  to  a  tendency  for  K^  to  increase  in  a  systematic 
manner  as  the  cracks  intersect  a  greater  number  of  grains.  Of 
these  two  possibilities,  only  the  first  appears  to  be  capable  of 
accounting  for  the  reversion  to  ideal  behavior  at  the  elevated 
temperature,  since  the  grain  geometry  remains  invariant  on 
traversing  the  Curie  point.  Thus  it  has  been  suggested  that  the 
room  temperature  data  are  consistent  with  the  existence  of  local 
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internal  stresses  associated  with  a  non-cubic  (tetragonal) 

2  4 

polycrystal.  '  These  stresses  are  felt  more  strongly  as  the 
cracks  become  smaller,  approaching  the  scale  of  the  micro¬ 
structure;  in  this  context,  it  is  noted  in  Fig.  2  that  the 
deviation  occurs  at  a  lower  level  of  P  in  the  1  pm  material. 

(iv)  At  high  loads  the  258C  data  for  the  1  pm  material  again 
deviates  from  ideal  behavior.  The  cause  of  this  anomaly  is  not 
understood  at  present,  although  analogous  behavior  has  recently 

g 

been  reported  in  a  La-doped  lead  titanate. 

ENVIRONMENTAL  EFFECTS 

The  access  of  water  to  crack  tips  in  brittle  ceramics  can  cause 
subcritical  growth,  i.e.  growth  at  stress  levels  less  than  that 
needed  to  maintain  the  equilibrium  configuration  implied  in 
Eq.  (1).  The  rate  of  subcritical  crack  growth  increases 
dramatically  with  stress  intensity,  approximately  to  some  power  N, 
where  N  >>  1  typically.  Accordingly,  specimens  loaded  at  low 
stressing  rates  are  subject  to  more  extensive  subcritical  growth, 
and  hence  to  a  greater  degradation  in  strength.  For  specimens 
with  Indentation  flaws,  the  "fatigue  susceptibility"  parameter  N 
can  be  measured  directly  from  the  (inverse)  slope  of  a  strength  vs 

9 

stressing  rate  plot. 


FIGURE  3  Strength  as  function  of  stressing  rate. 
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Accordingly,  such  a  plot  is  shown  in  Fig.  3  for  the  7  pm 
BaTiO^  tested  in  water  at  25°C.  The  strength  Indeed  falls  off 
significantly  as  the  stressing  rate  decreases  (or,  alternatively, 
as  the  time  to  failure  increases) .  The  value  N  ■  67  obtained  from 
these  data  is  considerably  higher  than  that  for  soda-lime  glass 
(N  =  18,  representing  one  of  the  most  susceptible  of  all  brittle 
materials) ,  but  is  within  the  range  found  for  a  broad  spectrum  of 
ferroelectric  ceramics.  Small  changes  in  composition, 
particularly  in  the  form  of  grain  boundary  impurities,  could  well 
give  rise  to  deleterious  decreases  in  the  N  parameter. 
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The  purpose  of  this  paper  is  twofold*  First  it  will  show  to  what  extent  the 
fracture  of  pisxoslsctric  materials  such  as  BaTi03  are  influenced  by 
microstructure  end  the  phase  transformation  from  the  paraelactrlc  to  the 
ferroelectric  stete*  It  will  be  demonstrated  that  multiple  toughening 
mechanisms  such  ee  twinning  and  mlcroeracklng  account  for  the  observed 
behavior  and  that  fracture  mechanics  principles  can  be  used  to  develop 
materials  having  an  optimum  resistance  to  cracking*  Second,  it  will  show  that 
subcrltlcal  crack  growth  in  capacitors  will  ultimately  lead  to  failure* 

Again,  fracture  mechanics  techniques  will  be  used  to  make  predictions  as  to 
the  likely  failure  times  for  actual  capacitors* 

FRACTURE  OF  BaTi03 

Pohanka  and  co-workers  (1976)  have  shown  that  the  strength  of  polycrystalllns 
BaT103  is  lower  in  the  ferroelectric  state  compared  to  that  in  the 
pereelectric  state  for  grain  sizes  between  1  and  100  vim.  Using  a  fracture 
mechanics  approach,  one  obtains  the  following  expression  for  the  strength  of 
the  material  above  the  Curie  temperature  where  the  material  is  cublcs 

by  iTi 

o  ■  Y  — ~  m 


where  E  is  Young's  modulus,  Yc  is  the  critical  fracture  energy  of  the  cubic 
materiel,  e  is  the  flew  size,  and  Y  is  a  measure  of  the.  flaw  geometry.  For 
the  ferroelectric  state,  this  equation  must  be  modified  to  take  into  account 
both  Internal  stresses  arising  due  to  the  phase  transformation  and  mechanisms 
leaching  to  an  increase  inYJ  1/2 

i  ETf 

0+  <Oj>  -  Y  —  (2) 

where  O^represents  the  average  of  the  tensile  and  compressive  stresses 
around  the  perimeter  of  a  flaw.  The  flew  that  propagates  to  failure  will  be 
the  one  under  the  largest  combination  of  size  and  tensile  stress.  Larger 
flaws  Intersect  many  tensile  and  compressive  components,  leading  to  low  values 
of  <Oj>.  The  contribution  of  <Oj>  to  failure  as  a  function  of  flaw  size 
follows  the  trend  shown  in  Figure  1*  In  order  to  completely  understand  the 
fracture  process,  one  also  must  know  the  effect  of  the  microstructure  of 
BaTi03  on  the  resistance  of  the  material  to  crack  growth  (i.e.,Y). 


In  BaTl.03  the  critical  fracture  energy,  y,  Is  a  function  of  grain 
•ize,  G,  in  the  ferroelectric  state  but  is  independent  of  grain  size  in  the 
paraelectrlc  state  (Pohanka  et  al.  1982)  (Figure  2).  The  150*C  data  fit  an 
empirical  expression: 

Y  -  -  0.02  Log  G  +  3. 46  (3) 

Because  the  crystal  structure  in  the  paraelectrlc  state  is  cubic,  this  type  of 
behavior  is  to  be  expected  (Rice  et  al.  1981).  The  room  temperature  fracture 
energy  of  BaT103,  having  a  tetragonal  crystal  structure,  remains  constant  up 
to  a  grain  size  of  15  ym,  increases  rather  sharply,  and  then  decreases 
(Figure  2).  This  behavior  can  be  explained  on  the  basis  of  a  combination  of 
twinning  and  microcraclting  enhanced  toughness  over  different  portions  of  the 
grain  size  range. 

1  to  5  Urn 

In  the  1  to  5  m  grain  size  regime,  ferroelastic  twinning  does  not  occur; 
therefore,  the  fracture  energy  of  cubic  and  tetragonal  BaT103  would  be 
expected  to  be  quite  similar,  as  is  observed.  These  are  the  lowest  values  of  y 
measured  experimentally. 

5  to  15  yrn 

Above  a  grain  size  of  5  ym,  90-degree,  ferroelastic  twins  ere  produced  by  the 
cubic-to-tetragonal  phase  transformation  (Dennis  1972).  Based  on  the  increase 
in  the  fracture  energy  of  single-crystal  BaTi03  from  0.8  to  1.4  J/m^  due 
to  the  interaction  of  the  crack  with  the  twins,  one  would  expect  the  y  of  the 
polycrystalline  materials  to  be  increased  proportionally.  However,  it  is 
observed  that  fracture  over  this  grain  size  range  is  predominantly 
intergranular.  This  means  that  the  cracks,  Instead  of  passing  through  the 
twin  structure-,  take  the  easier  path  around  the  grains.  Since  no  crack-twin 
interactions  take  place,  twinning  does  not  influence  the  toughness. 


Microcracking  near  the  crack  tip  does  affect  the  measrued  fracture 
energy.  That  microcracking  must  play  some  role  in  the  fracture  process  is 
shown  by  the  spontaneous  cracking  of  polycrystalline  BaTi(>3  at  grain  sizes 
between  150  to  250  ym.  Following  Rice  and  Freiman  (1981),  one  can  calculate 
the  contribution  that  microcracking  would  make  to  the  measured  fracture 
energy;  this  contribution  is  shown  as  the  dashed  line  in  Figure  2. 

Considering  the  number  of  assumptions  and  the  experimental  scatter  in  both  the 
measured  y  and  the  parameters  needed  to  calculate  it,  the  fit  is  quite  good. 

15  to  50  ym 

Over  this  range,  both  microcracking  and  twinning  contribute  to  the  measured 
fracture  energy.  A  maximum  in  the  y-grain  size  curve  at  'xAO  ym  is  calculated 
based  on  microcracking  contributions,  in  good  agreement  with  the  data. 


cc  0 


FLAW  SIZE  (Mm) 


FIGURE  1  Effective  internal  tensile  stress  in  fine-grained  BaTiOj, 

determined  from  measurements  of  flaw  sizes,  fracture  stress  and 
critical  fracture  energy,  as  a  function  of  measured-flow  size. 
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FIGURE  2  Fracture  energy  of  BaTi03  as  a  function  of  grain  size: 

_  best  fit  to  room  temperature  data;  -  contribution  due  to 

twinning,  -  -  -  combined  effect  of  twinning  and  microcracking 
predicted  by  model  of  Rice  and  Freiman  (1981). 


zoemxi*. 
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50  to  200  ua 


At  grain  sizes  greater  than  about  50  ym,  both  twinning  and  aicrocracking  still 
contribute  to  Increased  toughness,  but  linking  of  oicrocracks  begins  to  occur, 
leading  to  decreasing  fracture  energy  with  grain  size.  Eventually,  spontaneous 
failure  takes  place. 


DESIGN  DIAGRAMS  FOR  CAPACITORS ' 


The  above  discussion  relates  the  strength  and  fracture  toughness  of  a  primary 
capacitor  material  to  microstructurally  related  processes.  It  also  is  known, 
however,  that  in  the  combined  pressure  of  stress  and  moisture,  initially  small 
cracks  will  grow  slowly,  leading  to  eventual  failure.  In  this  section  it  will 
be  shown  how  fracture  mechanics  techniques  can  be  used  to  predict  the  time  to 
failure  for  capacitors  containing  cracks  under  stress. 

Cracks  have  been  observed  in  capacitors  and  to  some  degree  can  be 
associated  with  premature  capacitor  failure.  These  cracks  can  be  stressed  by 
a  variety  of  loading  conditions.  These  include  possible  mechanical  loads, 
stresses  due  to  the  differential  thermal  expansion  between  the  capacitor  end 
the  substrate  to  which  it  is  bonded,  and  internal  stresses  due  to  electrical 
field  concentrations  at  the  crack  tip.  Based  on  the  two  equations: 


V  -  A  Kj" 

and 


(4) 


a  -  Y  Kj 


(5) 


where  V  is  crack  velocity,  Kj  is  the  stress  intensity  at  the  crack  tip,  a  is 
the  crack  size,  and  A  and  n  are  parameters  that  are  dependent  on  the  material 
and  the  environment.  A  time  to  failure,  t,  can  be  calculated: 


t 


B  o'\  ("-2) 


where  B  is  a  constant  involving  A  and  Y  and  is  the  initial  strength  of 
the  material  (i.e.,  a  measure  of  the  crack  size  distribution  before  any 
stresses  are  applied).  The  constants  n  and  B  can  be  determined  by  measuring 
the  strength  of  the  capacitors  as  a  function  of  stressing  rate.  The  tests 
usually  are  performed  in  four-point  flexure  over  three  to  four  orders  of 
magnitude  in  stressing  rate  as  shown  in  Figure  3  for  a  typical  capacitor 
composition.  Approximately  100  specimens  are  required  to  obtain  accurate 
measures  of  n  and  B.  The  initial  flaw  size  distribution  is  determined  by 


*Detalls  of  the  analysis  techniques  for  tests  of  this  type  are  given  by  Jakus 
and  co-workers  (1978)  and  Ritter  and  co-workers  (1981). 
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FIGURE  5  Comparison  of  failure  time  predictions  for  two  capacitor  materials. 
Ceramic  1;  n  ■  19.5  and  InB  •  30.95*  Ceramic  2;  n  ■  42.2  and 
InB  —  27.0. 


measuring  the  strength  under  inert  conditions  (e.g.,  dry  M2  8ns)  at  a 
relatively  rapid  rate  of  loading.  The  distribution  in  Inert  strengths  Is 
analyzed  using  Welbull  statistics.  Weibull  theory  allows  us  to  assign  a 
probability  of  failure  for  a  given  level  of  stress.  Combining  these 
statistics  for  with  the  crack  growth  parameters  obtained  from  the 
stressing  rate  tests  allows  a  design  diagram  such  as  that  shown  in  Figure  4 
to  be  constructed.  The  three  sets  of  curves  in  this  diagram  represent  the 
same  material  with  three  different  initial  flaw  sizes.  The  central  line  in 
each  set  is  the  calculated  values  of  failure  time  for  a  given  applied  stress 
in  Equation  6.  The  bounding  lines  represent  the  95  percent  confidence  limits 
on  the  failure  times  and  are  based  on  the  errors  in  determining  B  and  n*  One 
can  see  that  errors  in  B  and  n  lead  to  large  uncertainties  in  estimating  the 
safe  lifetime  for  a  given  applied  stress.  One  can  reduce  this  uncertainty, 
however,  by  reducing  the  applied  stress  to  the  edge  of  the  confidence  band. 

Figure  5  shows  the  failure  time  data  for  two  different  capacitor 
compositions  plotted  in  a  way  that  normalizes  the  flaw  sizes.  It  can  be  seen 
that  the  material  having  the  larger  value  of  n  would  be  expected  to  survive 
for  much  longer  times  over  most  of  the  range  of  applied  stresses. 
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SUMMARY  AMD  CONCLUSIONS 


Based  on  the  above  discussions,  it  seems  clear  that  the  application  of 
fracture  mechanics  principles  can  contribute  significantly  to  an  understanding 
of  the  fracture  processes  in  capacitor  type  materials.  Both  the  chemical 
composition  of  the  material  as  well  as  the  grain  size  affect  the  material's 
resistance  to  both  subcritical  crack  growth  and  catastrosphic  failures.  By 
measuring  the  slow  crack  growth  parameters  for  a  material,  it  is  possible  to 
predict  lifetimes  under  stress. 

There  are  a  number  of  areas  in  need  of  further  work  however.  These 
include: 


1.  The  quantitative  dependence  of  crack  growth  parameters  on 
composition  and  microstructure, 

2.  The  effects  of  internal  stresses  and  flaw  and  grain  sizes  on  crack 
growth, 

3.  Crack  growth  in  complex  stress  fields,  and 

4.  More  accurate  techniques  for  analysis  of  static  fatigue  and 
stressing  rate  data. 
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